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The absorption in lead of the shower rays which produce 


the bursts of cosmic-ray ionization is measured by two 


methods. The first method consists in observing the 


ionization produced above and below a lead absorber 
placed across the center of an ionization chamber; the 
second is to observe the probability that a burst of ioniza- 
tion in a chamber is accompanied by a simultaneous 
discharge of three Geiger-Miiller counters over one of 


which has been placed an absorber. The results of the two 


phage the first experiments of Rossi,' in 1932, 
\7 on the measurement by means of Geiger- 


Miiller counters, of the absorption of cosmic-ray 


showers, many refinements have been made, 
particularly by H. Geiger? and his collaborators 
and by J. C. Street,’ which have led to a better 
understanding of the phenomena involved. How- 
ever, since these showers occur with a large 
range of sizes and degrees of complexity, addi- 
tional information can be obtained by observing 
the behavior of showers of definite size, as may 
be done with an ionization chamber, instead of 
observing only the integrated effect of showers 
of all sizes, as is done in counter experiments. 
The experiments here described were begun in 
collaboration with Professor W. F. G. Swann, 
and some preliminary reports have been pub- 
lished. The absorption in lead of large showers 
(of the order of a hundred rays, or more) has 


1B. Rossi, Zeits. f. Physik 82, 151 (1933). 

2H. Geiger and O. Zeiller, Zeits. f. Physik 97, 300 (1935) 
and earlier papers. 

3 R. H. Woodward and J. C. 
(1935) and earlier papers. 


Street, Phys. Rev. 47, 800 


methods are in good accord and may be stated in the form 
that the probability that a ray of a shower will penetrate 


a thickness of lead decreases linearly with the thickness, 


becoming zero at approximately 11 cm. The experiments 
serve to emphasize again the high energies that are 


involved in a large shower. The results are applied to 


observations on the effect of shielding on the ionization 


observed in the stratosphere 


been measured by two methods in which the 
ionization that a shower produces is measured 


before the shower is absorbed. 


First METHOD 


The principle of this ionization method is to 
measure the ionization which a given shower 
produces before and after passing through a 
thickness of material. Since, presumably, the 
ionization which a shower produces is propor- 
tional to the number of rays which it contains, 
the observed difference is a measure of the num- 
ber of rays which are stopped within the ma- 
terial. To measure this, a large cylindrical 
ionization chamber, 150 cm high and 90 cm in 
diameter, used. The chamber was 
structed of welded steel of approximately 18 mm 
thickness. It was divided into two halves, see 


was con- 


4\W.F. G. Swann and C. G. Montgomery, Phys. Rev. 
43, 782 (1933); 44, 52 (1933). W. F. G. Swann, ‘‘Report on 
the Work of the Bartol Research Foundation,”’ 1933-1934, 
J. Frank. Inst. 218, 173 (1934). 
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Fic. 1. Diagram of double ionization chamber with ab- 
sorber in place, showing the electrode arrangement. 


Fig. 1, each of 538 liters volume, and each 
provided with an independent electrode system 
and guard ring. The electrodes were concentric 
sheet iron cylinders with conical caps, separated 
by approximately 10 cm. To collect the ions, 
the electrodes were maintained at a potential of 
about 500 volts, and the time of collection of the 
ions was less than one second. These conditions 
insure that the statistical fluctuations of the 
cosmic-ray currents are small.' The potential of 
each electrode system was recorded photo- 
graphically by allowing the image of an illumi- 
nated slit reflected from a galvanometer mirror 
to fall upon a moving strip of photographic 
paper. Each galvanometer measured the changes 
in the plate current of an FP-54 pliotron in the 
usual manner. The chamber contained nitrogen 
at 6.8 atmospheres pressure. The capacities of 
the electrodes were determined by inserting 
known resistances in the circuits and measuring 
the time constants of the decay of potential 
differences between the electrodes. The capacity 
of each system was found to be about 500 cm. 
The thick steel walls of the chamber and the 
surrounding material served as a source of 
showers. The ionization which these showers 
produced was measured in each half of the 


5C. G. Montgomery and D. D. Montgomery, Phys. 
Rev. 47, 430 (1935). 
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chamber and the frequency distribution of the 
shower sizes obtained. Then, to absorb the 
showers, a circular plate of lead, 7 cm thick and 
80 cm in diameter, was inserted between the two 
halves of the chamber, and the frequency 
distribution of the sizes of the showers in each 
half redetermined. Only those showers which 
produced more than 410° ions in each half of 
the chamber were measured. This amount of 
ionization would be produced by the passage of 
about 150 high speed electrons through the 
chamber, if the value for the specific ionization 
is taken as 60 ions per centimeter.® The observa- 
tions were extended over two 35-hour periods, 
during each of which more than 200 showers were 
observed. Since the root-mean-square fluctuation 
in the cosmic-ray ionization in the time necessary 
for the galvanometer to deflect was only 6X 10° 
ions, the showers of the sizes measured were 
quite unaffected by these fluctuations. It is of 
interest to note that the fluctuations in the two 
halves of the vessel are not statistically inde- 
pendent, but show a correlation caused by the 
passage of cosmic-ray electrons through both 
halves. The correlation coefficient is approxi- 
mately 0.2.7 

The interpretation of shower observations is 
complicated by the fact that the showers which 
emerge from a piece of material are not all of 
equal size, but are of many sizes whose frequency 
distribution covers a wide range. In the present 
experiment, there is the further complication 
that showers of equal size produce different 
amounts of ionization, the amounts depending 
upon their points of origin and their path lengths 
in the chamber. Thus, we observe three kinds of 
showers: those which produce amounts of ioni- 
zation above the limit of measurement in the 
lower half only, those whose ionizations are 
above the limit in the top half only, and those 
which produce measurable amounts of ionization 
in both halves. Table I gives the numbers of 


TABLE I. Numbers of bursts observed in 35 hours. 





Botu UPPER LOWER 
Hatves Harr ONLY HALF ONLY Toral 
Without lead 49 89 78 216 


With lead 24 142 102 268 





®\W.F. G. Swann, Phys. Rev. 44, 961 (1933). 
7Cf. W. S. Pforte, Zeits. f. Physik 72, 511 (1931). 
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Fig. 2. Frequency distribution of bursts which occurred 
in the lower half only. 


bursts of these three classes observed in 35 hours, 
with and without a lead absorber. 

The effect of introducing the lead is twofold. 
As well as acting as an absorber of showers 
produced above it, the lead acts as an additional 
source of showers. The showers produced in it 
show themselves most prominently, of course, 
in the lower half of the chamber. Thus, if we 
consider the showers which occur in the lower 
half of the vessel only, we see that more showers 
are observed with the lead present than without 
the lead. Fig. 2 shows the distribution in size of 
the bursts of ionization which occur in the lower 
half only. However, the difference between the 
two curves does not represent all the showers 
which the lead has produced. Certain showers 
whose points of origin were in the upper half of 
the chamber but whose sizes, in the absence of 
the lead, were above the limit of measurement in 
the lower half only, are now absorbed by the 
lead and become too small in the lower half to be 
measured. Therefore, the number of these 
showers should be added to the difference be- 
tween the two curves in Fig. 2 to obtain the 
total number of showers produced by the lead. 

The phenomenon in which we are at present 
most interested is the absorbing effect of the 
lead upon the showers which pass through it. 
We see from Table I that the number of showers 
which produce more than 4X 10° ions in both the 


upper and lower halves of the vessel are con- 
siderably reduced by the presence of the lead. 
This reduction in number is to be regarded as 
the effect of the lead in reducing the size of the 
showers and making some of them fall below the 
limit of measurement in the lower half. We may 
obtain an estimate of the amount of this reduc- 
tion in size in the following way. Fig. 3 shows the 
frequency distributions of the ratios of the 
amounts of ionization produced in the lower half 
of the vessel to those produced in the upper half 
by the showers under consideration, with and 
without the lead absorber. When no lead is 
present, most of the showers occurring in the 
upper half of the chamber produce, in the lower 
half, an amount of ionization only slightly less 
than in the upper half. This slight difference is 
to be ascribed to the spreading of the shower 
rays from their point of origin and their passing 
out through the side walls of the vessel. When 
the lead is present, the maximum in the distribu- 
tion of ratios of sizes shifts to smaller values, 
since the sizes of the showers in the lower half 
of the vessel have been reduced without altera- 
tion of the sizes in the upper half. We may take 
the ratio of the positions of the maxima of the 
two curves as an estimate of the most probable 
value of the fraction of the rays of a shower 
which are able to pass through the lead.* This 
value is 0.35/0.85, or 0.4. 
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Fic. 3. Frequency distribution of bursts which occurred 
in both halves of the vessel, with respect to the ratio of 
the size in the lower half to that in the upper. 


8 It is conceivable that there are showers produced in the 
lead which proceed upward as well as downward. How- 
ever, the agreement between the two methods of estimating 
the absorption in lead indicate that such showers, if 
present, are relatively few in number. 
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Fic. 4. Frequency distributions of bursts which occurred 
in both halves of the chamber. 


Having obtained a numerical estimate of the 
reduction which the lead brings about in the size 
of a shower, we may use it to calculate the 
number of showers which we should expect to be 
reduced to a size below the limit of measurement 
in the lower half of the vessel. This expected 
reduction, by the lead, of the number of showers 
which are measurable in both halves of the 
chamber, may then be compared to the observed 
reduction. Fig. 4 shows the frequency distribu- 
tion of the sizes of showers which are above the 
limit of measurement in both halves of the vessel, 
with and without the lead present. The lead 
removes all showers which would have produced 
less than 4/0.4 10°, or 10 10° ions in the lower 
half of the chamber. Let us take this number of 
showers to be equal to the number below 10 X 10° 
ions which did occur when no lead was present. 
This number is 37. The observed decrease is 25. 
These showers, the bottom components of which 
have been absorbed, are still observable in the 
top half of the chamber, and thus the introduc- 
tion of the lead should cause a corresponding 
increase in the number of showers observed in 
the upper half only. Fig. 5 shows that there is 


Ange BD. DB. 
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such an increase and within the limits of the 
experimental error it is the number expected. 
Thus the observations may all be consistently 
explained if the absorbing effect of the 7 cm of 
lead is that it reduces the number of rays in a 
shower passing through it by a factor of 1/0.4. 
This estimate of the absorbing power of a lead 
plate may be compared with the absorption 
effects as measured by a second method, de- 


scribed below. 


SECOND METHOD 


A second method of measuring the absorption 
in lead of a cosmic-ray shower is based upon 
measurements of the probability that a shower 
of a given size will simultaneously discharge 
Geiger-Miiller counters. In previous 
it has been shown that, with reasonable 


is possible to compute and 


several 
papers,” 
assumptions, it 
measure this probability. Suppose now that an 
absorber is interposed between the source of 
showers and one of the Geiger-Miiller counters. 
The probability that a ray of the shower will 
set off that counter will then be decreased, and 
the probability that a ray of the shower can 
pass through the absorber can be derived. Fig. 6 
shows the experimental arrangement used. The 
source of showers consisted of a tray, L, 41.5 cm 
square, containing lead shot to a depth equiva- 
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Fic. 5. Frequency distribution of bursts which occurred 
in the upper half only. 


°C. G. Montgomery and D. D. Montgomery, Phys. 
Rev. 48, 786 (1935); J. Frank. Inst. 221, 59 (1936). 
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Fic. 6. Experimental arrangement used in the 
second method. 


lent to one centimeter of solid lead. The showers 
were detected by means of the magnesium 
ionization chamber, S, and the vacuum tube 
electrometer, E. The shower rays passing out of 
the chamber fell upon three Geiger-Miiller 
Over the 


” 


counters, C, arranged in a ‘“‘cradle. 
center counter was placed an absorber, A, of 
solid lead. A simultaneous discharge of the three 
counters made an imprint upon the same photo- 
graphic paper that recorded the bursts of ioni- 
zation caused by the shower. The fraction of the 
number of bursts of ionization which were 
accompanied by a discharge of the set of counters 
was determined for different thicknesses of the 
absorber, A. Table II gives the total numbers of 
bursts of ionization observed in the chamber and 
the number of bursts with which a simultaneous 
discharge of the counters occurred. The bursts 
of ionization are divided, according to size, into 
five groups. During the course of these observa- 
tions, the high resistances associated with the 
counters had two different values. The series of 
observations labeled I were taken with a lower 


value of these resistances than the series labeled 
II. The counter efficiencies in the series I were 
therefore higher than in II, but the data were 
treated in such a way, as explained below, that 
the efficiencies of the counters entered only as a 
small correction, and the final results are strictly 
comparable. 


The efficiencies of the counters were deter- 
mined in the usual manner. The counters were 
arranged so that they formed a telescope pointing 
in the vertical direction. The number of triple 
coincidences and the individual counting rates 
of each counter were measured with and without 
radium in the vicinity. From these data the 
recovery time of a counter was computed and 
hence the efficiency for any given value of the 
individual counting rate. The efficiency for a 
counter with no lead over it was 95 percent for 
the observations in series I, 85 percent in series 
II. 

In order to derive from the above observations 
a value for the probability that a ray of a shower 
will penetrate the absorber, we must apply the 
formulae previously calculated. If N is the 
number of rays in a shower, p the a priori 
probability that a ray of the shower will pass 
through a counter in the absence of the absorber, 
p’ the probability that a ray will penetrate the 
absorber, and Ey and E, the efficiencies of a 
counter without and with an absorber over it 
respectively, then it can be shown that the 
probability that the occurrence of a shower will 
be accompanied by a simultaneous discharge of 
the three counters is: 

P=EE.[1—2(1-—p)* —(1—pp’)*¥+(1—2p)4 
+2(1—p—pp’)* —(1—2p—pp’)* ]. 
For values of N and p which are here under 
TABLE II. Bursts of various sizes observed for different 
thicknesses of lead absorber. Rows labeled A give the 


number of bursts accompanied by counts. Rows 
labeled B give the total number of bursts. 


Size or Burs 1.1 to 20to 30to 5.0to10.0to 
x10 -*ion pairs 0 3.0 5.0 10.0 20.0 
A 117 38 18 16 2 
No Pb 
B 755 105 36 22 3 
Series A 92 27 16 17 5 
I 2.54 cm Pb 
B 632 97 44 28 5 
l 54 20 17 10 4 


5.08 cm Pb 
533 102 45 15 4 


1 62 30 24 12 8 
No Pb 
B 509 101 35 15 9 
Series A 32 18 12 4 3 
II 7.62 cm Pb 


B 412 92 46 15 6 
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penetrate a given thickness of lead. A, mean value un- 

corrected for efficiency; B, observations by first method. 

@ weighted mean values; A, 1.1 to 2.0X10® ion pairs; 

Vv, 2.0 to 3.0 10° ion pairs; X, 3.0 to 5.0X 10° ion pairs; 
], 5.0 to 10.0 x 10° ion pairs. 


consideration, it is sufficient to employ the 


approximate expression : 
P=E,?E,(1—e-*%”)*(1—e-*? ). (1) 
If we designate by P, and P» the observed 
fraction of the bursts of ionization which were 
accompanied by a simultaneous discharge of the 
three counters with and without the absorber, 
respectively, we can solve Eq. (1) for p’ thus: 


p’ =log (1—P/EaPo')/log (1—Po'/E,). 


We notice that the value of p’ so derived is 
independent of N and hence independent of the 
value of the specific ionization of a shower ray. 
The values of p’ were calculated in this way for 
the first four groups of shower sizes in Table II, 
and are plotted in Fig. 7. The values of p’ that 
we calculate are, of course, relative to the value 
of p’ for the condition of no absorber, which 
value has been taken as unity. Now, even when 
no absorber is present, a ray of the shower, to 
discharge a counter, must pass out through the 
wall of the ionization chamber and through the 
wall and shield of the counter. The amount of 
material which it must traverse is 2.9 g/cm, 
which is equivalent to 2.5 mm of lead. Hence 
the abscissae in Fig. 7 have been displaced by 
this amount. The probable error for each value 
of p’ was estimated in the usual way from the 
number of bursts observed and accordingly the 
weighted mean value of p’ for each thickness of 
lead was computed. These mean values with 
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their probable errors are also recorded in Fig. 7. 
The values of ~’ so obtained are largely inde- 
pendent of the efficiency of the counters. The 
point labeled A in Fig. 7 is the value computed 
assuming the efficiency of the counters to be 
100 percent. It differs very little from the value 
computed when the actual efficiencies were taken 
into account, although the error introduced by 
assuming 100 percent efficiency would be greatest 
for this particular thickness of lead. 

We see that, within the limits of the experi- 
mental error, p’ decreases linearly with the 
thickness of the absorber and is independent of 
the size of the shower. The point labeled B in 
Fig. 7 represents the value of p’ obtained by the 
ionization method described in the first section. 
The agreement between the two methods is thus 
very satisfactory. If we extrapolate the observa- 
tions to obtain the maximum range of a shower 
in lead, we find a value of 11 cm. 


DISCUSSION 


In the counter experiments on the absorption 
of showers, the counts observed* 
decreases very sharply with the introduction of 
a lead absorber a few millimeters thick, and then 
as the thickness of the absorber is increased, 
much more 


number of 


the number of counts decreases 
gradually. Of course, it is not possible to express 
the results of the counter observations directly 
in terms of the probability that a shower ray 
will penetrate the absorber, since the frequency 
distribution in size of the showers counted is not 
known, and hence the results of the present 
experiment may not be directly compared with 
the counter Unfortunately the 
ionization chamber used for our experiments is 
so thick that the sharp initial decrease in the 
probability of penetration, if it were present, 
could not be observed. The measured proba- 
bilities of a shower ray penetrating a given 
thickness of lead, as given in Fig. 7, are relative 
values calculated by choosing as unity the 
probability when no additional absorber other 
than the chamber itself was present. However, 
the very gradual falling off of the number of 
showers, observed in counter experiments, as 
large thicknesses of lead are added is in good 
accord with the experiments here described. 


observations. 
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These experiments emphasize again the fact 
that very large amounts of energy must be 
involved in the production of showers. Using 
Anderson's!” value of the energy loss of electrons 
in lead, we see that at least half of the rays of 
the shower which emerge from the chamber have 
energies greater than 310° electron volts, and 
that the total energy in the showers which we 
observe must often exceed 3 X 10" electron volts. 
It is also apparent that such penetrating powers 
as we observe are incompatible with the simple 
identification of the position of the maximum of 
a Rossi curve with the range of the shower 
particles. 

An interesting consequence of these experi- 
ments is the quantitative agreement with the 
observations on the ionization produced in 
shielded and unshielded chambers in the strato- 
sphere. As we go to higher and higher elevations, 
the contribution to the ionization of the ex- 
tremely penetrating (primary) cosmic rays be- 
comes smaller and smaller relative to the 
contribution of the secondary (shower) particles 
produced in the atmosphere. Thus we should 


10C, D. Anderson and S. H. Neddermeyer, Papers and 
Discussions, International Conference on Physics, London, 


1934, I, p. 171. 
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expect that the percentage decrease in the 
ionization brought about by shielding the 
chamber should approach the value observed in 
the present experiments for the probability that 
a shower ray will be stopped in the shield. 
Bowen, Millikan, and Neher" have published a 
curve showing the variation with elevation of 
the decrease in ionization caused by a 6.5-cm 
lead shield, derived from their measurements and 
those of Compton and Stevenson."? The per- 
centage decrease approaches the value of 68 
percent as the elevation increases. The proba- 
bility that a shower ray will be stopped in this 
shield thickness is 62 percent. Thus the observa- 
tions are in accord with the interpretation 
suggested here. 

The authors wish to express their gratitude to 
Professor W. F. G. Swann, in collaboration with 
whom a portion of these experiments were 
carried out, and whose interest and guidance 
throughout the course of this investigation have 
been invaluable. 


1]. S. Bowen, R. A. Millikan and H. V. Neher, Phys. 
Rev. 46, 646 (1934). 

2 A. H. Compton and R. J. Stevenson, Phys. Rev. 45, 
441 (1934). 
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Coincidence Counter Studies of Cosmic-Ray Showers 


R. H. Wooowarp, [Iarvard University 


(Received March 28, 1936) 


It has been found that the production curves (counting 
rate against thickness of lead) of electron-produced and of 
photon-produced showers are similar and that the average 
penetrating power of the rays from electron-produced and 
from photon-produced showers is the same. Shower 
production curves have been obtained at four elevations 
and analysis of the curves shows that the absorption per 
nucleus of the shower-producing radiation (photons) is 


1. INTRODUCTION 


HE cloud chamber photographs of Anderson 
et al.' and Blackett and Occhialini? have 
shown that cosmic-ray showers consist of groups 


1C. D. Anderson, R. A. Millikan, S. Neddermeyer and 
W. Pickering, Phys. Rev. 45, 352 (1934). 

?P. M.S. Blackett and G. P. S. Occhialini, Proc. Roy. 
Soc. A139, 699 (1933). 


approximately proportional to the square of the atomic 
number. Coefficients for lead (0.33 cm"), iron, and air 
have been determined. and a comparison with theory is 
given. Also the penetrating power of the shower rays 
emerging from a block of lead has been measured and 
found to be independent of the elevation and of the 
thickness of the lead from which the rays emerge. 


of from two to several score electrons which are 
accompanied by numerous low energy photons. 
The shower-producing rays are usually non- 
ionizing, although the photographs of Stevenson 
and Street* show that ionizing rays may produce 


3E. C, Stevenson and J. C. Street, Phys. Rev. 48, 464 
(1935). 
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showers (possibly through the intermediary of 
nonionizing rays). Many of the shower electrons 
radiate from a single point. There is a high 
probability, however, that the production of one 
shower will be accompanied by the simultaneous 
production of other showers, the several shower 
foci usually being unconnected by tracks of 
ionizing rays. Either in the production of a 
shower or in a primary process nonionizing rays, 
photons, are created which themselves are capa- 
ble of producing showers. 

Geiger and Fiinfer* have proposed a convenient 
shall adopt in this 


nomenclature which we 


discussion : 


A-rays. The primary charged particles which lose energy 
through ionization and the radiation of B-rays (photons) 
showers 


in nuclear collisions. The latitude effect of 


requires that they be related to a primary corpuscular 
radiation. 

B-rays. The shower producing photons. 

C-rays. Shower electrons created in groups by B-rays. 

D-rays. Low energy photons produced in the absorption 
of C-rays. 

E-rays. Low energy electrons from photo and Compton 
collisions of the D-rays. The D- and E-rays account for 
the back scattering reported by Fiinfer.? The distinction 
between A-, C-, and E-rays is not clearly defined. Their 
properties may be similar, and there may be overlapping 
of energies of different groups. Likewise the B- and D-rays 
may behave in similar fashions. But in general the classes 
are arranged in order of decreasing energy. 


the coincidence 


counter method has been made by Rossi* and 


The study of showers by 


many others, and the results of these experiments 
have been discussed in detail by Geiger.’ In 
general, curves of the type shown in Fig. 2 are 
obtained. Such a curve relating the counting 


rate N to the thickness XY of material can be 
roughly expressed as the difference of two 
exponentials, 

N=C(e* —e~"*). (1) 


Here » may be considered as the absorption 
coefficient of the primary radiation and » that of 
the secondary radiation or vice versa, and conse- 
quently there are two possible interpretations of 


4H. Geiger and E. Fiinfer, Zeits. f. Physik 93, 543 
(1935). 

5 T. H. Johnson, Phys. Rev. 47, 318 (1935). 

°T. H. Johnson and D. N. Read, Abstract 2, New York 
Meeting of Am. Phys. Soc., Feb., 1936. 

7 E. Fiinfer, Zeits. f. Physik 83, 92 (1933). 

* B. Rossi, Zeits. f. Physik 82, 151 (1933). 

*H. Geiger, Ergeb. d. Exakt. Naturwiss. 14, 42 (1935), 
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the curve. The initial rise may represent the 
coming into equilibrium of a secondary radiation 
(v~0.7 per cm of lead) with a more penetrating 
shower-producing radiation (u~0.3 cm~'), the 
absorption of which accounts for the fall of the 
curve beyond the maximum. The alternative 
interpretation is that the shower-producing radia- 
tion is rapidly absorbed (v~0.7 cm™), while the 
secondary radiation is the more penetrating 
(u~0.3 cm—'!). Direct measurements (see sections 
2 and 4) of the absorption of the secondary rays 
emerging from a block of lead show that most of 
these rays are not penetrating, being absorbed 
to a large extent in two millimeters of lead. 
Furthermore it does not seem plausible that the 
secondary radiation, which on the average must 
be considerably less energetic, should be more 
penetrating than the primary radiation. For 
these reasons we prefer the former interpretation 
and in our discussion shall disregard the possi- 
bility of the latter. 

2. A COMPARATIVE STUDY OF ELECTRON AND 
PHOTON-PRODUCED SHOWERS 


The production of cosmic-ray showers has been 
studied by the arrangement shown in Fig. 1. The 
counters in each of four groups were connected in 
parallel so that each group acts as a single 
counter. Iron blocks were placed below and at 
the sides of the lower counters to reduce the 
number of cross-scattered rays. The triple coinci- 
dence counting rates Nj23; between groups 1, 2, 
and 3, and Ne3s between groups 2, 3, and 4 were 
recorded simultaneously for various thicknesses 
of the lead in position A. A coincident discharge 
of groups 1, 2, and 3 requires ionization above the 
lead, whereas a discharge of groups 2, 3, and 4 is 
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Fic. 1. Arrangement of counters. 
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Fic. 2. Curves of shower production in lead. 


independent of the nature of the radiation above 
the lead. As shown in Fig. 2, the variation of the 
counting rate Noe3. with thickness of lead is 
similar to that of Nj23, the factor of five between 
the scales presumably depending on the geometry 
of the arrangement. Stevenson and Street! have 
shown that about 70 percent of the coincidences 
of type Ni23 due to showers produced in the lead 
actually arise from incident electrons. Most of 
the coincidences of type Nes, involving no selec- 
tion in favor of incident electrons are due to an 
incident photon radiation.! It is surprising, 
therefore, that the curves are so nearly similar, 
for in each case the shape of the curve and the 
position of the maximum are dependent upon the 
relative absorptions of the incident and secondary 
radiations. However, there is evidently a real 
difference in the curves at great thicknesses of 
lead. The significance of this difference will be 
discussed in section 3. 

The absorption of the C-, D-, and E-rays which 
actually discharge the lower counters has been 
measured by placing plates of absorbing material 
of various thicknesses in position B, so that they 
absorb the rays produced in the block A but do 
not themselves produce many showers capable of 
discharging the triple sets of counters. Fig. 3 
shows a schematic diagram of the arrangement 
together with a plot of four pairs of absorption 
curves which were obtained by placing various 
thicknesses of lead or aluminum in positions B 
under thicknesses 0, 0.476, 1.59, and 7.94 cm of 
lead at A. The thickness of aluminum is plotted 
on such a scale that equal distances along the 


10 E. C. Stevenson and J. C. Street, Phys. Rev. 49, 425 
(1936). 


abscissa represent equal numbers of extranuclear 
electrons per cm? in the lead and aluminum. The 
nature of the absorption curves is nearly similar 
for the two types (Vi23 and Ne34) of coincidences 
discussed above, and only one set of curves is 
given. This indicates that there is no distinction 
between the photon and_ electron-produced 
showers as regards to penetrating power of the 
secondaries. The interpretation of the curves is 
complicated by the fact that the lead and 
aluminum in positions B are not pure absorbers ; 
i.e., there are coincidences due to the production of 
showers and also there is possibly an increase of 
the efficiency of the photons (D-rays) in dis- 
charging the lower counters. These effects account 
for the peaks in the absorption curves, and 
corrections for them are unreliable. However, the 
initial rapid fall of the curves, which probably is 
not greatly influenced by these effects, shows 
that most of the shower rays are absorbed in two 
millimeters of lead. Also the curves show that 
the penetrating power of the shower rays does 
not change appreciably as the thickness of the 
lead block A is increased. The absorption per 
extranuclear electron is less rapid for aluminum 
than for lead, indicating that the absorption of 
shower rays depends not only on extranuclear 
electrons but also on the nuclei. 


3. THE VARIATION OF THE RATE OF SHOWER 
PRODUCTION WITH ELEVATION 

The coincidence rate with no lead in positions 

A and Bisaconsiderable fraction of the maximum 

change due to the introduction of lead, and it is 
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Fic 3. Curves of absorption in lead and aluminum of 
shower rays from lead. The thickness of aluminum is 
expressed in equivalent cm of lead on a scale of equal 


numbers of electrons per cm*. 
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lic. 4. Arrangement of counters 
not clear how corrections can be made for the 
variation of this ‘‘background’’" rate with the 
introduction of lead at A and B. This difficulty, 
together with the fact that lead at B does not 
behave as a perfect absorber renders a quanti- 
tative treatment of the data obtained with the 
arrangement of Fig. 1 uncertain. The arrange- 
ment of counters shown in Fig. 4 was designed to 
minimize these difficulties. Quadruple coinci- 
dences between the four pairs of counters were 
recorded for various thicknesses of lead placed 
above the counters at A. The zero reading is too 
low to distort the measurements and the quad- 
ruple coincidence counting rate due to showers 
produced in the absorbing plates at B is negli- 
gible; i.e., no distortion of the absorption curves 
is observed. With this arrangement the previous 
measurements have been repeated at Cambridge 
and at three higher elevations. The results with 
the probable errors are tabulated in Table I and 
the principal curves are plotted in Fig. 5. All 
measurements were taken under canvas, and 
corrections were made for barometric fluctuations 
and for change of the efficiencies” of the counters 
with altitude and with thickness of lead. The 
accidental counting rate was negligible (order of 
0.0001 per min.). The four curves are similar 
within the limits of experimental error. The 
increase of the rate of shower production with 
elevation is in satisfactory accord with other 


1! Supplementary tests have shown that not more than 
ten percent of this rate is due to accidental coincidences. 
The major part is no doubt due to showers from the air or 
light wooden roof above the counter set. 

27. C. Street and R. H. Woodward, Phys. Rev. 46, 
1029 (1934). 


WOODWARD 


coincidence counter observations.*: '*» '* However 
it is in serious disagreement with the ionization 
chamber and D. D. 
Montgomery." They find that the frequency of 
cosmic-ray bursts increases by a factor of 26.6 
from sea level to Pike’s Peak (4300 m), whereas 
our coincidence rates increase by a factor of only 
8.5 from sea level to Mt. Evans (4300 m). As a 


observations of C. G. 


possible explanation for the discrepancy, we 
suggest that the small showers (an average of 3 
or 4 rays per shower) recorded by the counter set 
increase with altitude less rapidly than the large 
bursts (estimated at 100 or more rays per burst) 
recorded by their ionization chamber. The results 
obtained at the same stations (Cambridge, 
Denver, Echo Lake, and Mt. Evans) by R. T. 
Young'® with a small ionization chamber seem 
to support this explanation. He finds that while 
the small bursts (10 to 15 rays) increase by a 
factor of 9.0 between sea level and Mt. Evans, 
the large bursts (greater than 30 rays) increased 
by a factor of 22. 

Gilbert!’ has reported that the peak of the 
shower production curve shifted from 1.6 cm of 
lead at sea level to 2.2 cm at an elevation of 
3500 m. Although there is an indication of such a 
shift of the peak in our curves, it is certainly 
smaller, if it exists at all. Here again the expla- 
nation may lie in the different types of showers 
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Fic. 5. Shower production curves at Mt. Evans (4300 
meters), Echo Lake (3250 m), Denver (1620 m), and 
Cambridge (60 m). 


13B. Rossi and S. de Benedetti, Ricerca Scient. 5-I], 
379 (1934). 

4]. A. Priebsch, Zeits. f. Physik 95, 102 (1935). 

1 C, G. Montgomery and D. D. Montgomery, Phys. 
Rev. 47, 429 (1935). 

1 R. T. Young, Abstract 10, New York Meeting Am. 
Phys. Soc., Feb. 1936. 

17 C, W. Gilbert, Proc. Roy. Soc. Al44, 559 (1934). 
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raBLeE I. Quadruple coincidence counts for various thicknesses of lead placed above the counters at A and for various altitudes. 


? 


Column 1 gives the position A in cm of lead; column 2, 


the position of B in cm of lead; columns 3, 6, 9, and 12 


give the coincidence counts per minute, NV; columns 4, 7, 10, and 13 give the probable error as defined by 
0.67[=(AN)*?/n(n—1)}*; columns 5, 8, 11, and 14 give the probable error as defined by 0.67¢1/T. 


BAROMETRIC 





PRESSURE 76 cm 64 cm 52cm 44.om 
ELEVATION 60 meters 1620 meters $250 meters 4300 meters 
1 4 3 4 3 6 8 9 10 1 1 12 13 14 
0.00 0.00 0.014 0.0017 0.0010 0.066 0.007 0.004 0.138 0.005 0.006 0.206 0.012 0.008 
.63 .00 39 004 012 1.28 017 033 
95 .00 1.23 007 .005 
1.27 .00 320 .0062 .0048 75 019 014 2.82 037 035 
1.59 .00 388 0064 .0054 .89 017 012 1.84 032 .028 3.35 .020 040 
1.91 .00 .398 0072 .0051 91 O11 012 1.98 .040 025 3.38 .034 027 
2.22 .00 398 .0065 .0064 85 O11 012 2.00 .026 025 3.44 038 030 
2.54 .00 345 .0075 .0061 1.97 019 .027 3.35 .040 044 
2.85 .00 42 018 013 
3.18 .00 1.59 017 .034 2.86 .054 .037 
4.77 .00 197 0037 .0047 43 014 012 0.93 .007 028 1.84 .039 .030 
7.30 .00 1.02 017 034 
9.85 .00 102 0021 0036 20 .006 007 46 030 O19 0.79 009 O18 
12.40 .00 77 007 030 
0.00 63 043 006 
00 5.08 O11 0017 030 006 071 009 105 008 
63 Be | 62 024 
63 1.59 39 021 
63 5.08 19 010 
1.91 16 .259 .0057 .60 O15 
1.91 32 161 .0042 37 012 85 028 1.54 042 
1.91 .63 115 .0047 9 008 .60 029 
1.91 1.11 .069 .0042 
1.91 1.59 .063 .0033 14 .007 39 .024 .63 041 
1.91 5.08 013 .0015 .035 .004 17 .010 .29 .016 


recorded. Gilbert’s arrangement of counters re- 
quired three rays to produce a_ coincidence 
whereas ours required only two. It is evident 
from a comparison of the results of sections 2 and 
3 that there is a shift of the peak due to the use of 
different geometrical arrangements. Probably 
this shift is not significant, for the distortion 
caused by the variation of the initial counting 
rate with the introduction of lead may well 
explain the entire difference. It is to be expected 
that the superposition of an initial decreasing 
curve upon the production curve will shift the 
peak to a smaller thickness. 

The intensity of the single corpuscular radia- 
tion was measured by recording double coinci- 
dences between one of the upper pairs of counters 
in the arrangement of Fig. 4 and the pair below 
it. Corrections have been made for accidental 
coincidences. However, the coincidence rate due 


to showers probably is considerable, and, since 
the solid angle subtended by the counters is large, 
the observed counting rate is not truly repre- 
sentative of the vertical intensity. The ratio of 


the shower intensity to this measured corpuscular 
intensity increases markedly with altitude as 
shown in Table II, even though the factors 
mentioned above tend to make this increase less 
evident. 


TABLE II. Ratio of showers to corpuscular intensity at 
diffe rent altitudes, 


Barometric pressure (cm) 44 51 64 76 
Relative shower intensity 8.5 50 2.25 1.00 
Relative corpuscular intensity 3.6 2.5 1.46 1.00 
Ratio of showers to corpuscular 

intensity a4 28 83 1.0 


Although the theory as so far developed is 
unsatisfactory in treating multiple pair formation 
by photons, it is of interest to compare the 
theoretical predictions for single pair formation 
with the experimentally determined absorption 
of photons in shower production. Oppenheimer'* 
has presented an argument which would lead to a 
lower limit for the nuclear cross section for pair 
formation by photon encounter. 


18 J. R. Oppenheimer, Phys. Rev. 47, 44 (1935) 
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1 


~ (28/9)aZ*p? In (1/a*Z'e), (2) 
: all Z) 
o~0 for e>a?*Z-3, 

where 


a=e/hc, p=e?/mc*, e=Energy of photon/mce?. 


As a basis of comparison we shall consider the 
absorbing power of a material as a function of its 
atomic number. Since the effect of the logarithmic 
factor is not large, a first comparison can be made 
by neglecting it. Then the absorption per nucleus 
of absorbing material should be proportional to 
the square of the number, and for 
comparison, the atmospheric 
stations have been expressed in equivalent cm of 
lead on the Z? scale. Fig. 6 is a logarithmic plot of 
the counting rate against the total equivalent 
thickness of lead. Each curve rises approximately 
linearly, falls approximately linearly, and then 
levels off. Presumably the linear portion of each 
curve beyond the peak represents the absorption 
of the B-radiation in lead. It is parallel to the 
corresponding portion of each of the other curves 


atomic 


layers between 


and approximately parallel to the line drawn 
through the peaks of the curves and representing 
the absorption of the B-radiation in air.'® This 
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Fic. 6. Logarithmic plot of curves of Fig. 5. The atmos- 
pheric layers are expressed in equivalent cm of lead on 
the basis of an absorption per nucleus proportional to Z?. 


19 W. F. G. Swann, Phys. Rev. 48, 641 (1935), has pro- 
posed an alternative interpretation, based on the assump- 
tion that the increase of shower production with elevation 
may be attributed to an increase (with energy) of the 
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Fic. 7. Shower production curve at sea level extended to 


20 cm of lead. 


indicates that on the scale of the square of the 
atomic number the coefficients of absorption are 
roughly the same for lead and air. However, we 
do not wish to attach too much significance to 
the facts that the portions of the curves appear to 
be linear on a logarithmic plot and that the 
coefficients depend on the square of the atomic 
number. The shape of a logarithmic plot is very 
sensitive to the addition or subtraction of a 
constant value and the departure of the curves 
from linearity at thicknesses of lead 
requires some explanation if this analysis is to 
apply. We considered the possibility that this 
departure might be due to rays which fall 
obliquely on the sides of the lead block and 
produce showers without traversing the entire 
thickness of the block. To test this we built up 
lead wedges on the sides to absorb rays coming in 


great 


obliquely. Since the lead wedges produced no 
appreciable decrease in the counting rate, we 
concluded that this effect small to 
account for the departure of the curves from 
linearity. However, consideration of the experi- 
ments of Ackemann*’ and Hummel*! leads to an 
explanation of the difficulty. By an experiment 
similar to theirs we have verified the existence of 


was too 


a second maximum at 18 cm of lead as shown in 
Fig. 7. The experiments of Kulenkampff” have 
given an explanation of the second maximum. He 
finds that it arises from a very penetrating 
B-radiation which produces a penetrating C- 
probability that a primary ray (either directly or indirectly 

will produce a shower rather than to an increase in abun 

dance of shower-producing rays. 

20M. Ackemann, Naturwiss. 22, 169 (1934). 


217. N. Hummel, Naturwiss. 22, 170 (1934). 
* H. Kulenkampff, Physik. Zeits. 22, 785 (1935). 
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been subtracted from the curves of Fig. 6. 


radiation. We suggest that the showers of the 


second maximum arise from this C-radiation, 
and that the variation with thickness of lead is 
represented by the dashed curve of Fig. 7.7% We 
have assumed, as is strongly indicated by the 
curves of Fig. 5, that similar maxima would be 
found at higher elevations, bearing approxi- 
mately the same ratio to the first maxima. Thus 
from each set of data we subtract the counting 
rate due to showers of this type, and obtain the 
curves of Fig. 8. Beyond the maxima the plots 
are now nearly straight and all have approxi- 
mately the same slope, corresponding to an 
absorption coefficient of about 0.33 per cm of 
lead. This is to be compared with a coefficient of 
0.53 per cm of equivalent lead for the absorption 
of the atmosphere. We estimate the probable 
error to be 0.04 cm for each, which is con- 
siderably less than their difference. The ratio of 
the coefficients (u air/u lead =1.6) may be com- 
pared with the ratio of the logarithmic term for 
air and for lead given by Oppenheimer (Eq. (2)). 
The computed value for an energy of 300 MEV 
is 1.4. More recently an expression for the 
absorption coefficient for high energy photons 

*% The curves Nj; and Noss (of 


difference between 


section 2, Fig. 3) can now be explained as due to the con- 
tribution of showers of the type which cause the second 
maximum to the rate Nosy, but not to the rate Nj, since 
this requires an ionizing ray above the lead. 


RAY 


ic. 8. Curves similar to the dashed curve of Fig. 7 have 





717 


SHOWERS 


has been suggested by Nordheim.** 


N2*r,2 28 1377137mc?\ 3 
p= In ; : (3) 
Z 6 


137 9 € 


w 


where 79>=e?/mc? and N=number of atoms per 
cm’. This expression gives a ratio of 5.0 for the 
logarithmic terms for air and lead. Although the 
experimental errors are large, the indication is 
that the dependence of the nuclear cross section 
on Z is somewhat stronger than that given in 
than 
Nordheim’s. The experimental coefficient 0.33 


Oppenheimer’s expression and weaker 
per cm of lead is to be compared with 0.77 cm™ 
predicted by Oppenheimer and 0.20 cm~' pre- 
dicted by Nordheim. The computations were 
made for an estimated average energy of 300 
MEV.” The uncertainty in the estimated energy 
is large and might possibly account for the 
discrepancy either way. 

To substantiate our conclusion that the ab- 
sorption of high energy photons depends at 
least to a first approximation on the square of the 
atomic number, we have made measurements of 
the absorption of the B-radiation in iron at sea 


level. For this purpose we have placed at A 
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Fic. 9. Curves similar to those in Fig. 8 of the absorption 
of B-rays in lead and in iron (on top of 3.2 cm of lead). 


“1. W. Nordheim, Phys. Rev. 49, 189 (1936). 

*% This estimate is made from a consideration of the 
average number of rays per shower, reference 10, the 
average energy of the shower rays (C. D. Anderson and 
S. H. Neddermeyer, Int. Conf. Physics, London, October 
1934), and the importance of secondary photons, reference 
9, in the phenomenon. 
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bring the C-radiation into equilibrium with the 
B-radiation. Iron placed on top of the lead acts 
essentially as an absorber of the B-radiation. 
The counting rates for lead and iron corrected 
for the second maximum and plotted against the 
total equivalent thickness of lead should lie 
along the same line, and it is found that within 
experimental error they do (Fig. 9). However, 
the correction for the second maximum here is 
rather uncertain. We take its position at 30 cm 
of iron as given by Kulenkampff” and assume 
that the magnitude of the first and second 
maxima for iron are in the same ratio as those for 
lead. For the case of iron and lead the deviation 
from the Z? relationship (i.e., the logarithmic 
term) is too small to be detected by this ex- 


(Fig. 4) 3.2 cm of lead, an amount sufficient to 


periment. 


4. ABSORPTION OF THE RADIATION PRODUCED IN 
Cosmic-RAay SHOWERS 


By placing plates of lead and aluminum in 
position B of Fig. 4 we have measured the 
absorption of the shower radiation from lead 
(probably a mixture of C-, D- and E-rays) in 
these two materials. Small corrections have been 
made for the initial counting rate (with no lead 
at A). It is found that within the experimental 
limits of the investigation the penetration of the 
rays is independent of the thickness of the lead 
at A and of the elevation. Since this is so, we 
have plotted curves for lead and aluminum (Fig. 
10) which represent averages of all the data on 
absorption. Two curves are given for aluminum. 
For one the thickness is taken equivalent to lead 
on the Z scale and the other on the Z? scale. 
Since one of these curves lies above and the other 
below the lead curve, it is seen that the de- 
pendence upon Z is between the first and second 


power. 
The logarithmic plot shows that at least 
empirically the shower radiation can be con- 
sidered as made up of two components with 
coefficients 5.0 and 0.4 per cm of lead. But, if the 
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EQUIVALENT CM OF LEAD 


Fic. 10. Absorption of shower rays in lead and aluminum 
The thickness of aluminum is expressed in equivalent cm 
of lead on the basis of an absorption per nucleus pro 
portional to Z and 2?. 


shower production curve is to be represented by 
Eq. (1), the absorption coefficient v of the shower 
radiation is determined by the value of the 
coefficient 4 of the B-radiation and the position 
of the peak, and is thus found to be about 0.7 per 
cm of lead. This does not agree with the coeffi- 
cient of either of the components deduced from 
direct measurements. It will be noted that since 
the absorption takes place in both plates at B, 
the measured coefficient is too large by a factor of 
two. However, the absorption probably enters in 
the same way in block A and in Eq. (1). The 
coefficient 0.4 may possibly apply to a C-radia- 
tion capable of producing large numbers of soft 
D- and E-rays which are more effective in 
actuating the counters. Further cloud chamber 
studies are required to clear up this point. 

In conclusion the author wishes to express his 
gratitude to Dr. J. C. Street, whose advice and 
inspiration have made this investigation possible. 
Also it is a pleasure to acknowledge the coopera- 
tion of Professor J. C. Stearns and of the Uni- 
versity of Denver in affording laboratory facilities. 
Transportation of the apparatus was generously 
provided by the city of Denver through the 
courtesy of Mr. G. E. Cranmer and Mr. R. 
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On the Geomagnetic Analysis of Cosmic Radiation 


G, LEMAITRE AND M. S. VALLartA,* University of Louvain and Massachusetts Institute of Technology 
(Received March 30, 1936) 


The results of an extensive study of trajectories asymp- 
totic to a certain family of periodic orbits in the earth’s 
magnetic field, carried out by means of Bush’s differential 
analyzer, are presented in this paper. The theory of the 
region of full light, or main cone, is fully discussed. Atten- 
tion is then restricted to the section of the main cone in 
the plane of the geomagnetic meridian and it is shown 


1. INTRODUCTION 


ECENT experiments by the 
multiple coincidences carried out by John- 
son! in Mexico (geomagnetic latitude 29°N) and 
by Clay? in Java (geomagnetic latitude 18°S) 
have shown that in the northern geomagnetic 


method of 


hemisphere the intensity of cosmic radiation in 
the geomagnetic meridian is, for equal zenith 
angles, greater from the south than from the 
north; conversely in the southern hemisphere it 
is greater from the north than from the south.® 
That this is a consequence of the action of the 
earth’s magnetic field on the motion of charged 
particles was pointed out by the present authors 
as early as 1932.‘ Shortly afterwards Bouckaert® 
was able to calculate this north-south asymmetry 
for geomagnetic latitudes up to 20° and moderate 
zenith angles. Considerable difficulties stand in 
the way of extending these calculations to higher 
latitudes and larger zenith angles. The use of 
Bush’s differential analyzer,'' which was made 
available to us to carry out the investigation 
reported here, has made it possible to include in 
the present analysis latitudes as high as 40° and, 
in most cases, zenith angles as far as the horizon. 
Thus with the completion of these researches 
the problem of determining the allowed main 





* At present Visiting Professor of the C.R.B. Educational 
Foundation at the University of Louvain. 

'T. H. Johnson, Phys. Rev. 47, 91 (1935) 

* J. Clay, Physica 8, 867 (1935). 

3A northern excess in Eritrea (geomag. lat. 11.5°N) at 
a zenith angle of 45°, slightly larger than the rather large 
experimental error, has been reported by B. Rossi, Ricerca 
Scient. 5, 583 (1934). 

*G. Lemaitre and M. S. 
(1933). 

5 L. Bouckaert, Ann. de la Soc. Sci. de Bruxelles, A54, 
174 (1934). 


Vallarta, Phys. Rev. 43, 87 


that the north-south asymmetry furnishes the most direct 
approach to the analysis of the energy spectrum in a wide 
region, independently of the particles’ sign. Further it is 
shown that the general shape and the minimum of the 


north-south asymmetry discovered by Johnson in the 


course of his Mexican experiments are fully accounted for 
by the action of the earth’s field. 


cone, or region of full light, for the latitudes 
mentioned is essentially solved. 

The chief results of our present research are, 
first, that the north-south asymmetry, already 
discussed in a previous paper,® depends for each 
latitude on a narrow band of energy. Therefore, 
as further shown in the sequel, the experimental 
study of this asymmetry provides a direct 
workable method for the analysis of the spectral 
distribution of corpuscular cosmic radiation, 
independent of the latter’s sign. Second, we 
that the the 
asymmetry found by Johnson at a geomagnetic 
latitude 29°N and at about 45° zenith angle, as 
well as the general features of his experimental 
results, are fully accounted for by the action of 
the earth’s field. In fact the theory predicts a 


show minimum of north-south 


slight reversal of the sign of the asymmetry for 
zenith angles between 45° and 55° which depends 
on the existence of a very narrow energy band. 
This reversal begins at latitudes around 20° for 
70° and eventually 
broadens out so as to include the whole region 
between the zenith and 40° at latitudes around 
35°. Contrary to our earlier® suggestion it is not 


for zenith angles about 


necessary to invoke atmospheric absorption in 
order to account for this minimum. 

The results we wish to present here must be 
considered as preliminary insofar as a critical 
examination of their precision is concerned. 
There is in fact ground for suspecting that our 
method of finding asymptotic trajectories leads 
to slightly too large a value for the aperture of 


®G. Lemaitre, M. S. Vallarta and L. Bouckaert, Phys. 
Rev. 47, 434 (1935). Full references bearing on our applica- 
tion of Liouville’s theorem to the present problem are 
given in the footnote on p. 435. 
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the main cone. The whole question as to the 
precision of our present results cannot be finally 
settled calculations of 


asymptotic trajectories in which we 


before more rigorous 


are now 
engaged are carried through to completion; we 
hope to return to this question shortly. 


2. OUTLINE OF THE THEORY OF THE 


ALLOWED CONE 


Let us now review briefly the statement of our 
problem and our method of attack. We start 
from the assumption of an isotropic distribution 
of charged particles at very large distances from 
the earth. Because of the action of the earth’s 
magnetic field some of these particles are shot 
back to infinity, others are allowed to reach the 
earth. We have then that at any 
point a given distance from the earth's magnetic 


shown?!: 7 


center all particles of a given energy reaching 
that point must come within a cone, generally of 
very involved shape, which forms the boundary 
the region where all directions are 
the full light region, to borrow from 
and the region where 


between 
allowed 
the terminology of optics 
only some or no direction are allowed. In any 
allowed direction, however, we have shown that* 
the intensity is precisely the same as would 
exist in the absence of the field and hence the 
same as it is at infinity. 

The problem is thus to determine the genera- 
tors of the cone. These generators, as shown by 
one of us’ following ideas already outlined in our 
first paper,‘ may be of either of two kinds. The 
trajectories of the first kind are asymptotic to a 
periodic orbit, in fact asymptotic to periodic 
orbits forming a family first studied by Stérmer.* 
This family has been calculated in detail by one 
of us,’ as well as trajectories asymptotic to the 
members of this family in the immediate vicinity 
of the equatorial plane.’ It was there shown, 
moreover, that some of these trajectories are 
doubly asymptotic, i.e., asymptotic to the same 
periodic orbit in the past as well as in the future. 
It follows further that asymptotic trajectories 
in the vicinity of these doubly asymptotic orbits 
separate again into trajectories going to infinity 

7G. Lemaitre, Ann. de la Soc. Sci. de Bruxelles A54, 
165 (1934). 

8 C. Stormer, Zeits. f. Astrophys. 1, 237 (1930). 


*G. Lemaitre, Ann. de la Soc. Sci. de Bruxelles A54, 
194 (1934). 


AND 


M.S. VALLARTA 


and trajectories returning to the earth. The latter 
may return infinitely near the first asymptotic 
trajectory; further among them is an infinite 
number of other doubly asymptotic trajectories. 
One of us has carried out the analysis of this 
maze of trajectories in the immediate vicinity of 
the equator.’ The consequence is that the full 
light region is surrounded by an infinity of 
infinitely narrow bands of extremely complicated 
structure. This region we have called, borrowing 
again from the language of optics, the region of 
penumbra. 

We have thus been led to distinguish among 
the trajectories of the first kind those which 
directly reach the earth, i.e., without returning 
to the vicinity of the periodic orbit. These 
asymptotic trajectories form the generators of 
what we. have called the main cone or region of 
full light. As mentioned above the latter is 
surrounded by the penumbra. The distinction 
between full light and penumbra is of course 
only a matter of convenience and in no wise is 
to be looked upon as fundamental. 

In addition to the trajectories of the first kind 
which are asymptotic to a periodic orbit there 
are also trajectories of the second kind which can 
form part of the boundary of the main cone.’ 
These trajectories are tangent to the earth before 
reaching the point of observation. They limit a 
part of the main cone because of what we might 
call, again to borrow from the vocabulary of 
optics, the shadow of the earth. 

We have finally to recall briefly some of the 
properties of the differential equations of motion 
of a charged particle in the earth’s dipole field to 
which the above discussion applies. These equa- 
tions have been first set up, studied and made 
generally well known by Carl Stérmer to whom 
the results we are about to recall are also due. 

First by suitably defining the units of time and 
length one may study at once trajectories of 
different particles and correspondingly different 
energies. For a given particle and for each energy 
the standard earth's radius is represented by a 
certain length which measures the energy of the 
particle. We suggest that this unit of length be 
called in the future a Stérmer unit, or briefly a 
Stérmer, in honor of its discoverer. For high 
energy particles this length is very nearly pro- 
portional to the square root of the energy. 
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Second, the motion of a particle may be re- 
solved into two components: one in the meridian 
plane and the other a rotation of the meridian 
plane around the magnetic axis. The calculation 
of this last motion has no bearing on our problem 
and consequently is left out in the sequel.'° If r 
and \ are polar coordinates fixing the position 
of the particle in the meridian plane then the 
space trajectory intersects this plane at an angle 
# given by Stérmer’s equation 


sin 6= —2y,/(r cos \)+(cos A)/Pr*. (1) 


The constant y; is a constant of the motion; it 
is the component of momentum of the particle 
conjugate to the ignorable coordinate of the 
meridian. 

For purposes of calculation it is convenient as 
noted by Stérmer to introduce a new variable x 
by the transformation 


et =2y19. (2) 


If now x and X are looked upon as rectangular 
coordinates a conformal map of the meridian 
plane is obtained where the radius vector r 
becomes a line parallel to the ’ axis and the 
circles r=const. become lines parallel to the 
\ axis. As already explained above the earth is 
represented on this map by a line parallel to the 
\ axis for each value of the energy r in Stérmer 
units. The angle » between a trajectory and the 
zenith direction is therefore read off directly on 
the conformal map since it is the slope defined by 


tan n=dd/dx. (3) 


Finally in terms of an independent variable oc, 
function of the time, the equations of motion 
become 
d*x/do*=30P/dx, d*r/do?=30P/dd, (4) 
where 
P= (dx/da)*?+(dd/da)* 
= e**/16y,;'—(e~-* cos \—1/cos A)? (5) 


and from (1) and (2) it is readily seen that the 
locus of the points for which @=0 is given by 


1°The motion of the meridian plane is important in 
connection with the variation of intensity due to the 
rotation of the galaxy which has been discussed by A. H. 
Compton and I. A. Getting, Phys. Rev. 47, 817 (1935). 
This point has been partially studied by L. A. van Wijk 
and H. Zanstra, Physica 3, 75 (1936). See also C. Stérmer, 
Astrophys. Norv. 1, 157 (1935). 


x= 2 log cos A (6) 


and represents on the conformal map those 
points of the cone lying on the meridian plane. 


3. THE DETERMINATION OF ASYMPTOTIC 
TRAJECTORIES 


The capacity of the differential analyzer" now 
available at the Massachusetts Institute of 
Technology is just sufficient to integrate the 
system of differential equations (4). An ingenious 
set of machine connections was devised for this 
problem by S. H. Caldwell, but for lack of space 
the details of operation must be omitted here. 
Besides a preliminary attempt in the winter of 
1932 the work was carried out at two different 
periods. During the first period (December, 
1933) we found about a score asymptotic tra- 
jectories for 1/16y;'=0.07, y;=0.972. The initial 
slopes were calculated from the first terms of an 
analytical expansion of the asymptotic trajectory 
valid in the vicinity of the periodic orbit. This 
method of attack provided valuable information 
on the general appearance and properties of 
asymptotic trajectories and gave one of us the 
incentive for the introduction of harmonic ex- 
pansions such as those used later by himself? 
and by Bouckaert.® It suffers from the serious 
disadvantage that the trajectory is started from 
a region of instability particularly sensitive to 
small errors either of the machine or in the 
initial data. During the second period (spring of 
1935) we preferred, for this reason, to start the 
trajectory outwards from a suitably chosen 
point on the equator and then to adjust the 
initial slope so that the trajectory would neither 
intersect nor fall short of the periodic orbit, 
which was already accurately known.’ It was 
found possible to adjust the initial slope after a 
few runs with a precision usually around one- 
thousandth of a radian. When the initial slopes 
had been found for several points on the equator 
the initial data at intermediate points were 
determined by graphical interpolation and a 
score trajectories were started inwards. Later it 
became obvious that certain interesting tra- 
jectories would be missed in this way and so, 
instead of starting from the equator, appropriate 


“For a full description of the differential analyzer and 
its operation see V. Bush, J. Frank. Inst. 212, 447 (1931). 
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points were chosen on the line 6=0 and tra- 
jectories were run from these points outwards. 
The initial slope was then adjusted so that the 
trajectory would be tangent to the envelope of 
the asymptotic family, which had been previ- 
ously determined from those asymptotic tra- 
jectories already known. As an additional verifi- 
cation a few trajectories were run as far as the 
periodic orbit. After the initial data at inter- 
mediate points had been found by graphical 
interpolation as in the case of the equatorial 
runs, a score Or more asymptotic trajectories were 
continued inwards. We have thus studied the 
0.81, 0.83, 0.85, 0.87, 
and 0.99. Runs from 
made for the following 


following values of y;: 
0.89, 0.91, 0.93, 0.95 
equatorial points were 
intermediate even values of y,: 0.88, 0.90, 0.92, 
0.94, 0.96, 0.98. In all three hundred 
asymptotic trajectories were traced. 

As mentioned in the introduction, a critical 
study of the precision of the results thus obtained 


some 


must in the final analysis await the outcome of 
calculations now in progress. We should mention 
at the outset that fairly systematic checks based 
on the use of the energy Eq. (5) do not reveal 
serious systematic errors. Only in one instance 
did we find a small systematic error which was 
later traced to incorrect plotting of one of the 
known functions the machine. 
When this mistake was corrected and a part of 
the previous work was repeated, the observed 


introduced in 


changes in the trajectories were so small that it 
was not deemed necessary to repeat all of the 
previous runs, some twenty-five trajectories in 
all. One can, however, mostly judge of the 
accuracy of the present material from the point 
of view of its self-consistency. Asymptotic tra- 
jectories possess, as a matter of fact, a compli- 
cated system of envelopes to be more fully 
discussed below and even small accidental errors 
are immediately visible. It should be borne in 
mind that the final precision of a curve traced 
by the differential analyzer depends not only on 
the mechanical trustworthiness of the machine 
itself, but also on the accuracy of the laborious 
and often complicated numerical calculations 
required to introduce the initial conditions into 
each of the members of a coupled system of six 
integrating units and five input tables. It even 
depends on the care and faithfulness with which 
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a crew of five assistants keep in step their 
respective cross hairs on the known functions 
plotted on metallic plates on the input tables. 
Our guiding principle for the reduction of the 
data obtained with the help of the differential 
analyzer has thus been to accept only such 
conclusions as are based on the study not of a 
single, but of a good many consistent trajectories. 
4. GENERAL DESCRIPTION OF AsymPTOTIC TRA- 


JECTORIES AND THEIR ENVELOPES 


After these considerations we pass on to a 
description of some general features of asymp- 
trajectories. To begin with we _ should 
mention the trajectories calculated by Carl 
Stérmer™ by the method of numerical integration 
These 

G@=() 


totic 


which reach the pole r=0 or x=— <=. 
trajectories are osculating to the line 
already mentioned above. Although most of 
them belong to a family different from the one 
we are concerned with here, there are a few 
which are members of both, in particular one 
for y, between 0.9313 and 0.9314 and a few 
which are of importance for the study of the 
penumbra. Stérmer has shown that for sufh- 
ciently small values of x all other trajectories 
oscillate about the nonoscillating orbit going to 
the pole. This is precisely the general character 
of the trajectories we have obtained. They 
oscillate about Stérmer’s nonoscillatory orbit; 
for decreasing values of x the loops come closer 
together and the trajectory oscillates outwards 
in the same manner as inwards. An example of 
a family of asymptotic trajectories is shown in 
Fig. 1, where for the sake of completeness the 
locus of the vertices of the family of periodic 
orbits and the line 6=0 are also drawn. 

A family of asymptotic trajectories possesses 
an extensive family of envelopes which we are 
now about to describe (Figs. 1 and 2). Beginning 
at the periodic orbit there is an envelope denoted 
by E» which first proceeds in a direction fairly 
parallel to the x axis and then rises rapidly 
towards a cusp Cy where it meets a corresponding 
envelope proceeding upwards, which we denote 
by Fy. The system of envelopes continues in a 
similar manner; corresponding branches are 


2C, Stérmer, Terr. Mag. and Atmos. Elec. 37, 
(1932). Astrophys. Norv. 1, 145 (1935). 
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Fic. 1. The family of envelopes and periodic orbits, and a 
family of asymptotic trajectories (+,;=0.972). 


denoted by Ci, C2 --+, Fi, Fe, «++. These enve- 
lopes are tangent to the line P=0 (not drawn 
in the figure) at points which are denoted by S 
for the E envelopes and 7 for the F envelopes. 
At any S or 7 there is an asymptotic trajectory 
turning back on itself, i.e., a self-reversing doubly 
asymptotic trajectory. On either side of S or T 
an asymptotic trajectory touches the envelope 
and then proceeds in directions away from S or 
T. In particular the point Sp is the vertex of the 
periodic orbit. The corresponding vertex for 
negative value of \, and in general the corre- 
sponding points of tangency below the x axis 
will be denoted in the following by primed 
symbols. 

The trajectories are divided by their points of 
contact with their envelopes into sections (Figs. 
1 and 2). According to a known theorem in the 
theory of kinetic foci'* any trajectory infinitely 
near an asymptotic orbit intersects each one of 
the sections once only in any finite domain of 
variation of ¢. All asymptotic trajectories tend 





81t is well known that the equation of the distance 
between trajectories infinitely close to some trajectory and 
this trajectory, as a function of the arc measured along 
the trajectory in question, is linear, homogeneous and of 
the second order, cf. E. T. Whittaker, Analytical Dynamics 
(Cambridge, 1904), p. 382. By a known result due to 
Sturm the zeros of such an equation separate one another, 
cf. M. Bocher, Lecons sur les méthodes de Sturm (Paris, 
1917), p. 46. See also H. Poincaré, Méthodes nouvelles de la 
mécanique céleste, Vol. 3, p. 261. 
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; 
Fic. 2. Northern and southern sections of asymptotic 
trajectories, 


towards periodic orbits while performing an 
infinite number of oscillations between the 
envelope Ey and its symmetrical counterpart Ey’. 
The trajectory in its infinitely close vicinity will 
therefore intersect each of the sections into 
which it is divided by its points of contact with 
its envelope. One may therefore choose the 
neighboring trajectory sufficiently close to the 
asymptotic orbit that the former may intersect 
the latter’s sections as many times as desired 
but finally the former will pull away from the 
latter either by cutting the periodic orbit or by 
returning inwards. We thus have to distinguish 
between two kinds of sections (Figs. 1 and 2): 
northern sections (dashed) are those for which 
an infinitely close trajectory intersecting the 
asymptotic trajectory at an angle smaller than 
its own slope at that point finally cuts the 
periodic orbit. They form the northern boundary 
of the full light region or main cone. Southern 
sections (dotted) are such that an infinitely close 
trajectory intersecting the asymptotic trajectory 
at an angle greater than its own slope finally cuts 
through the periodic orbit. They form the 
southern boundary of the main cone. We have 
concluded from a study of the trajectories drawn 
by the differential analyzer that the sections 
coming from Zp» are northern sections, those 
coming from E,’ are southern sections. If we 
denote their slopes by n and ¢, respectively, we 
have the result that all trajectories included 
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Fic. 3. Schematic representation of main cone. 


between » and ¢ come from infinity. The region 
between » and ¢ lies within the main cone of 
which, in agreement with the results outlined in 
a preceding section, asymptotic trajectories are 
the generators. Each time that an asymptotic 
trajectory has a point of contact with its envelope 
the sections bounded by this point of contact 
are of different kinds. Sections of asymptotic 
trajectories diverging from the first cusp Cy are 
therefore northern sections. After having touched 
the envelopes ~; or F, they become southern 
sections and therefore the sections diverging 
from the cusp C; are southern sections. In the 
region between the two cusps Cy) and C; we may 
therefore have three trajectories 
through the same point. We denote by 7 the 
slope of the trajectory coming from the cusp Co, 
by m the slope of the trajectory tangent to FE, 
and by ¢, that of the curve tangent to F,. At 
the next cusp C, all three slopes n, m and & 
become the same. Moreover, even before reach- 
ing C, our trajectories may have touched other 
envelopes £2 or F2 returning with slopes m2 and £2. 
What conclusions are to be drawn from these 
facts regarding the structure of the main cone? 

Let us assume for instance that we have 
obtained, for a certain point x, \ and a given 
value of y; the directions 7, m, f1, m2, f. The 
latter belong to the given value of the latitude 


asymptotic 


\ and are related to the energy r and the angle @ 
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Fic. 4. Typical trajectories of second and third kinds 
by (2) and (1), respectively. For positive values 
of X, i.e., in the northern hemisphere and for 
positive particles @ will be taken positive east- 
wards. We have adopted for our representation 
of the cone (Fig. 3) the orthogonal projection on 
the horizontal plane of its trace on the unit 
sphere centered at the point in question. In 
accordance with the conventions stated above we 
take the coordinates of our representative point, 
i.e., sin @ and cos @sin » as positive eastwards 
and northwards, respectively. Then the cone 
will have the general appearance shown Sche- 
matically in Fig. 3, the full light region being 
to the west. Let us now see the significance of 
the envelopes and cusps. The envelope £2 is 
such that on it m and m coincide. Generally 
envelopes are double points on the trace of the 
cone and there is a tangent parallel to the 
meridian plane. Cusps are triple points where 
the trace of the cone has a point of inflection 
with a tangent parallel to the meridian plane. 
These general principles suffice to interpret 
asymptotic trajectories in order to determine 
the main cone or region of full light. 

It remains to show how the shadow of the 
earth limits a part of the main cone. In order 
that the preceding considerations remain valid 
it is necessary of course that all the asymptotic 
trajectories we have considered remain outside 
the earth. Hence if any of them has anywhere a 
tangent parallel to the A axis it retains its 
previously discussed significance only for suffi- 
ciently small values of the energy, less than that 
which corresponds to the position of the tangent. 
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At a point such as A (Fig. 4) the direction 
gives a generator of the main cone, but at B the 
latter is already bounded by a trajectory of the 
second kind. At a point such as O the main cone 
is bounded by a trajectory which is at once of 
the first and of the second kinds. Such points 
indicate where the region of shadow, or region 
of the second kind, begins. A trajectory which is 
at once of the first and of the second kinds we 
shall briefly refer to as a trajectory of the third 
kind. 

The complete main cones, as well as a fuller 
discussion of the azimuthal effect than has been 
hitherto possible, are reserved for a forthcoming 
companion paper. 


5. CALCULATION OF THE NORTH-SOUTH 
ASYMMETRY 


Directions in the meridian plane correspond 
on the conformal map to the line +=2 log cos \ 
Eq. (6) for which 6=0 by Stérmer’s formula (1). 
It is quite feasible to measure for each asymptotic 
trajectory cutting through this line the slope 7 
and the corresponding value of the latitude \. 
These values of » and A can then be plotted on 
a diagram, for each value of y;, and a smooth 
curve drawn through the plotted points (Fig. 5). 
As is to be expected in accordance with our 
previous discussion a few of these points fall 
clearly outside the graph indicated by the other 
points, but in general the plotted curves are de- 
termined without ambiguity. For regions where 
sufficient data are missing it is easy to interpolate 
between neighboring well-determined curves. 
Particular attention should be paid in connection 
with this diagram to the significance of the 
envelopes and cusps. The intersection of an 
envelope with the line @=0 is shown as a point 
where the tangent to the plotted curve is parallel 
to the » axis. These are the ‘“‘loci of maximum 
latitude”’ shown on the diagram. A point where 
a cusp lies on the line 6=0 is a double point of 
the curves plotted. Its characteristic feature is 
that two branches of the locus y;=const. inter- 
sect, one of them having a tangent parallel to 
the axis. This is clearly so because on one side 
of the cusp, as already discussed elsewhere, 
there are three asymptotic directions while on 
the other there is but one. The value of y, for 
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Fic. 5, The zenith angle in the meridian plane as a function 
of latitude. 


which the first cusp C) lies on the line 6=0 is 
somewhere near y,=0.94. The following cusps 
C,, Co, ete., lie practically on it for y,=0.93. 
Further, Stérmer’s results quoted above show 
that these values form a progression the limiting 
value of which is 0.9313. 

The trajectories of the third kind are important 
because they mark the point where the earth’s 
shadow begins to appear. It is usually possible, 
after examination of the family of asymptotic 
trajectories drawn by the differential analyzer, 
to decide with good approximation where a 
trajectory of the third kind appears and to 
interpolate among neighboring curves to find the 
angle at which it intersects the line @=0. The 
locus of the values of zenith angle and latitude 
for which a trajectory of the third kind appears, 
for different values of y;, is shown in Fig. 5. It 
should be carefully borne in mind, in interpreting 
this locus, that a region may be in shadow for a 
given value of y; and not for others. Thus for 
instance a point in shadow for y,=0.93 is still 
in the region of light for y,=0.81 at \= 20°. 

The relation connecting y, and the energy r 
for points on the meridian plane at a geomagnetic 
latitude X is, from (1) 


r=(cos* \)/27. (7) 


By using this formula and the results of the 
previous analysis the least energies required for 
a particle to arrive at any given zenith angle in 


the meridian plane, at any point on the earth up 
to latitudes of 40°, may be readily found. It is 
simply necessary to read off the A diagram 
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Fic. 6. North-south asymmetry. 


(Fig. 5) the value of y: corresponding to the 
given » and X, interpolated if necessary. Typical 
curves relating the least energy of arrival with 
the zenith angle for latitudes 0°, 10°, 20°, 30° 
and 40° are reproduced in Fig. 6. It should be 
remembered that the energy (in millistérmer 
units) is measured with respect to the standard 
earth’s radius (6370 km) and hence for observa- 
tions at a given point on the earth account must 
be taken of the eccentricity of the earth's 
magnetic center.'* Several important features 
are to be emphasized in connection with these 
curves. First it will be noted that, irrespective of 
the particles’ sign, the energy region explored 
by the north-south asymmetry varies consider- 
ably with the latitude. The geomagnetic field 
thus provides for the separation of energies 
required for a spectral analysis of the corpuscular 
cosmic radiation. Experiments on this asym- 
metry thus furnish the most direct approach to 
the study of the energy spectrum. For this 


4 M.S. Vallarta, Phys. Rev. 47, 647 (1935). 
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reason it is earnestly to be hoped that they will 
be continued particularly in the region between 
35° and the geomagnetic equator, both at high 
and low altitudes above sea level. In this way 
valuable information will be obtained not only 
in view of the purpose mentioned above but also 
for the study of atmospheric absorption and 
scattering of particles of different energies. 

Secondly it will be seen that the theory gives a 
full account of the minimum found by Johnson! 
in his Mexican experiments at a zenith angle of 
around 45°. It goes even further and, in a roughly 
quantitative way commensurate with present 
experimental precision, accounts for his whole 
result. 

Lastly the effect of the earth’s shadow should 
be carefully noted. It results in a large north- 


south asymmetry at high zenith angles for 
intermediate latitudes. For low latitudes the 
asymmetry near the horizon should therefore 


be reduced. 
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A Note on the Possible Effect of Screening in the Theory of Beta-Disintegration 


M. E. Rose, Institute for Advanced Study, Princeton, N. J. 
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Che present paper considers the possibility that screening by the atomic electrons has an 


appreciable effect on the energy distribution of 8-particles emitted by heavy nuclei. A general 


formula is derived from which it is possible to conclude that the effect of screening is negligible 


rhis result is also derived by another method 


Finally, an explicit calculation based on the 


model of a charged sphere is shown to lead to the same conclusion. 


i. 

N a recent paper Konopinski and Uhlenbeck! 

have remarked that, on the basis of the 
Fermi theory of §-disintegration with their 
modified interaction, one finds for heavy ele- 
ments such as Ra E an energy distribution which 
predicts too many slow electrons as compared 
with the experimental results.* An attempt by 
these authors to remove this difficulty by other 
modifications of the interaction between the 
heavy particle and the electron-neutrino field was 
not successful in diminishing the number of slow 
electrons without at the same time losing the 
degree of asymmetry in the distribution curve 
required by experiment.* Although there may be 
some question as to the existence of a real 
discrepancy between experiment and the theory, 
as modified by Konopinski and Uhlenbeck, it is 
of some interest to inquire whether there may be 
some influence hitherto not considered which 
could appreciably affect the energy distribution. 

In this connection it is plausible to consider the 
possible effect of screening by the atomic elec- 
trons. It is evident, from a qualitative stand- 
point, that the screening will alter the distri- 
bution curve in the direction of the “free 
particle’ distribution’ which, incidentally, does 
fit the experimental curve for Ra E as given by 
Sargent. In the following we shall determine the 
effect of screening quantitatively. 
‘sphere’ (whose 


By enclosing the system in a 
size is afterwards allowed to increase without 
limit) we may consider the wave function y of 


1 E. J. Konopinski and G. E. Uhlenbeck, Phys. Rev. 48, 
7 (1935). 

2B. W. Sargent, Proc. Camb. Phil. Soc. 28, 538 (1932). 

’ This feature of the experimental results is apparently 
more reliable than the results for the low energy electrons. 
See also reference 11. 

‘That is, the distribution obtained by neglecting the 
Coulomb, field of the nucleus. 


the 8-particle to be of a discrete type. The 
screening can be represented by a Thomas- 
Fermi potential® 

V=—(aZ/r)®(yr), y=a(128Z/9r)). (1) 


It is to be noted however, that the main argu- 
ment given below depends on the use of the 
Thomas-Fermi potential only secondarily. 

We now select a point 7) such that for r>r a 
solution of the wave equation of the W. kK. B. 
type will be valid,® and for r<rp the argument of 
® in (1) is small enough to allow the approxi- 
mation? 

@=1-—y7r. 


Then for r<7r the wave functions will be 
Coulombian with the energy W replaced by 
W—V> where the constant screening potential 
Vo=aZy. The phase of the wave function at 1% 
will be denoted by x and the increase in phase 
from 7% to the classical limit of motion (or to the 
boundary of the sphere) by @. Further we label 
each state with an index nm (radial quantum 
number). Then 


Xnet— XntOng1—-8,=m0 or d(x+60)/dn=r7. (2) 


We then apply adiabatically an infinitesimal 
perturbation at small r producing an energy 
change 

SBW=)AS F(r) |W) Ad tance. (3) 


where F(r) is an arbitrary function of r confined 
to nuclear dimensions and \ is very small. The 
integration in (3) is to be taken over the volume 
of the nucleus. 


5 We use rational relativistic units, length being measured 
in h/mc, energy in mc? and momentum in mc. In these 
units e? =a (fine structure constant). 

* Only differences in phase as given by Eq. (2) below 
rather than the exact values of the phases have to be given 
correctly by the W. K. B. formulas. 

7 For Z~80, y =0.036 while the singularity in the W. K. 
B. solution is at rl. 
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In general x and @ could depend on \ and on W. 
Since the W. K. B. method is supposed to apply 
for r>7 and since F=0 in this region one may 
consider @ as depending only on W and not on d 
while x naturally depends on both \ and W. The 
perturbation formula (3) is supposed to apply to 
a constant v, i.e., to such changes of x, @ that 

6x+60=0. 
Since @ is a function only of W and not of A this 
means that 
(dx/0W+00/0W)bW= —ddx/OAX. (4) 
Substituting for 6W and 0(x+8@)/dW by means 
of (2) and (3) we have 
(dn/dW) S Fir) WW) *dtnuc. = —(1 


The right side of this equation is the same for a 


w)dx/OX. (5) 


given F as long as the wave equation is the same. 


Thus it is the same for an electron of energy W 


in the field of a screened nucleus and an electron 
of energy W— 1% in the field of an unscreened 
nucleus. Such states of screened and unscreened 
nuclei may be called corresponding states and 
the essence of the discussion lies in the possibility 
of establishing such a correspondence. The left 
side of the equation is the product of the number 
of states per unit energy range dn/dW and a 
weighted mean of the density |y/*. Since F is 
arbitrary (dn/dW)|y\* at any point in the 
nucleus is the same for corresponding screened 
and unscreened states. This latter quantity we 
shall refer to as the electron density per unit 
energy range. 

It should be noted that du/dW and |y/* 
change by large factors if the standard normaliza- 
tion to the same value of r*| y |? at infinity is used. 
Also the function F can be generalized to an 
operator which may involve the spin indices of 
the function ¥ without altering the above proof. 
It thus follows that 


cr 


(dn/dW) Xy,*Guw, (2°) 
is the same for corresponding screened and 
unscreened states for any point inside the 


nucleus and for an arbitrary Hermitian G,,. 
The essential quantities which enter in the 

theories of 8-decay are all of the form (5’)* and 

can be 


7 


the above consideration 


the result of 
expressed without formulas by saying that the 
electron distribution is always such as though the 
~ 8 Cf. E. Fermi, Zeits. f. Physik 88, 161 (1934). 
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nucleus were not conscious of the screening and 
as though it emitted electrons into its immediate 
vicinity always in the same way; the only effect 
of the screening is then to accelerate the electrons 
by the Coulombian repulsion of the screening 
electrons and to change the kinetic energy W— V 
into W by this acceleration. From this point of 
view the result has superficially an appearance of 
being self-evident. Actually it 
because it is possible to find examples in which it 


requires proof 


is not correct. One such example is given in 
section 3 of this note. In fact conclusions reached 
in (5) and (5’) may be inexact because inter- 
ference effects in the space between the nucleus 
and the screening charge can change the density 
y *. These interference effects are neglected in 
the W. K. B. method used in obtaining (5). In 
the applications to 8-decay it is reasonable to use 
such an approximation because the screening 
potential changes little within a wave-length of 
the electron. 

Proceeding to the disintegration formula we 
write the Konopinski-Uhlenbeck result for the 
probability of emission of an electron with 
energy between W and W+dW as 


Pxru(W)dW=W(W,—W)*u(W)dW, (6) 
W, being the maximum emission energy. Then it 
follows from the above considerations that the 
effect of screening is to change this to 
P(W)dW 
= (W— V.)(W.— W)4u(W— Vo)dW. 
Since for Z~80, Vy=0.02, the change in P is of 
Hence one may 


(6’) 


less. 


the order 3 
conclude that the effect of screening is neg- 


percent or 


ligible.® 
2. 


the conclusion 
expressed in (5) by another method in which 
explicit use of the W. K. B. solution is made. We 
denote by ¢)/r and ¢2/r the radial functions of 
the Dirac y¥ which can be made to fulfill the 


normalization condition 


R 
f (| ¢1 eo ¢2 *\dr=1, (7) 
0 


where R is the radius of a large sphere in which 
the system under consideration is enclosed. By 


It is possible to arrive at 


*I am indebted to Professor G. Breit who pointed out 
to me the derivation of this result. 
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the substitution 
ve=(1+W-—V)' oo 

one finds that v2 satisfies the differential equation"® 
d*v2/dr*>+f"(r)ve=0. (8) 


To (8) the W. K. B. method of solution may be 
applied. With 


Ri gi(2)|)*+ | oof 2) 2)|=1 
we find 
go( ~)=[(T+1)/2RT]! sin (kr+phase), (9) 


where 7=W-—V and the wave number & is 
given by 


T?=k?+1. (10) 


In a region where the W. K. B. solution just 
begins to be applicable we have 
T’+1 kk} 
go= sin (k’r+phase), (11) 
2RT k’ 
where k’ is the wave number and T” the kinetic 
energy in this region. 
Since R can be made much larger than any 
other length entering the problem, the phase 
shift in the wave function is essentially given by 


RAk =r. (12) 
From (10), (11) and (12) we find that 
go|?/AT=(1+T7")/2rk’, (13) 


which is independent of k and 7. Hence keeping 
k’ the same it follows that |y|?/AT is inde- 
pendent of the screening. Since the screening has 
the effect of changing k’ in accordance with the 
energy shift W—-W-— V4, it is evident that the 
electron density per unit energy range is given 
by the corresponding expression in the case of no 
screening with the energy W replaced by W— V». 


3. 


An alternative way in which we could repre- 
sent the screening is to use a charged sphere 
potential. 


V=—aZ(1/r—1/r.) for r<rz, 
V=0 for r>r.; 7,.=1/7. 


_™Cf. G. Breit, Phys. Rev. 38, 463 (1931). The function 
f*(r) is defined on p. 470 of this paper. 
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With this model the calculation of the disinte- 
gration probability may be carried out with 
exact wave functions and thus it is possible to 
take into account the interference effects arising 
from the reflection of the electron waves by the 
charged sphere. The result one obtains is not the 
disintegration probability (6) but a more 
complicated expression which contains (6’) as a 
first approximation. 

For r<r, the wave function is again the 
Coulombian one with the energy displaced by an 
amount V» while for r>r, linear combinations of 
Bessel functions are used. The normalization 
constant entering in |y/* will depend on the 
wave functions at the point 7,. For large Z the 
“atomic radius” 7, is large compared to unity 
(Z ~80, r,= 28) sothat we may use the asymptotic 
expansions of the wave functions. The dis- 
integration probability thus obtained is 


P.(W)dW=P(W)[1+0(r.2) WW, (14) 


where P(W) is given in (6’). The correction term 
in (14) denoted by O(r,-*) arises from terms 
beyond the first in the asymptotic expansions of 
the wave functions. These higher terms are of an 
oscillatory nature and an average over a range of 
large values of the radius r, has been carried out. 
Actually the term O(r,-*) contributes a correction 
much smaller than the relative difference be- 
tween P(W) and Pxu(W) so that to a first 
approximation the screening is represented in 
the case of the present model in the same manner 
as before." 

I wish to express my sincere thanks to Pro- 
fessor G. Breit for his many courtesies to me. 


1 Since the completion of these calculations there has 
appeared a notice of a new investigation of the continuous 
8-ray spectrum of Ra E by A. I. Alichanow, A. I. Alichanian 
and B. S. Dzelepow, Nature 137, 314 (1936). In these 
experiments particular attention has been given to the low 
energy portion of the distribution and it was found that a 
large number of slow electrons are emitted. These authors 
have made measurements starting at an energy (not 
including the rest energy) of 30 kv and have obtained a 
distribution curve which appears to be in agreement with 
the Konopinski-Uhlenbeck theory. The appearance of 
their curve suggests that the probability does not tend to 
zero at the origin, as is the case with the Sargent curves, 
but rather to a finite value comparable with the maximum 
probability. With these new measurements as a criterion 
in the case of heavy nuclei and with the results for light 
nuclei, see especially F. N. D. Kurie, J. R. Richardson and 
H. C. Paxton, Phys. Rev. 49, 368 (1936), it seems that the 
Konopinski-Uhlenbeck theory is capable of accounting for 
all the features of the continuous 8-ray spectrum. 
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The Spectra of Xenon in the Extreme Ultraviolet 
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Increased dispersion and resolving power have made 
possible a revision and extension of the work of previous 
investigators on the spectra of xenon in the extreme 
ultraviolet. Lists of classified lines are given for Xe I, Xe II 
and Xe III in the range from 42000 to \600. No new lines 
were found for Xe | but the accuracy of the measurements 
Revised identifications are 


is considerably improved. 


HE extreme ultraviolet spectra of xenon 

have been investigated by using the two- 
meter normal incidence vacuum spectrograph! 
of the Carnegie Institution of Washington, which 
is located in the Spectroscopy Laboratory of the 
Institute of The 
the plates and the 


Massachusetts 
method of reduction of 
estimated precision of the been 
discussed with the results on the 
standards used have been published.’ The lines 
newly or differently identified in the present 
investigation are denoted by an asterisk. The 


Technology. 


have 
and 


data 
neon” 


presence of a Roman numeral in the intensity 
column denotes that the line is a blend with a 
line of that spectrum. Previous measurements on 
xenon in this spectral range were made by 
Abbink and the 
present investigation was purchased from the 
Air Reduction Company. Exposures were made 
different total 
pressures of a mixture with neon. 
Lines due to different stages of ionization were 
distinguished by the variation of their intensity 


Dorgelo.t The gas used in 


of electrodeless discharges at 
of xenon 


with pressure. 
Xel 

Classified lines are given in Table I. As this 
type of discharge does not favor the excitation of 
first spectra only a few lines were observed. The 
notation used for the upper levels from which 
these lines arise is that given by Humphreys and 
Meggers® in their term table. The ground term 


> 


1K. T. Compton and J. C. Boyce, Rev. 
(1934). 

2 J. C. Boyce, Phys. Rev. 46, 378 (1934). 
3]. C. Boyce and C. A. Rieke, Phys. Rev. 47, 653 (1935). 


Sci. Inst. 5, 218 


4J. H. Abbink and H. B. Dorgelo, Zeits. f. Physik 47, 
221 (1928). 

5C. J. Humphreys and W. F. Meggers, Bur. Standards 
J. Research 10, 139 (1933). 


suggested for Xe II in the light of new knowledge of the 
limits of that spectrum, the low states of Xe III. 128 lines 
of Xe III have been identified in collaboration with C. J 
Humphreys, who will publish elsewhere with his analysis 
of the longer wave-length portion of the spectrum the 


consolidated term table for this spectrum. 


is now located at 97834.4 cm~ as related to the 
other term values given in that table. This 
makes no appreciable change from the earlier 
results and the ionization potential, 12.078 volts, 
already published® remains the same. 


Xe Il 

A partial term table for this spectrum has 
already been published by Humphreys, de Bruin 
and Meggers.’ Present knowledge of the limits 
the Xe ITT, 


indicates that the previous assignment for the 


of this spectrum, low states of 


TABLE I. Xe J classified lines. 


X INT v CLASSIFICATION 
1469.621 5 68044.7 Se Is°4 
1295.560 &d 77 186.7 Ss 1 
1250.199 ) 79987.2 So-3 
1192.04 2 83890.2 Ss—3 

TABLE II. Xe JJ classified lines. 

X INT y CLASSIFICATION’ 
1244.756 5 80337.0 p> 2p. p* S 
1183.053 7111 $4527.1 st psp - P) 6s 41 
1169.63 2 85469.9 p> 2P% P) 5d 4D 
1158.474 5 86320.5 e p> 2p) ‘P) 6s 41 
1100.432 10 90873.4 p>2P s ps s 
1083.860 5 92262.8 2 ps 2p — (P) 6s *Py 
1074.476 15 93068.6 Pps2P°n — (8P) 6s 4Poy 
1051.920 10 95064.3 2 pe 2p? — (§P) 6s *Py 
1048.272 & 95395.1 s? p52P°, — (3P) 5d 4D 
1041.306 9 96033.3 s? ps2P°, — (§P) 5d 4D 

*1037.680 6 96368.8 s? ps2p% — (§P) Sd 2D 
1032.438 4 96858. 1 2 pi2P° — (P) 6s *P 
*976.678 6 102388 2 ps2P%, — (1D) 6s *Dy 
972.769 7 102799 ps 2p? — (P) 6s *Pij 
*935.405 2 106906 p® 2P°n — @P) 5d2*Dy 
*925.866 5 108007 2p? — ‘P) Sd 2D2 
*885.54 3 112925 p> 2P°1 — OD) 6s *Dy 
*880.802 5 113533 p> 2P°1, — (OD) 6s *Dy 
803.066 3 124523 P°s S) 6s 2S; 

740.406 3 135062 Pen S) 6s 2S; 


®W. F. Meggers, T. L. de Bruin and C. J. Humphreys, 
Bur. Standards J. Research 3, 731 (1929). 

7C. J. Humphreys, T. L. de Bruin and W. F. Meggers, 
Bur. Standards J. Research 6, 287 (1931). 
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ance Ill. Xe J//] classified lines 


CLASSIFICATION 


Pa 
7 


1978.702 3 50538.2 ‘S) 5d 4D°. — (2P) 67 3D, ? (20 
1921.630 2 5 2039.2 ta 1P°. — (4S) 6p 3P 
*1896.904 5 52717.4 ‘S) 5d3D°,; — (2P) 6p °Ds ? (18 
*1874.907 5 53336.0 ‘S) 5d3D°. — 26 

*1854.365 6 53926.8 ‘S) §5d5D°, — Ss 
*1834.254 4 54518.1 4S) 5d 5D°. — 82 
*1826.477 7 §4750.2 4S) 5a5D°; — (2D) 6p 2 (6 
*1825.858 3 54768.8 ‘S) 5d 51°, — Re 
*1819.845 6 54949.7 ‘S) 5d5D°; — (*D) 6p ¥F ? (10 
*1817.381 2 55024.2 4S) 5d °D°2 — 122 
*1804.109 2 §5429.0 p P°, — (D) 6p 2D, 
*1750.749 2 57118.4 p P°; — (2D) 6p *D2 
*1579.498 5 63311.2 p P?, — 12 
*1562.563 4 63997.4 p sP°s — (2D) 6p *F 
*1554.438 1 64331.9 p P°s — (2D) 6p *D 
*1549.975 2 64517.2 p® 3P?, — (2D) 6p *P i 14) 
*1511.121 3 66176.0 p> o — 2s 
*1509.454 2 66249.1 p 3P°. — (2D) 6p 'P1 4 
*1468.180 4 O8111.5 s p> Pp — 82 
*1457.356 3 68617.4 p +P — 122 
*1432.204 Shi 69882.4 p P — (2—)) 6p *P1 141) 
*1393.496 1 71762.0 sp P?, — (*P) 6p*Di ? (20, 
*1356.364 3 73726.5 sp 3P°s — (2D) 6p'De ? (162 
*1253.645 2 79767.4 p P — 26 
*1232.074 25 81164.0 ss? p* De —s f° P 
*1225.089 3 81626.7 2 pt Ss —sp a 
*1183.053 SIl 84527.1 p* So — (48S) 5d 3D% 
*1156.480 9 86469.3 > ps De. —sp P 
*1130.344 30 88468.6 5s? pt P, —sp P 
*1066.391 12 93774.3 stp aP; -sp sp; 
*1058.128 2 94506.6 2 ps ID» 4S) 5d 5D 
*1055.328 5 94757.3 sp 1D. — (4S) 5d 5p? 
*1048.754 3 95351.3 2 ps De — (8S) 5d5D*, 
*1047.801 10 95438.0 s2 pt Po —sp5 Ps 
*1017.680 35 98262.7 s?p* %P,y —s ps 1P?s 
*1010.376 361 C II 98973.1 s* p4 So — (2D) 5d iP 
*1003.370 35 99664. 1 5? ps Py -sp I Po 

*98 1.088 7 101928 pt 61D, — sp PP) 
*974.124 & 102656 5? ps Py — (S) 5d5D° 

*97 1.818 8 102900 2 ps Py — (CS) 5d51% 

965.540 10 103569 s? ps ap —sp PrP?) 
960.325 2 104131 p* ID, — (4S) 5d 3p? 

958.585 4 104321 s2 ps 3Po — (4S) 5d Dp 1 
*953.975 3 104825 2 pt 1De — (48S) 5d 3p° 
*917.257 $ 109021 2 ps De — (2D) 5d 4F°s 2? (1 
*915.488 3 109231 s? pt Pi -sp 1p; 
*901.746 7 110896 stp! 3%Po —s pd P?; 

*ROR R73 & 111250 2 pt De — (2D) Sd 3F°. 205 
*896.003 20 111607 2 ps PP: — (4S) 5d 5p° 
*895.406 4 111681 pt 3Py — (4S) 6s 58 
*893.989 20 111858 ps P: — (4S) 5d 5D 

*29 1.833 9 112129 stpt 3Py — (48) 5439 
*889.276 15 112451 stp %P2 — (4S) Sd 5D, 
*878.790 8 113793 s2 ps Po — (S) 5d3pD°, 
*870.346 6 114897 stpt 3P, — (4S) 5d 3% 
*863.386 8 115823 s? ps 1Py — (4S) 6s 38% 
*861.071 5 116134 spt De — (2D) 6s 3D, 
*852.950 25 117240 stpt So — (2P) 6s 1P%, 
*851.147 8 117489 stpt %Po — (4S) 6s 38% 
*850.572 5 117568 st pt 'De — (2D) 6s 3D%s 
*840.162 7 119025 stp* %P, —s p> Ip? 
*838.449 3 119268 s? ps IDs — (2D) 5d '1P°, 
*826.134 4 121046 s? ps ‘(De — (2D) 5d 3p, 
*824.881 30 121230 s? p* 3Ps — (48S) 5d 4p; 
*823.210 25 121476 ps P2 — (4S) 6s ®S°2 


(‘D)6s 2D term was in error. It would now seem 
better to identify these terms as (*P)5d 2D. 
Another term has been found which is more 
satisfactory for the first level. Table II lists the 
20 classified lines of Xe II which can be identified 
with some certainty. Publication of additional 
lines belonging to this particular spectrum of 
xenon, as well as of any revised term table, is 
reserved pending further investigation. There is 
no reason to alter the approximate value, 21.1 
volts, given for the ionization potential by 


Humphreys, de Bruin and Meggers. 


d In CLASSIFICATION 
*822.647 4 121559 2 ps De — (*D) 6s *D 
*820. 166 4 121926 stpt 3P_ — (4S) 5d 2D 
*810.119 7 123430 st pt Py — (2D) 6s *D% 
*808.860 3 123631 ps De — (2D) 5d 4S% 
*801.980 15 124691 2 ps Pe — (4S) 5d4D°%s 
*800.819 3 124872 stpt 3P, — (2D) 6s *D°s 
*800.228 2 124964 s? pt De — (2D) 5d 4D 
*799.338 ~ 125104 stp! Po — (2D) 6s *D% 
*796.070 12 125617 s? pt 3P>s — (48) 6s 38% 
*793.977 & 125948 st ps IDe — (2D) 6s 1D°%s 
*793.292 8 126057 2 pt De — (2D) 5d 5D 
*792.896 15 126120 pt %P, — (2D) Sd 3F° 1 
*790.064 5 126572 stpt 3Py — (2D) Sd Ps 
*780.030 7 128200 stp! ID, — (2D) 5d'1D 
*779.781 5 128241 s2 ps Po — (2D) 5d'P* 
*77¢ ? 5 R40 s? pt *Py - (2D) Sd #D* 
shear 2 128349 <2 pt 3P, — (2D) $d 3I < 
*769.143 10 130015 stpt 8Po — (2D) 5d 4D% 
*765.120 7hl 130698 s? ps De — 13% 
*763.736 7 130935 tpt 3Py — (2D) 5d 4S 
*761.790 5 131270 s? p* De — 15 
*761.540 10 131313 2 ps De — 17° 
*756.031 10 132270 stp* 3Py — (2D) 5d5p 
*754.144 5 132601 s? ps ‘Po — (2D) $d 48%: 
*750.447 8 133255 stp! 3Py — (2D) 6s 1D 
*750.155 8 133306 s? ps D 19 
*744.141 6 134383 st pt De — (*P) 6s *P* 21° 
*742.566 15 i 34668 s? ps Ps — (2D) 6s *D®: 
*737.979 7 135505 s? pt Py 2—)) 5d 'D*% + 
*733.314 10 136367 ps Ps — (2D) 5d'P* 
*731.030 15 136793 s? ps De — (*P) 6s *P 
*727.058 9 137541 s? p* De — (@P) 6s 'P* 
*724.623 3 138003 stp! %Pi — 13% 
*723.873 5 138146 s? ps 3P,. — (2D) 5d 4D 
*721.630 4d 138575 s? p* Py — 15 
*721.199 10 138658 s? ps P: — (@D) 6s §D 
*7 15.986 4 139668 s? ps Po — 13° 
*711.190 + 140609 stp! 3P, — 19 
*7 10.677 5 140711 sips 3P, — (2P) 6s 3P 
*710.576 5 140731 s? p* Ps — (2D) S§d 4S 
*705.783 5 141687 st pt » —(2P) 6s 3P% 1 
*705.096 12 141825 tpt IDs — 23% 
*703.906 9 142064 stp* 3P, — (2D) 5d 4D 
*702.799 x 142288 stp! ID, — 25% 
*699.070 5 143047 s? ps iP, —(*D) 6s 'D 
*698.541 20 143156 s2 ps P - (°D) 5d 4p 
*697.526 ~ 143349 s? pt Po — (CP) 6s *P% (21 
*693.972 10 144098 stp! 3P, — (@P) 6s §P 
*691.036 7 144710 st?p* =1De — 27% 
*690.397 7 144844 stp! %Py — (2P) 6s Ps 
*688.231 $ 145300 stp! %P, — (2D) 5d'D 
*682.564 7 146506 s? p* Po — (*P) 6s Ps 
*676.606 9bl 147796 s? pt IP. — 13°, 
*673.996 9 148369 s? p4 iP. — 15 
*673.813 9 148409 s? pt P: 17 
*670.55 2 149132 apt  3Py — 23 
*668.473 149595 st pt Py — 25° 
*664.877 6 150404 stpt 3P, — 19% 
*661.124 4 151258 8 =s?pt =3Po — 25° 
*660.124 8 151487 s? pt ip — (*P) 6s *P*; 21°) 
*657.828 8 152015 spt 8, = 27's 
*646.667 5 154639 stpt 8P, — (°P) 6s 'P%1 
*629.217 7 158928 s? p* Ps —23°s 
"627.393 4 159390 s2 ps 3P2 — 25% 


Xe III 

The principal group of lines was located by 
analogy with the s* p'*P—s p®*P° groups of 
Ne III, A III and Kr III. The ground P term is 
partially inverted as in tellurium I and iodine 
II.*» ° Other multiplets having differences char- 
acteristic of the ground state were then found to 
give connections with many odd terms, including 
two terms already found by de Bruin.’® An 


8. Bartelt, Zeits. f. Physik 88, 522 (1934). 

® P. Lacroute, Ann. de physique 3, 5 (1935). 

1 T, L. de Bruin, Pieter Zeeman Jubilee Volume (Mar- 
tinus Nijhoff, The Hague 1935), p. 413. 
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analysis by C. J. Humphreys of the visible and 
near ultraviolet data, carried out simultaneously 
with the present investigation, made possible a 
rather complete identification of these odd terms. 
The exact description of many of the 128 lines 
classified in Table III was made possible by the 
very generous cooperation of Dr. Humphreys. 
All of the low terms of Xe III have been located 
and many intersystem combinations have been 
found. The consolidated term table for Xe III 


will be published by Humphreys elsewhere"! 

1 C, J. Humphreys, Bur. Standards J. Research (in 
press). 
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together with his analysis of the longer wave- 
length portion of this spectrum. The ionization 
potential is 32.0 volts. 

The experimental part of this work was made 
possible by a grant to Dr. K. T. Compton from 
the Carnegie Institution of Washington. Some 
of the technical assistance necessary in plate 
measurement and reduction has been provided 
by a grant to the writer from the Permanent 


Science Fund of the American Academy of 
Arts and Sciences. The careful work of Dr. 
Carol A. Rieke in this capacity is gratefully 


acknowledged. 
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The Christiansen Filter Effect in the Infrared 


R. BOWLING BARNEs AND LYMAN G. BONNER,* Palmer Physical Laboratory, Princeton University 


(Received March 23, 1936) 


A brief discussion is given of the history, construction, 
and applications of Christiansen light filters for use in the 
These 


regions of the spectrum. 


visible and ultraviolet 
filters, which consist of a powder of some crystalline or 
glassy substance suspended in a transparent medium, 
depend for their action upon the equality of refractive 
indices of powder and medium for some one wave-length of 
light. The theory of the effect in the infrared is briefly 
discussed, and methods are given for preparing and 
mounting the powders. The principal factors affecting the 
shape and size of the transmission peak, namely, the angle 


at which the two dispersion curves cross, the size of 


ITHIN the past few years, several papers 

have appeared in the literature in which 
various uses of Christiansen light filters have 
been discussed. One commercial company has 
advertised in the Review of Scientific Instruments 
that it can furnish various crystals in a powdered 
form suitable for use in the construction of 
these filters. In Germany the firm of Schott 
and Genossen-Jena offers for sale such filters 
under the name ‘‘monochromatic universal 
filters.”” Although the Christiansen filter effect 
has been known since 1884 it has lain dormant 
until the appearance of these papers except for 
brief descriptions which have appeared in each 


* National Research Fellow. 


particle used, the thickness of the powder layer, and the 
solid angle included by the light striking and leaving the 
filter Experimental results include trans 
mission 
maxima for the following crystalline powders dispersed in 
air: quartz, MgO, calcite, marble, CaSO,, eleven alkali 
thallium 


are discussed. 


curves and wave-lengths of the transmission 


halides, and three monovalent halides. Curves 
are also given for quartz, MgO, and NaCl dispersed in 
organic liquids. The peaks measured lie 
between 3u and 90x. In conclusion are discussed possible 


applications of the Christiansen effect in the infrared. 


transmission 


edition of Wood’s Physical Optics. Even today 
the effect as such is not known to the average 
physicist. It is the purpose of the present paper 
to review the subject very briefly and then to 
describe the construction and several uses of such 
filters in the near and far-infrared. It will be 
shown that these filters may be used (1) to 
produce monochromatic beams of infrared radi- 
ation, (2) to eliminate the undesired short wave- 
length radiation in a grating spectrometer, (3) to 
determine the refractive indices of various fluids 
in the near-infrared. 

If curve a of Fig. 1 represents the dispersion 
curve of a crystal, and curve 0} that of a given 
transparent medium, then a suspension of small 
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FILTER EFFECT 


particles of this crystal in the medium should 
possess a maximum of transmission at the wave- 
length corresponding to the point of intersection 
of the two curves a and b. This wave-length we 
shall designate as Acnr. The reason for this 
transmission obviously lies in the fact that only 
at this wave-length is the sample optically 
homogeneous. Light of this wave-length passes 
through the cell unaffected, while light of 
all other wave-lengths suffers refraction at 
every particle-medium interface and is scattered 
strongly. In the case that Acn, lies in the visible 
region of the spectrum, directly transmitted light 
has the color corresponding to this wave-length. 
The direct image is surrounded by a haze having 
the color complimentary to Acur. 

C. Christiansen! first noticed this effect in 1884 
while studying particle-liquid mixtures in an 
attempt to determine the refractive indices of 
various crystals. He explained the effect correctly 
and showed that light of practically any color 
could be obtained by varying the curve 0 of 
Fig. 1 so that it crossed a at the desired wave- 
length. He further pointed out that since the two 
curves a and 6 respond differently to changes in 
temperature, the point of intersection, would 
shift if the temperature of the cell were varied. 

In 1885 Lord Rayleigh? published a paper 
entitled “On an Improved Apparatus for 
Christiansen’s Experiment.’”’ From then until 
1927 the filter effect was unused. At that time 














| Acar. 





A- 


Fic. 1. Curve a represents the dispersion curve of a 
crystal, and 6 that of a transparent medium. The two 
intersect at an angle @ at the point A=Achr. 


1 C, Christiansen, Ann. d. Physik 23, 298 (1884); 24, 439 
(1885). 
* Lord Rayleigh, Phil. Mag. 20, 358 (1885). 
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Weigert’ and his co-workers published the first 
of a series of papers in which attention was 
directed to the fact that with proper temperature 
control and regulation such filters could be used 
,juantitatively for the production of monochro- 
matic light. 

By correct choice of the substances involved, 
Ache may be shifted from the visible to the 
ultraviolet or infrared regions of the spectrum. 
Christiansen, Staude, and others had previously 
suggested this possibility, and in 1932 Kohn and 
von Fragstein* published an account of the 
construction of a Christiansen filter for the 
ultravialet. By using particles of fused quartz 
and a mixture of 56 percent benzol and 44 percent 
ethyl alcohol they obtained suitable filters having 
Achr=3600A at 20°C. With filters of this type 
they were able to isolate any desired region from 
3000A on up into the visible. Individual lines of 
the mercury spectrum could be filtered success- 
fully. These filters were also used by von Frag- 
stein in his work on the determination of the 
reflection power ol metals in the visible and 
ultraviolet regions. 

In a still more recent paper, McAlister® has 
presented a rather thorough account of ‘‘The 
Christiansen Light Filter: Its Advantages and 
Limitations.’’ Practical suggestions are given 
concerning the construction of a filter which may 
be used in an intense beam of light for the 
production of monochromatic light. The filters 
were made of borosilicate crown glass and a 
solution of carbon disulphide in benzene at 20°C. 
A series of filters covering the entire visible 
range, each having a half-width of about 100 to 
200A were prepared. Fig. 2 is taken from 
McAlister’s paper. Since a complete description 
is given by McAlister of the various choices of 
materials, particle size, etc., which are most 
advantageous for filters designed for use in this 
part of the spectrum, no further discussion of 
these points will be presented in this paper. 

Other recent papers on the Christiansen effect 
include one by Abbot® on the applications of 
~ 3F. Weigert and H. Staude, Zeits. f. physik. Chemie 
130, 607 (1927); F. Weigert, H. Staude and E. Elvegard, 
ibid. B2, 149 (1929); F. Weigert and J. Shidei, ibid. B9, 329 
(1930). 

*H. Kohn and K. von Fragstein, Physik. Zeits. 33, 929 
(1932). 

5E. D. McAlister, Smithsonian Misc. Collections, Vol. 


93, No. 7 (1935). 
®C. G. Abbot, J. Opt. Soc. Am. 25, 42 (1935). 
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Fic. 2. Transmission curves of a series of five filters con- 
sisting of borosilicate crown glass particles iramersed in 
mixtures of carbon disulphide and benzene (McAlister). 


Christiansen filters to stellar energy spectra and 
one by Denmark and Cady’ on the optimum 
grain size for these filters for use in the visible 
spectrum. 

In order to demonstrate the principles involved 
in the construction of Christiansen filters the 
authors have found the following procedure both 
simple and expedient. A glass absorption cell 
having plane parallel windows is first filled with 
powdered fluorite (CaF2) whose refractive index 
Np» = 1.434. Care should be taken to insure that 
the particles are clean and dry. Benzene, whose 
refractive index mp=1.501, is now poured into 
the cell until the latter is roughly half full. 
Again powdered fluorite is added to refill the 
cell. To this ethyl-alcohol mp=1.362 is added 
and the cell thoroughly shaken. This addition is 
continued until the light transmitted by the cell 
takes on a reddish hue. Each additional drop of 
alcohol the refractive index of the 
liquid causing Acnr to shift shorter 
wave-lengths. Sufficient alcohol should be added 
to cause the lamp 
filament to be clear yellow. The scattered light 
will be the complementary color to deur. 

It is easy to see in the case of our fluorite- 
benzene-alcohol cell that the value of Ac, is a 
function of the temperature. If the cell be cooled 
in a stream of tap water, Acnr will shift to a deep 
red. Upon returning to room temperature Achr 
will again be yellow. Placing the cell over a 


decreases 
towards 


transmitted image of a 


7H. S. Denmark and W. M. Cady, J. Opt. Soc. Am. 25, 
330 (1935). 
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steam radiator will cause Acne to be blue or 
violet. A temperature either too high or 
low will cause the color to disappear Acar being 


either in the ultraviolet or the infrared. 


too 


THE CHRISTIANSEN FILTER EFFECT IN 
THE INFRARED 


Many investigators* have studied the trans- 
mission of various powdered crystals mounted in 
such a way that the surrounding medium was 
air. These powders have invariably shown a 
rather high transmission peak on the short 
wave-length side of the characteristic frequency. 
No attention was paid to this abnormally high 
transmission although the curiosity of the in- 
vestigators should have been stirred by the fact 
that a layer of powder shows a trend just oppo- 
site to that shown by a thin single crystal of the 
same material. Fig. 3 shows the general type of 
curve found in the Pfund showed 
that, in case the average particle size was 
larger than X», the wave-length of the character- 
istic absorption, this high transmission peak 
occurred at the wave-length at which the re- 
flection of the crystal was a minimum, and 
explained the effect on this basis. The reflection 
power of a crystal is given, however, by the 


literature. 


formula 

R=[((n—1)?+ m2? ]/[(n+1)?+ n°? ]. 
At a wave-length well removed from Xo, since « 
is nearly equal to zero, R=(n—1)*/(n+1)*. We 
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Fic. 3. Typical transmission curves of plate and powder. 
a is the transmission of a crystalline plate and b that of a 
powder of the same substance. 


$A. H. Pfund, J. Opt. Soc. Am. 23, 375 (1933); A. H. 
Pfund and S. Silverman, Phys. Rev. 39, 64 (1932); S. 
Tolksdorf, Ber. Diss. 1928; O. Reinkober, Zeits. f. Physik 
3, 1 (1920). 
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Fic. 4. Typical dispersion curve of a crystal for wave- 
lengths in the neighborhood of A». At Acnr the index of the 
crystal is equal to that of air. 


see then that this transmission peak also corre- 
sponds with the wave-length at which the re- 
fractive index of the crystal is unity, hence that 
at this point the dispersion curve of the crystal 
crosses that of air. It was pointed out by one of 
the authors® that this abnormal transmission is 
to be explained by a lack of refraction and not 
of reflection, and is in reality a true Christiansen 
filter effect. In Fig. 4 is shown a characteristic 
dispersion curve for a crystal in the neighborhood 
of Ay. In the case of particles even a thin layer 
will be opaque due to refraction over most of 
the range. At some point, however, the indices 
of crystal and medium, which is air in this case, 
cross as is shown in the figure and a rather high 
transmission results. 


PREPARATION OF FILTERS 


The powder to be studied is usually dusted 
onto some material transparent in the desired 
wave-length region. Rocksalt plates may be used 
for the base material in the near infrared, while 
for wave-lengths >30u thin sheets of paraffin 
(~0.5 mm in thickness) have been found quite 
satisfactory. Nitrocellulose films may be used 
over the whole infrared although some difficulties 
arise as a result of the characteristic absorption 
of the nitrocellulose, and trouble is frequently 
experienced in making the powder stick to the 


*R. B. Barnes, R. R. Brattain and R. S. Firestone, Phys. 
Rev. 47, 792 (1935). Since in regions removed from the 
characteristic frequency the reflection power of the crystals 
studied is usually about five percent, reflection losses alone 
could not account for the opacity found for powder layers 
only two or three particles thick. 


IN THE INFRARED 


film. This last difficulty may sometimes be 
avoided by making a thick suspension of the 
powder in water, pouring this onto the film, and 
allowing the water to evaporate. In the near 
infrared a second rocksalt plate may be put over 
the particles if necessary, while in the case of 
paraffin as base material the powder may be 
pressed gently into the surface of the paraffin. 

As will be shown later, the size and degree of 
separation of the powder is important in deter- 
mining the characteristics of the Christiansen 
peak. For best results the particles should be all 
larger than XAcpr, to increase the short wave 
opacity. On the other hand, if they are too large 
the true absorption of the crystal will cut down 
the size of the transmission peak. The particles 
should therefore be only slightly larger than cu, 
in diameter, and the separation should be sharp. 
In the experience of the authors, this separation 
may be most easily performed in one of two ways. 

If the particles are insoluble in water they may 
be separated by allowing them to fall through a 
known height of water. Stoke’s law, assuming 
spherical particles, may be used to calculate the 
velocity of fall for the various particle sizes 
desired. The procedure is simple. The material 
is ground to as fine a powder as possible in an 
agate mortar and then placed in a cylinder of 
water high enough to give a convenient time of 
fall. After a thorough shaking they are allowed 
to settle for a time equal to that calculated for 
particles slightly larger than those desired to fall 
the full height of the vessel. At the end of this 
time the water and smaller particles are quickly 
decanted and the process repeated several times 
on the remaining powder. The particles that 
have been poured off are now taken and the 
time of fall increased to that for particles slightly 
smaller than those desired. The decantation 
process is continued until a sharp line of de- 
markation is seen between powder and clear 
liquid during the time of fall. A check with a 
microscope has shown in several cases that the 
diameter calculated, even with the assumption 
of spherical particles, was not far different from 
that obtained. Samples have been prepared in 
which the diameter ranged from 10 to 12y, 124 
to 15u, etc. For quartz particles in water the 
times to fall one inch for 1p and 100u particles 
are 28,200 and 2.82 seconds, respectively. 
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Fic. 5. An apparatus for separating powders according 
to particle size by an air stream method. 


Particles soluble in water may be separated by 
the same method if a suitable liquid can be found. 
Otherwise the separation may be effected by 
allowing them to settle out of an air blast under 
suitable conditions. In Fig. 5 is diagrammed a 
method which we have found convenient in 
separating some of the alkali halides. A steady 
stream of air is blown into the tube through the 
nozzle. Particles are carried out of the jet by the 
air stream and enter the box, which serves to 
reduce cross currents. Receptacles placed at 
various points catch particles of different sizes, 
the smallest landing the greatest distance from 
the jet. The separation may, of course, be 
regulated also by the strength of the air blast. 
The separation here is not as sharp as that 
obtained by the other method, and the necessary 
conditions must be determined empirically for 
each size and variety of particle, but even so 
separations satisfactory for most purposes may 


be made. 


DETERMINING THE CHARACTERISTICS 
OF THE PEAKS 


FACTORS 


There are, in all probability, three principal 
factors regulating the sharpness and height of 
these Christiansen transmission peaks and these 
will be taken up in order. 

(1) The width of the peak will obviously be 
some function of the angle (@ of Fig. 1) at which 
the dispersion curve of the crystal crosses that 
of the medium. If the angle is small refraction 
at the surface becomes small for a wave-length 
appreciably removed from Ac», and the peak is 
broad, while for a large angle the peak will 
tend to be sharper. 

(2) Particle size, as was mentioned above, and 
the thickness of the powder layer play quite 
important parts in determining the character- 
istics of the transmission peak. In case any large 
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fraction of the particles are of diameter p<)cur, 
the layer will not be completely opaque for 
\<Acnr due to a simple powder effect.'® On the 
other hand since, for the case of particles in air, 
Achr is always in the neighborhood of and on the 
short wave side of \y, the true absorption of the 
crystal is appreciable. If, then, any large fraction 
of the particles are of diameter p> cn, this true 
absorption may reduce the percentage trans- 
mission at the center of the peak considerably. 
For the same reason the thickness of the powder 
layer is important. Microphotographs of filters 
which we have used have shown that a layer 
only 2 or 3 particles deep will be quite opaque 
in the region \ <Acur. 

(3) Since the action of these filters depends, 
not on absorption of the undesired light, but 
simply on scattering, the efficacy of the filter 
will be in part determined by the angle at which 
the light strikes the powder layer, and by the 
solid angle included in the observation method. 
The maximum efficiency will of course be ob- 
tained if only light emerging perpendicular to 
the face of the filter is used. We have found the 
filters to be quite satisfactory, however, even 
when used at the focus of an optical system of 
aperture f : 2. Some improvement may be ob- 
tained in such a case by using two filters spaced 
such a distance apart that the light scattered by 
the first is well removed from the light beam 
before it strikes the second. 


EXPERIMENTAL RESULTS 


Apparatus 

Two types of apparatus were used in per- 
forming the experiments reported in this paper. 
For the region 1—15u a rocksalt prism spec- 
trometer was used. In this region were investi- 
gated principally the Christiansen peaks resulting 
from anomalous dispersion in the neighborhood 
of the inner vibration frequencies of complex 
crystals, and also the peaks for crystals embedded 
in media other than For investigation 
between 30 and 130u a wire grating spectrometer 
of the type described by Barnes'' was used. The 


air. 


spacing constant of the grating was 0.4 mm, but 


1 A. H. Pfund, Phys. Rev. 36, 71 (1930). 
u R. B. Barnes, Rev. Sci. Inst. 5, 237 (1934). 
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since in order to facilitate measurements wide 
slits were used the full resolution of the spec- 
trometer was never obtained. In this long wave 
region were found principally the peaks due to 
lattice vibrations of the alkali halide crystals. 
In all cases percentage transmission curves were 
obtained by taking readings, at each wave- 
length, of the energy with the sample in, and 
out of, the light beam. The samples were placed 
at the focus of a strongly converging beam. 


Results 


(1) NaCl. Of the systems investigated in the 
present work, the one about which the most data 
from other sources have accumulated, and which 
gives complete confirmation of the theory of the 
Christiansen effect that of NaCl 
immersed in CS,. The peak for rocksalt in air 
should occur at 32y, as is shown both by the 


is 


powder 


dispersion and reflection curves, but because of 
experimental difficulties we were unable to 
verify this. The peak for the NaCl-CS, system, 
however, falls in a readily accessible region, and 
in one in which the dispersion curves of both 
and medium are well known. The dis- 
persion of CS: has recently been determined by 
Pfund,” while the values for rocksalt have long 


crystal 
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Fic. 6. NaCl powder immersed in CS». The transmission 
curves are a, the dry powder; 5, the powder wet with CS»; 
c, pure CS». At the top are shown the dispersion curves of 
d, CS., and e, NaCl, showing that they intersect at \=Achr- 


? A.H. Pfund, J. Opt. Soc. Am. 25, 351 (1935). 
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Fic. 7. Effect of particle size on the quartz Christiansen 


peak. The curves are for the following sizes: a, 5-10,; 


b, 10-20u; c, 30-40u; d, 50-60u; e, unseparated powder 
Curve f shows the reflection power of a quartz plate. 


been known from Paschen’s work on the disper- 
sion of a prism of the material. The cell for this 
investigation was made by putting a thin layer 
of finely powdered NaCl between rocksalt plates, 
sealing three of the edges with water glass, and 
flowing CS, into the space between. Fig. 6 shows 
the curve obtained. a is for the dry NaCl powder, 
b for the powder wet with CS», and c is the 
transmission curve of pure CS», taken from 
Pfund. At the top are shown the dispersion 
curves for the two materials, and it may be seen 
that they cross exactly at the wave-length, 5.45y, 
at which the Christiansen peak occurs. 

(2) Quartz. Powdered quartz was available in 
great quantity, and since the material is quite 
insoluble in water and so could be readily 
separated by the rate of fall method, most of 
the experiments on the effect of particle size 
have been made on this crystal. One of the 
fundamental of quartz, which 
probably an inner vibration of the SiO, tetra- 


vibrations is 
hedron occurring in the structure, lies at about 
9u. The Christiansen peak, expected just on the 
short wave side of this falls experimentally, for 
the powder in air, at about 7.35y. Fig. 7 shows 
the transmission curves for approximately equal 
layers of quartz particles of various thicknesses, 


mounted on nitrocellulose films. The small 
irregularities on the sides of the curves, which 
become particularly marked for the 5-—10u 


particles, result from the true absorption of the 
nitrocellulose. The effect of particle size in 


narrowing and decreasing the height of the peak 
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is quite well shown by these curves. A curve is 
also shown for quartz particles which had under- 
gone no separation whatsoever, and it may be 
seen from this that even a random selection of 
particles shows the Christiansen effect quite 
markedly. The dotted reflection curve shows 
that the wave-length of the Christiansen peak 
agrees with that of the reflection minimum. In 
Fig. 8 the effect of immersion in CCl, and in a 
50 percent by volume mixture of CCl, and CS. 
on the quartz peak is shown. The positions of 
these peaks agree well with the known dispersion 
curves for the substances involved. 

(3) MgO. Fig. 9 shows the positions of the 
Christiansen peaks for unseparated MgO powder 
in air and in CCl,. The dispersion curve of CCl, 
was also obtained by Pfund,” and so the estab- 
lishment of these two values, m=1 at 12.2u, 
1.43 at 9.0u, together with previously known 


Percent Transmission 


| 


‘-wwewvwT?e &FfVetTt eS 
A wy 





Fic. 8. Shift of the Christiansen peak for quartz powder 
in liquids. a, quartz in a 50 percent by volume mixture of 
CS, and CCl; 6, quartz in pure CCl,; c, quartz in air. 


TABLE I. Positions of the Christiansen peaks for various 
crystals, 


Wave-LeNGTH (in y») 
RESTRAHLEN 








CRYSTAL CHRISTIANSEN FUNDAMENTAL 
LiF 11.2 26 32.6 
NaCl (32) 52.0 61.1 
NaBr 37 74.7 
Nal 49 85.5 
KCl 37 63.4 70.7 
KBr 52 81.5 88.3 
KI 64 94 102.0 
RbCl 45 73.8 84.8 
RbBr 65 114.0 
RbI 73 129.5 
CsCl 50 102.0 
CsBr 60 134.0 
TIC! 45 91.9 117.0 
TIBr 64 117 





values at shorter wave-lengths, permit a rough 
dispersion curve to be drawn for MgO. This 
curve also is shown in the figure. 

(4) Carbonates. From reflection measurements 
it is known that all crystalline carbonates have 
three fundamental frequencies, characteristic of 
the CO;~ radical, in the near infrared. Experi- 
ments were performed both on powdered calcite 
and on marble and the results, though not 
conclusive, indicated the presence of two peaks 
of 6.2 and 11.14 with possibly a third at about 9u 

(5) Sulfates. Only one sulfate was investigated, 
namely, plaster of paris (CaSO,-1/2H.O). A 
small but definite transmission peak, character- 
istic of the sulfate radical, was located at 7.9u 

(6) Alkali- and thallium-halides. The Chris- 
tiansen peaks in air were located for a total of 
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Fic. 9. a, MgO powder in CCl,, 6, MgO powder in a 
At the top is shown a rough dispersion curve for MgO as 
determined by these two points and by measurements 
the visible spectrum. 
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Fic. 10. Curves showing the Christiansen transmission 
peaks for powdered KCI, KBr, and KI. These curves are 
typical of those obtained for the alkali halides. 
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Fic. 11. Here the wave-lengths of the Christiansen, reststrahlen, and fundamental frequencies 
are plotted, for each halogen separately, as a function of the atomic weight of the metal ion. 


eleven alkali halides, and also for three mono- 
valent thallium halides, which also crystallize in 
a cubic lattice. With the exception of that for 
Lik these peaks all fell in the region of the grating 
spectrometer and so were investigated as unsepa- 
rated powders mounted on 0.5 mm paraffin. 
this thickness has quite a high 
transmission in the infrared beyond 30. Crystals 


Paraffin of 


of this type possess only one characteristic 
frequency, and hence only one Christiansen 
peak. Table I gives the positions of the peaks 
for all the crystals of this type investigated. 
Included in the table for purposes of comparison, 
are the wave-lengths of the fundamental fre- 
quencies and reststrahlen frequencies for the 
same crystals, as given by Barnes.'* Fig. 10 
gives the actual measured transmissions for the 
powders of KCI, KBr, and KI, while in Fig. 11 
the three Christiansen, 


types of frequency, 


13 R. B. Barnes, Zeits. f. Physik 75, 723 (1932). 


reststrahlen, and fundamental are plotted, for 
each halogen separately, as a function of the 
atomic weight of the metal ion. 


THE CHRISTIANSEN EFFECT IN 
THE INFRARED 


USES OF 


Below we shall discuss briefly several inter- 
esting and important applications of the Chris- 
tiansen effect in connection with infrared spec- 
troscopy. 

(1) In the use of grating spectrometers one of 
the most difficult problems that arises is that of 
the elimination of overlapping orders. With most 
gratings, in order to make measurements at a 
wave-length \ all energy of wave-length less than 
r/2 
selective sources and receivers, filters, prism 


must be removed. Many devices such as 


monochromators, etc., have been used. In a good 
many cases the desired effect may be attained 
more simply by the use of a suitable Christiansen 
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filter. In the near infrared powders of quartz, 
LiF, MgO, and doubtless many other crystals 
have sharp and conveniently located peaks. In 
the far infrared the alkali halides form a series, 
together with the monovalent thallium halides, 
covering quite a range of wave-lengths. By 
proper choice of particle size these filters may be 
made quite opaque in the short wave regions 
but will still possess a high transmission in the 
neighborhood of the peak. 

(2) The production of comparatively mono- 
chromatic beams of light by means of these 
filters immediately suggests their possible use 
in the construction of a survey spectrometer for 
the far infrared, similar to the reststrahlen 
method introduced by Rubens. Since only a 
source, two mirrors, and a receiver would be 
required the apparatus would be simple. Whereas 
the reststrahlen method requires one or more 
good sized crystals of each material, which are 
often difficult or impossible to obtain, only a 
small amount of the crystal powder is needed 
to make one of these filters. Using only the 
substances investigated in the present work, 
some fifteen or twenty different wave-lengths are 
available between 6u and 90u, and it is very 
likely that this list will be extended by further 
work in the field. 
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(3) Since the very existence of the Christian- 
sen effect depends upon the equality of the 
refractive indices of particle and medium for 
some wave-length, we have here a method for 
determining dispersion curves in the infrared. 
In case the dispersion of either crystal or medium 
is known, the measurement of the Christiansen 
peak gives one point on the dispersion curve of 
the other. Since it is usually easier to obtain the 
dispersion curves of crystals in the infrared than 
of liquids, a series of such crystal powders could 
be used to establish points on the dispersion 
curves of liquids. In the inverse case, the dis- 
persions of some liquids are known, and in the 
previous section were shown curves of the 
transmissions of MgO, NaCl, and quartz im- 
mersed in CS. and CCl. 

(4) It has recently been pointed out to us by 
Professor E. U. Condon, that since the measure- 
ment of the Christiansen peak in air gives one 
point on the dispersion curve of the crystal, 
namely, that at which n=1, this value may be 
used in the classical dispersion equation, together 
with values from the visible and ultraviolet 
regions of the spectrum to determine the in- 
tensity constant characterizing the infrared ab- 
sorption band of the crystal. This problem will 
be discussed in a later paper. 
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The spectrum of boron fluoride as excited by the 
electrodeless discharge has been photographed at both high 
and low dispersion. Sixteen new bands were observed in 
the region between AA3712 and 2580A. These bands 
degrade to the violet; have five head-forming branches, 
and resemble very closely the third positive bands of CO, 
a result to be anticipated from the fact that CO and BF 
both possess ten extranuclear electrons and should therefore 


INTRODUCTION 
HE BF molecule is isoelectronic with Ne 
and CO, hence considerable interest at- 
taches to its band spectrum. Largely because of 
difficulties in exciting and isolating the band 


have similar electronic configurations. Two systems were 
found with a common lower electronic level: an A system 
(v9. =32,020.5 cm™') of 12 bands and one isotope band, 
and a B system (v'°:® =37,991 cm™) of three bands. The 
constants obtained from the vibrational analysis are: 
we’, 1631.40 cm; w,’’, 1323.64 cm; x,.’w’, 23.10 cm 
xe'w"", 9.40 cm. 


spectrum of BF, it has been inadequately in- 
vestigated. The chemically stable compound of 
boron and fluorine, BFs, is a gas at room tempera- 
ture, and a Geissler discharge in BF; has been 
used by several investigators in searching for BF 
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bands. Johnson and Jenkins! attempted to excite 
the BF spectrum in BF; gas but were troubled 
with the appearance of SiF bands. Under better 
conditions of purity, Tawde and Johnson? ob- 
tained bands in the visible region of the spectrum 
which they ascribed to BF. However, the bands 
did not show the expected resemblance to CO 
bands. 

Experience with the electrodeless discharge in 
this laboratory* showed that the band spectrum 
of a diatomic fragment of a polyatomic molecule 
is brought out strongly in this source, hence the 
project was undertaken of subjecting BF; to this 
type of discharge for the purpose of obtaining 
BF bands. This technique also offered to simplify 
the problem of eliminating impurities. 

While this investigation was in progress, Dull‘ 
reported obtaining three bands of BF in the 
ultraviolet, showing a marked similarity in 
structure to bands of the third positive group of 
CO. We have observed the same bands, and 
sixteen others of the same type. This report deals 
with the vibrational analysis of the observed 
bands. 


EXPERIMENTAL 


Preliminary experiments were made with gas 
generated by the method of Germann and 
Booth.* The generator was sealed to an elec- 
trodeless discharge tube and a pumping system. 
A Féry type quartz spectrograph was used to 
photograph the spectrum, and many bands were 
obtained, including those of CO, BO, SiF, OH, 
Ne, and some which could not be fitted into any 
known band system. Visual observation with a 
spectroscope showed the presence of the bands 
of Tawde and Johnson.* Upon seeking conditions 
for the maximum intensity of the latter, it was 
found that the presence of a liquid-air trap in the 
pumping system weakened the intensity of these 
bands very rapidly. Lowering the liquid-air bath 
on the trap resulted in releasing BF; back into 
the discharge tube. In a very short time, how- 

1R. C. Johnson and H. G. Jenkins, Proc. Roy. Soc. 
A116, 327 (1927). 

2N. R. Tawde and R. C. Johnson, Phil. Mag. 13, 501 
mit "M. Strong, Thesis (M. S., Ohio State University, 
; cE B. Dull, Phys. Rev. 47, 458 (1935). 


5 A. F. O. Germann and H, S. Booth, J. Phys. Chem. 30, 
369 (1926). 
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ever, the bands of Tawde and Johnson again 
disappeared, and CO bands became very intense. 
This always occurred with a tube that was not 
baked out, suggesting that the CO was released 
from the walls of the tube, while the BF or BF; 
apparently was ‘“‘cleaned up”’ under the influence 
of the discharge. 

Through the kindness of Professor H. S. Booth, 
and Mr. Karl Willson at Western Reserve Uni- 
versity, we obtained some samples of BF; gas 
which they had dried and purified very carefully. 
The technique adopted in transporting and 
exciting the gas was as follows. A discharge tube 
of about 30 cm length and 5 cm diameter was 
fitted with a window transparent to the ultra- 
violet, either of quartz attached by a quartz-to- 
Pyrex graded seal, or of glass blown into a very 
thin re-entrant bubble. The other end of the 
discharge tube was provided with a side arm, and 
an exhaust tube. A short section of the exhaust 
tube was constricted, and contained a pellet of 
Wood’s metal. The tube was baked for several 
hours at 500°C while exhausting, and then sealed 
off from the pumps. The Wood’s metal was 
melted into the constriction in the exhaust tube, 
closing the latter so that it could be attached to 
the BF; generator without readmitting air to the 
discharge tube. The Wood’s metal was then 
melted out of the constriction, and pure BF; gas 
was admitted at a pressure of about half an 
atmosphere. Thereupon the tube was sealed off 
once more, and was ready to be transported to 
this laboratory. The precautions described were 
found to be essential in order to eliminate 
impurities, of which the most persistent was CO. 
With the technique described, the only remaining 
bands appearing besides those of BF were the 
SiF bands. Undoubtedly, some SiF, is formed by 
reaction of BF; with the hot glass when sealing 
off, as well as with the walls of the tube duringthe 
discharge. 

The discharge tube was placed within the helix 
of a high-frequency oscillator, and the pressure 
of the gas within the tube was brought down to a 
suitable value by refrigerating the side arm with 
liquid air. It was found that with a Dewar flask 
placed to bring the liquid-air level about one cm 
below the tip of the arm, the pressure could be 
maintained indefinitely by occasional adjust- 
ments as the liquid air evaporated. In this way, 
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the vapor pressure in the discharge was sub- 
stantially constant in spite of ‘‘clean-up”’ pro- 
duced by the discharge.® 

bythe 100-watt oscillator previously described’ 
was used for some of the preliminary work. The 
three-meter oscillator of the Central Scientific 
Company was also used at moderate dispersion. 
For exposures at higher dispersion, a more 
powerful oscillator was constructed. This con- 
sisted of an 852 triode in a Hartley circuit, the 
output of which was condenser coupled to a 
tuned circuit using an 849 triode biased for class 
C amplification. The plates of the tubes were 
supplied with direct current at 1200 volts by 
means of a motor-generator. A frequency of 
approximately 3000 kilocycles per second was 
obtained with the circuit elements employed. 

A Hilger E,; Littrow mounted quartz spectro- 
graph was used to photograph the region from 
4000 to 2300A. This instrument has a dispersion 
of 14A/mm at 4200A, and 5A/mm at 2300A. 
Suitable exposures were obtained in 1} hours. 
Longer exposures brought out more general 
blackening, rather than more band features. 

A 21-foot grating was also used, the region 
6000 to 2500A being covered in the first order. 


~ 6HLM, Strong and H. P. Knauss, Phys. Rev. 47, 642A 


(1935). 
7H. P. Knauss and J. C. Cotton, Phys. Rev. 38, 1190 


(1931). 
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The grating has 30,000 lines to the inch, giving a 
dispersion of 1.2A/mm, and is mounted on a 
Paschen-Runge mounting. The maximum ex- 
posure time attempted was 36 hours. 

The spectrum of the iron arc was superimposed 
on the band spectrum for comparison. A Geneva 
Society comparator was used to make wave- 
length measurements. 


RESULTS 


Nineteen bands, degraded to the ultraviolet, 
which were observed in the region 4000 to 2300A 
on plates taken with the Hilger E; spectrograph 
have been ascribed to BF. Of these, six were 
sufficiently intense to appear on the 21-foot 
grating plates in the first order, including three 
bands reported by Dull.* 

The probable error of measurements made on 
the 21-foot grating plates is 0.05 cm™+. The other 
measurements, on Hilger E, plates, have a 
probable error varying from 1 to 6 cm™ in going 
from higher to lower frequencies in the region 
occupied by the observed bands. 

In Fig. 1 the band at 3122A is reproduced 
together with a comparable band of CO at 
2833A. The scale of enlargement for the latter is 
only about 0.6 as large as that for the former, 
chosen in order to match similar features of the 
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Fic. 1. Above. (0,0) Band of the A system of BF at 3122A, with frequencies of heads shown 


in cm~!, Below 


(0,0) Band of the third positive system of CO at 2833A. Designation of heads 


from rotational analysis of Dieke and Mauchly.* The lower strip is enlarged about 0.6 as much 
as the upper one in order to illustrate similarities of band structure. 
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two bands. The designations for the sub-heads 
of the CO band are taken from the rotational 
analysis of Dieke and Mauchly.* The features of 
the BF band are not sufficiently clear-cut for 
purposes of rotational analysis. Three strong 
subheads could be distinguished in most of the 
bands, corresponding to the P3, Pz and P, heads 
of CO. Weaker heads, observable in some of the 
bands, appear to correspond to O and Q branch 
heads of the CO bands, and have been designated 
accordingly. The edge of the Q branch, while not 
strictly a head, is listed wherever it appears as a 
distinct feature of the band. 

Twelve of the observed bands belong to the 
same system, designated as the A system. One 
band at 3560A is the isotope band of B'® F’, 
corresponding with the (0,3) band of B" F'® at 
3548A. The intensities and frequencies of the 
heads of the A system bands are listed in Table I. 
Three other bands have final levels in common 
with the A system, but an initial level which is 
different ; these are designated as members of a B 
system, and are tabulated in Table II. In Table 
III are listed three bands not assigned to either 
system. 

The vibrational spacing of bands of the A 
system is represented by the equation: 


vp, = 31,869.98 + 1631.40(0' + 3) — 23.10(v’ +3) 
— 1323.64(v’’+3)+9.40(0’ +3)? 


where », represents the frequency in cm™ of the 
sub-head assumed to be nearest the band origin, 
the one corresponding to the P; sub-head of the 
CO band, Fig. 1. The frequency separation of the 
sub-heads varies from one band to the next, and 
in the absence of a rotational analysis of the 
bands, no attempt is made to show this in the 
equation. The B system bands are represented 
by: 


vy, = 38,651.4 — 1323.64(0” +4) +9.40(0" +3)? 


where », again refers to the P; head. 

The constants in the above equations are based 
on the high dispersion plates. The residuals (obs.- 
calc.) thus obtained for the heads measured at 
low dispersion are preponderantly positive, 
which may indicate a slight displacement of the 


iron comparison lines relative to the bands. 


*G. H. Dieke and J. W. Mauchly, Phys. Rev. 43, 12 
(1933). 
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rasB_e |, A system bands of BF. 


Oss Ons 


Cate Cal 
’ Int. » Ons. (cm cm INT. » Ons ! 
9 32058.93 0, 3 33640.7 O 
9 *32038.70 P, 4 *33623.59 P 
00 9 *32029.67 P, 1 0 4 *33613.86P 
7 *32020.49 P; +0.06 4 *33605.69 P +0.06 
2 *31995.18 O 2 33580.0 O 
5 30751.6 0, 0 32334.6 P 
7 *30734.70 P, 0 32320.1 Py» 
01 5 *30724.80 P, 110 32301.7 P +0.9 
5 *30715.41 P; —0.18 1 32287.2 O 
4 30696.82 O 
5 29467.5 Q, 1 31052.1 P, 
6 *29449.30 P, 1 31029.6 Py, 
02 5 *29438.67 P» 121 310195 P; +4.7 
3 *29429.70 P; +0.15 
4 29413.2 O 
4 28199.2 Q, 3 29782.1 P 
5 *28182.66 P, 4 29765.4 P» 
03 5 *28172.23 Pe 13 2 29756.2 P 8.7 
4 *28162.20 P; —0.11 
3 *28144.07 O 
0 3**0 28087.5 0, 2 28534.8 0, 
0 28072.1 P, 1 *28518.60 P, 
1 28057.4 Py» 0 *28509.07 P. 
0 28039.3 P; —4.3 |1 4 0 *28499.36P; +0.29 
1 28482.3 O 
0 26952.4 QO, 0 27305.6 P 
0 26938.4 P, 0 272904 P 
04 O 26922.3 P; +84 |1 5 0 272780 P 8.6 
0 26898.4 O 0 27255.9 O 


P; +06 


* Frequencies computed from wave-length measurements made on 


21-foot grating plates are marked by an asterisk. Other wave-length 
measurements and all intensity estimates are based on F; plates 


** Isotope band, BF, 


Since the residuals are not large in comparison 
with the probable error, no adjustment of the 
data seemed necessary. 


DISCUSSION 


With the aid of constants of the A system of 
the B'" F'® molecule, positions were calculated 
for possible bands of what we may call the A 
system of the less abundant molecule B’® F’. 
Only one measured band could be assigned to the 
A‘ system. In all other cases, the failure to 
observe the A‘ bands could be accounted for 
either by the low intensity to be expected or by 
confusion with an A band. The excellent agree- 
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TABLE II. B system bands of BF. 


IN vy OBS CALe 

6 38028.3 O 

7 38009.7 P, 

0 0 8 37999.1 P» 
7 37991.4P —0.5 

6 36722.1 QO; 

4 36704.7 P 

0 ] 6 36700.5 P» 
4 36687.9 P +0.8 

2 36660.3 O 

4 35435.8 O 

3 35418.4 P; 

0 2 4 35409.4 P» 
3 35400.8 P; —0.3 

2 35390.7 O 


ment between the observed and calculated 
position of the isotope band removes any doubt 
which otherwise might attach to the identifi- 
cation of the emitter of the bands. 

An estimate of the interatomic distance, 7,, 
can be made with the aid of the relationship 
proposed by Allen and Longair :° 


rw! = 10.6 X 10-*, 


where yu is the reduced mass, and w, is the fre- 
quency (in wave numbers) for vibrations of small 
amplitude about the equilibrium position, ob- 
tained from the spectrum. On this basis, the 
value of r. for the common (lower) state of the 
A and B systems of BF is 1.33A; that for the 
upper state of the A system is 1.24A. 

On the basis of the predicted similarity be- 
tween BF and CO spectra, and the observed 
resemblance illustrated in Fig. 1, it is natural to 


°H. S. Allen and A. K. Longair, Phil. Mag. 19, 1032 
(1935). 
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TABLE III. Unassigned dands. 


IN vy Oss. INT vy OBS INT vy OBS 
0 38833.1 | 0 38544.5 | 0 37224.8 
1 38801.4 0 38518.3 | 0 37207.1 
2 38787.3 | 1 38505.4 | 1 37196.0 
1 38775.1 | 0 38493.7 | 0 37187.1 
0 37118.6 


rasLe [V. Comparison of lower level of BF bands with aT 


level o7 CU. 


we(H vewe(H) D re I 

cm (cm volts i t 
CO 1739.3 14.5 4.0 1.204 1.686 
BF 1323.64 9.40 ? 1.33* 1.36" 


* Computed from empirical relationship of Allen and Longair.’ 


suppose that the A system bands correspond to 
the third positive bands of CO, and the B system 
to the 3A bands of CO. If this is correct, the 
lower level of the bands is itself an excited state, 
and the molecular constants obtained for it do 
not apply to the unexcited state of the molecule. 
As suggested by Dull,‘ the lower state may be 
either *II or *= depending on whether the trans- 
ition is *II—*> or *Y —II. In Table IV a compari- 
son is made with the a ‘II level of CO, on the basis 
of the assumed correspondence. 
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rhe empirical band spectrum rho-values, p*, as com 


puted from the constants of LiH and LiD depart markedly 
from the ratio of the reduced masses, p=(u/u')!. The 


isotope discrepancy, (p*—p) is positive for the 11 band 


constants measured for the ground state and negative for 
10 of the 13 in the upper state. By using Dunham's 
equations, the anharmonicity constants a), dy-+-ds are 


computed for the ground state and the effect of anhar- 


monicity on the isotope ratio calculated. In only 2 cases 


out of 5 does this reduce the isotope discrepancy and then 
by amounts from 6 to 22 percent. In the upper state the 


anharmonicity constants must be got by a laborious 


INTRODUCTION 


HE 'S—'Z band spectra of LiH and LiD are 
of particular interest for testing the finer 
points of the theory of molecular spectra for two 
reasons. In the first place the molecules are 
relatively light and therefore the isotope shifts to 
be expected are larger than for any known 
molecules except H, and HD. In the second place 
the upper electronic state is quite abnormal in 
that the molecules appear to shrink in size as well 
as to vibrate more rapidly as the amount of 
vibrational energy increases. Similar effects are 
present in the case of the hydrides of Na and K 
though they are there less pronounced in magni- 
tude than in the present case. In addition the 
simplicity of the structure of 'S—'Y bands is of 
great assistance in the measurement of those 
molecular constants on the precise determination 
of which the present discussion so largely depends. 
Preliminary accounts of results and the data 
on which the present paper depends have been 
given earlier.' 

It is now well known experimentally that 
the molecular constants of hydrides and deu- 
terides (JH and MD) are not in general re- 
lated precisely by powers of the reduced mass 


‘Crawford and Jorgensen, Phys. Rev. 45, 737 (1934); 
Part I, Phys. Rev. 47, 358 (1935); Part II, Phys. Rev. 47, 
932 (1935); Phys. Rev. 48, 475 (1935). In Parts I and II 
complete references to the earlier literature are given. 


approximation process. It is shown that the resulting set 
of a values accounts quite satisfactorily for both the sign 
and value of the 9 band constants, which have been 
previously reported as abnormal. This well-known ab 
normality can therefore not be ascribed in any significant 
measure to L uncoupling. The Dunham corrections further 
reduce the isotope discrepancy in 4 out of 5 cases by an 
average of about 25 percent. The dissociation energies 
(Morse) are found to be for LiD, Do” =2.53, Dy’ =1.15 
volts and for LiH Do’ =2.54 and Do’=1.18 volts (all 
+0.2 volt and measured from the v=0 level). 


ratio p=(u/u')?, where p= MH/(M+H) and 
u'= MD/(M+D), of the simple isotope theory. 
As here defined p = 0.750486 and is therefore <1 
and we find in the ground state that the eleven 
constants which can be measured with reasonable 
certainty all give ratios (say p*) which are 
slightly nearer unity. 

Thus (p*—p) which we may call the isotope 
discrepancy is in every case positive. In the first 
excited state, on the other hand, we find (p* —p,) is 
positive for only three of the constants and 
negative for the other ten. It seems probable 
although sufficient data are not at hand for 
certainty—that the first result, i.e., (e*—p) >0, 
is the normal behavior for the ground states of 
most molecules. 

As in Part II we shall write the complete 
energy expression (in cm™') for a _ rotating- 
vibrating diatomic molecule in a singlet state as 
W W, 

—- © £ 
he he t=0,1-++ j=0, 1--- 

x Yi,;*(v+3) LI(J+1)]! (1) 


The empirically determined “rho-value,”’ »p,;* 
referred to in the paragraph above, is then given 
by (¥,;'*/Y,,*)“!/"*2 and in general varies with 


both | and j. Since the elementary isotope theory* 


>See, for example, Jevons, Report on the Spectra of 
Diatomic Molecules, p. 212. 
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predicts that (p,;*—p)=0 the departure of this 
difference from zero may be taken as an indi- 
cation that the empirical molecular band con- 
stants Y,;"* and Y,;* have lost their (assumed) 
precise mechanical significance. This implies for 
example that Y,9* (identical with the vibration 
constant w,*, where the asterisk is used to 
emphasize its empirical nature) does not give 
exactly the rate of vibration for small amplitude 
motions about the equilibrium position. This 
latter mechanically simple constant we shall 
designate by w,. Likewise Yo,\*=B,*#B,, where 
B,=h/8r'cI,, and therefore the actual moment 
of inertia at the equilibrium position, J,, cannot 
be computed directly from h/87*cB,*. 

Now although various approximations in the 
solution of the wave equation for a rotating- 
vibrating diatomic molecule have been pointed 
out as partly responsible for these differences, it 
seems desirable to call attention to the classifi- 
cation of neglected effects given by Van Vleck® 
in a recent careful reexamination of the whole 
problem. The neglected energy terms, although 
not all capable of specific evaluation in forms 
suitable for immediate comparison with experi- 
ment can be considered as due to four effects: 


1. Anharmonicity terms due to the interaction 
of rotation and vibration. This question has been 
considered in great detail by Dunham‘ who, 
writing the molecular potential energy function 
in the general form: 

V =aof[ 1 +a, + ack +a 3é+a,t'+---] (2) 
obtained small correction terms to the usual 
energy expression. These correction terms (actu- 
ally only five were evaluated, see Eqs. (7) to (11) 
below) involve the potential coefficients a), de, 
etc., and a factor ~p?=(B,/w.)? and in general 
may be positive or negative depending on the 
precise form of Eq. (2). Since their values depend 
on p*, these terms will be very small for heavy 
molecules and of importance only for hydrides. 
Here p* varies from ~10~‘ for certain states of 
Ho, LiH, etc., to 10-* and 10~* for a large number 


3Van Vleck, Phys. Rev. 49, 417A (1936); J]. Chem. 
Phys. (in press). Note that the energy expressions used 
by Van Vleck are in absolute energy units while ours are 
in cem7!, 

4 Dunham, Phys. Rev. 41, 721 (1932). 
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of the electronic states of molecules formed from 
the first two rows of the periodic table. 

2. Terms due to the interaction of electronic 
and nuclear motions (neglecting the interaction 
of neighboring electronic levels). The use of a 
potential function of the type of Eq. (2) is 
equivalent to considering the frequencies of 
electronic motions so high that the nuclei in a 
nonvibrating molecule are essentially fixed. To 
this approximation, terms in the original com- 
plete wave equation for a diatomic molecule® 
which vanish when averaged over the frequencies 
connected with the electronic motions, are neg- 
lected. In higher order approximations of the 
energy of the system these neglected terms do, 
however, make contributions. The form of these 
contributions has been considered by Kronig,*® 
Dieke’ and Van Vleck* where simple approxi- 
mations for a few of these terms are given, 
suitable for electronic states of symmetrical 
molecules with L=0 or for those states (with 
I #0) which can be satisfactorily approximated 
by considering the nuclei coalesced. The precise 
evaluation of these terms in the general case 
requires more detailed knowledge of the molecu- 
lar wave-functions concerned. They do, however, 
contain p? as a factor and probably tend to 
increase both (py*—p) and (po:\*—p). The re- 
maining corrections arise from the interaction of 
other electronic states with the one under con- 
sideration and are of two sorts. 

3. L uncoupling terms.* These are due to the 
interaction of electronic states whose A values 
differ by +1, are zero for o'S states and have no 
appreciable effect on vibrational constants, i.e., 
the coefficients Y,9*. Their sign is + or — ac- 
cording as the perturbing state is above or below 
the one in question. They vary as },/v, where 
v=the electronic separation of the states 
involved. 

4. Terms arising from the interaction of states 
of the same A value. These terms are of a 
distinctly new type and do not affect the rota- 
tional constants Y»);*. They also are of the order 
of B,/v and may be + or —. 


5 Van Vleck, Phys. Rev. 33, 467 (1929), Eq. (11). 
6 Kronig, Physica 1, 618 (1934). 

7 Dieke, Phys. Rev. 47, 661 (1935). 

5 Hill and Van Vleck, Phys. Rev. 32, 269 (1928). 
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A detailed evaluation of the effects (3) and (4) 
requires of course a fairly complete knowledge of 
the electronic states of the molecule in question. 
Since the corrections of these four types (unless 
an actual large scale “‘perturbation”’ is present) 
are usually small, we can treat them as inde- 
pendently additive. Consequently we can write: 


Y,;*= Yi 1+ (Ad ij+ (Ae) + (Adit (Adis) (3) 


where Y, is the mechanically significant value 
Y,,* would assume if the contributions (4,);;, 
etc., of the above effects were zero or experi- 


mentally undetectable. Here then 


Y,;**/ Y1;* = (o*) **/, (4) 


and since p is usually near unity we may write 
to the same degree of approximation: 


(p1;* —e) =Arpijt+ Acpij t+ A391; +Api;, (5) 


where in general Aypyo and A4po; will be zero. 

From the above it is apparent that only the 
corrections due to anharmonicity [(A;),; and 
Aip.; | can be computed in a detailed manner at 
the present time. We therefore propose in this 
paper to deal primarily with these. 


EFFECTS OF VIBRATION AND ROTATION 
INTERACTION 


Dunham's calculations of the explicit form of 
the Y,,’s in terms of the anharmonicity constants 
a, to ag of Eq. (2) indicate that each coefficient of 
Eq. (1) is of the form: 


Y,*=fij(w., B.)[Aifai, ae-+-) 
+ pC) ;(a), de-+-) ] (6) 


and that therefore by neglecting the last three 
terms of Eq. (3), Yij=fij/Ai; and (A));, 
= p°C);/A1;. Here A,; and C,; are functions of the 
anharmonicity constants d@;, d2:-- and are there- 
fore the same for each isotope, while f,; and p 
both contain factors dependent on mass and are 
different in the two cases. In general f,; decreases 
rapidly with both / and 7. Dunham's calculations 
were carried far enough to obtain five of the C); 
terms. Since it will be necessary to refer re- 
peatedly to these equations we shall give them 
here in a form slightly more suitable for nu- 
merical reckoning.** 


_** For convenience fractions have been divided out and 
given uniformly to four decimals. This is quite accurate 
enough for any experimental data. 
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Vio* =a.* =w.[1+97Ci0 ], (7) 
Yo.*=B,*=B[1+p'Co], (8) 
Yo.* =D,.* = —4B.p'[1+p*Coe]), (9) 
Y),*=a,*=B,p(6(1+4,)+7°Ci1], (10) 
Yoo* = x,* =B,/2[3(a2—1.25a,7) + p?Coo |, (11) 


Yoo* = B,/8[3a2—1.75a,°], (12) 


_ 

. 

— 
\| 


Xo= B.p 2[ 10a. om 35a 1d 
2 (56.25a;* — 8.5a2)de om 22.031 3a,‘ ], ( 13) 


VY so* = x2 =5B.p*(3.5ag—15.75aja2 
+ (33.9375a;* —8.125a2)a,4 
— (61.0313a,3 — 60.375a;d2+7.875a3)a3 
+ (90.8791a;4—77.95314a;"ae 
+4.6875as")a2—22.6084a,;°], (14) 


Y2\* =a2=6B,p*{ 5a3— (3+ 13a,)dae 
+(10+7.5a,+7.5a)ai:+5], (15) 


Y3:*=a;3=20B,p*[ 7a; — (4.5+25.5a1)a4 
+ (14+ 26.254; —22.5a2+54.375a;")as 
— (8.5+ 45a; —6.375a2+59.0625a;’ 
— 51.375a,;a2+94.3125a;*)ae 
+ (214+ 28.125a,+35.25a;’ 
+ 29.7422a,°+29.7422a;*)ai+7 ], (16) 


Y..*=6,= 12B.p*[4a2.— (9+4.5a)a,—9.5 ], (17) 


Yoo = B= 24B,.p*[ 15a, — (30+47.5a;)az 
+ (61+117a,+103.5a,*—26a2)de 
— (125+123.75a,+90a,*+45a,*)a,—65], (18) 


Yo;* = 7 .* = 16B,.p*(a,+3), (19) 


Yi3* = 71 = 12B.p*[ 26.6667a3— (120+ 88a,)a2 
+ (279+ 1894a,+63a,*)a, +233], (20) 


Yo,* = 64B..p*[ ae _ (9+ 2.25a l ja yom 13], (21 ) 
where the C;;’s are given by 


Cy = 4[25a;—47.5a,a3+ (57.3750; 
—16.75a2)a2—18.0469a;*], (22) 


Co) — 4[ 15a3— (9+ 23a,)de 
+(14+10.5a,+10.5a,2)a,+15], (23) 


Coo = 3[ —45a4+ (90+ 102.5a,)as3 
— (119+ 207a,+ 166.5a," —46a2)a» 
+ (199+ 173.25a,+126a," 
+-62a,*)a,+ 163], (24) 
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Cy =([175a;— (112.54+397.5)as 
+ (190+ 356.254; — 357.542 


+643.125a;*)a;:—(1 
— 125.625a2+.602.4375a;? —584.625a,a2 
+891.1875a;*)ao+ (285+ 286.8754, 
+ 288.75a;°+ 243.6328a;' 
+ 243.6328a;*)a,+175], (25) 


67.5+459a, 
3 


Coo= + 256a¢ — 682.5a,a;+ (1066.875a;° 

—442.5a2)a4—(1491.5625a, 

— 1833.75a,a2+271.25a3)a3 

+ (1873.3203a,4— 1940.1563a;"ao 
+213.375a9")a2—408.3106a 1° | (26) 


From the equations it is at once apparent that 
if all of the Y*’s could be determined with 
comparable precision most of the a’s could be 
determined in two or more ways. Actually Yoo* is 
not measurable directly, Yos* is seldom obtained 
and only an extremely complete and careful 
analysis gives the a’s, 6’s or y’s with enough 
accuracy to be of much value. The detailed 
procedure to be employed therefore depends on 
the data available and in any event the uncer- 
tainty in a given a increases with its subscript. 


GROUND STATE 


In the ground state of both LiH and LiD the 
quantities w,, B,, D, and H7/, change so regularly 
with v that relatively good values of most of the 
smaller constants are obtainable. It is desirable 
therefore to compute the a’s from the smaller 
constants and to use these for determining 
independently the correction factors C); of Eqs. 
(22) to (26). This is open to the objection that the 
corresponding correction terms are not known 
for Eqs. (13) to (21) and may be relatively 
important. Since their contributions to the total 
energy expression will be small if not experi- 
mentally undetectable this means merely that to 
represent data of a given precision we need the 
a’s very much less accurately. 

A glance at Eqs. (13) to (21) shows that a; can 
be obtained from (19) and thence az from (17). 
Then a; can be obtained from either (15) or (20) 


and a, from (13) or (18). Finally a; and ag can be 


obtained from (16) and (14), respectively. The 
results of these calculations are given in Table I 
for both LiH and LiD where because of the 
larger number of bands observed in the latter 
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TABLE I. Potential coefficients for ground state. 


Fq Li?H Li?D 
(19) a — 1.84 — 1.858 
(17) a + 2.33 +2.580 
(15) ) — 2.69 — 3.446 
(20) as —2.98+ — 3.456 
(13) +5.51T +5.79+ 
(18) aes +4.35 +6.53 
(16) a — 14.0 —12.2 

(14) a ~ +33 + 20.7 

t Values so marked are the more reliable of a pair. 


case the results for LiD are more trustworthy. 
When a given a is obtained in two ways, one 
usually can be recognized as the more reliable. 
Thus x2 contributes 26 percent of the final value 
5.79 of ay in Eq. (13) while 82 contributes 23 
percent towards 6.53 in Eq. (18). Since x» is 
probably more certain than the very small 82 the 
first value is the more dependable. Such values 
are indicated by a dagger. It should be further 
remarked that since a3(= Y2:*) was too small to 
be measured for either isotope, theoretically a 
could not be determined. Actually, however, 
solution of Eq. (16) for a; shows that unless | a; 
exceeds considerably its rough theoretical value 


of 1.3 percent |a2|, a reasonably good value 
should result from setting a;=0. Likewise in the 
case of ag determined from Eq. (14) we require 
x3(= Y4o*) which could be measured only for 
LiD. Since this occurs only in the practically 
dimensionless ratio x«;/5B,p? in the solution for 
ds, the LiD value of this term was employed in 
finding ag for LiH. The coefficients are seen to 
alternate regularly in sign and to increase almost 
linearly in absolute value with the subscript 
This as Dr. H. M. James, in this laboratory, has 
found, is apparently a characteristic of normal 
electronic levels. 

With a reasonably good set of a values we are 
now in a position to do several things. We may 
first determine from Eqs. (7) and (8) better 
approximations to the mechanically significant 
constants w, and B,. These then enable us to 
determine what we might call the Dunham 
approximation to D,*, x,* and a,*, and finally to 
determine A,p,;; the Dunham contributions to the 
isotope discrepancy. The first step gives us 


Vio* =w,'[1+1/10,042], Yiot=o.[1+1/5660], 
Yo:'* =B,*(1—1/7376], Yo*=B.[1—1/4157]. 
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TABLE II. 
VALUES OF ¥;;** (ror LiD) VALUES OF Y;;* (ror LiH) 
C1 Observe Dunham Kratzer Pekeris Observed Dunham Kratzer Pekeris 
02 — 35.69 — 2.856-107* —2.72510-* —2.727 —8.62-10-* —8.58510-* —8.594-10-* - 
11 — 235.6 0.09198 — (0.0875 —0.1060 —0.0780) —0.2132 — 0.2070 — 0.2480 —0.1823 
20 — 837 — 13.23 — 11.046 12.70 — 23.20 — 19.64 22.54 - 
TABLE III. Ground state isotope effects. 
BAND 
CONSTANT l 1). Cy; p1;* pi; p aT Sis pas” 6 
w,* 10 1 +6.18 0.750628 +0.00014, — 0.00006 —42Q% 
B.° 01 1 —8.42 .750685 + .00020 + 00004 +22% 
D.* 02 1 — 35.69 .7520 + 0015 + 000084 + 6% 
a;* 11 6(1+<a;) <1 — 235.6 .7556 + .0051 — 000143 —3Y 
x,* 20 3(a2—1.25a,?) <1 — 837 7551 + .0046 — .000754 —16% 
sign in the last column indicates that the Dunham correction decreases the isotope discrepancy 
TABLE IV. Upper state isotope effects. 
BAND 
CONSTANT l 1); Pj pij* e1;"—P S10); S19); p14" p 
w,* 10 l — 53.25 0.770918 +0.02043 +0.002618 +13% 
B.* 01 1 —46.22 754561 + .004075 + 001136 +28% 
dD 02 1 +214.4 .7424 — .008 — 00264 +30% 
a* 11 +-0.83538 +3550 .2036 — .1622 — .06963 +43% 
x* 20 + 18.064 — 8840 6951. — .0553, + .01203 —22% 


These results give p*=2.859X10° and p 
= 1.611 X10~° and enable the observed values of 
D.*x,* and a; to be compared with the Dunham 
approximations as well as those of Kratzer® and 
in the case of a,*, of Pekeris.'"° The numerical 
values are collected in Table II, where it appears 
that the difference between Kratzer’s and 
Dunham's values for D,* are inappreciable and 
fairly close to the observations. In the case of 
a;* Dunham’s values are the best, though 
differing by more than the experimental uncer- 
tainty, with Kratzer’s too high and Pekeris’ as 
much too low. The agreement in the case of x;* is 
not as good and Kratzer’s expression in this case 
lies closer to the observations. 


IsOTOPE RATIOS 


Turning to the isotope relations a little manipu- 
lation of Eq. (6) shows that we can write 


Aipij= —[o(p?— p”)/(1+2/)Ay JC, (27) 


® Kratzer, Zeits. f. Physik 3, 289 (1920). 

1° Pekeris, Phys. Rev. 45, 98 (1934), also Jevons, re- 
ference 1, p. 27, and Birge, Bull. Nat. Res. Council (1926), 
p. 237. 


where p=0.750486. We find also that pj9* 
=().750628 and po;*=0.750685. Hence (pi9*—p) 
= +(0.00014. and (p»;*—p) = + 0.00020. But from 
Eq. (27) we have Ajpio= —0.00006 and Apo; 
= +0.00004. Therefore the residual discrepancy 
(p1;*—p—Asp.;), in the first case is increased 
to +0.00020 and in the second decreased to 
+0.00016. Since Z uncoupling is ruled out for the 
ground o'S state the first discrepancy is pre- 
sumably due to effects (2) and (4) while the 
second must be due to effect (2) (the Kronig 
effect) alone, since effect (4) does not influence 
rotational constants (and particularly B,*) ap- 
preciably. Similar calculations have been carried 


* and the results are 


out for poo*, pii* and poo 
collected in Table III where it is seen that the 
Dunham effect is inappreciable in the case of a,* 
and varies from 16 to 42 percent of the total 
discrepancy in the other cases. In only two cases 
does it decrease the discrepancy. 


EXCITED STATE (CALCULATION OF @ VALUES) 


The attempt to determine the @ values for the 
excited state in the same way as for the ground 
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state immediately leads to quite discrepant values. 
Thus Eq. (19) gives for a; roughly —2 whereas 
Eq. (10) (neglecting ~°C;,) gives —0.814. Now 
the discrepancy may be due to the uncertainty in 
HT .* as measured, or to the neglect of p?Co3 in Eq. 
(19), of p?Ci, in Eq. (10) as well as of effects (2), 
(3) and (4) above. Whatever the difficulty, a 
reasonable value for a; cannot be obtained by 
averaging such divergent values—and since the 
value of a2 depends to such an extent on the 
value of a, using Eqs. (11) and (17) we have four 
widely different values for a2 with no basis for 
preferring any particular one of them. Asa result, 
beyond a2 even the sign of the higher a’s is in 
doubt. After following down all the possibilities 
it soon became probable that the neglect of the 
Dunham correction terms was not the primary 
difficulty—but that the real trouble arose from 
the difficulty in extrapolating the measured 
values of D, and H,, respectively, to obtain D,* 
and H,*. This uncertainty was such as to 
produce even greater errors in 61, 71, 82, etc., and 
consequently to render the values of a, a2, a3 and 
a, as obtained respectively from Eqs. (19), (17), 
(15) and (13) quite untrustworthy. It seemed 
accordingly necessary to resort to a method of 
successive approximation, using only those ex- 
perimental constants which could be determined 
fairly reliably from least squares representation 
of the remaining data, viz., the values of w, and 
B,. Since the variation with vibration of w, and 
B, was much smaller for LiD than for LiH, the 
calculations were carried out for the LiD data 
(actually the LiH data gave comparable results— 
though the actual numerical values depended 
uncomfortably on the precise number of obser- 
vations used in the least-square reckonings). The 
necessary expressions for LiD are: 


w, = 180.711 + 13.9874 — 1.1784? +0.07932u' 
— 0.00353u4+0.00006940u° (7 points) 


and 


B, = 1.6048, +0.015940u —0.0023732u? 
+0.00012783u* —0.00000417u* (6 points). 


An examination of the equations from (10) to 
(17) then shows that the a’s of odd subscript can 
be obtained from the rotational constants a, a, 
and a; while the a’s of even subscript can be 
obtained from the vibrationai constants x), x2 
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and x;. The general procedure was to determine 
a first approximation to a), say @\, from Eq. (10) 
by neglecting ~?Cy. With this a first approxi- 
mation to de, say d2;, was then obtained from 
Eq. (11) neglecting »*C,0. Then by using Eqs. 
(15), (13), (16) and (17) in order, @3;, a4), @5; and 
as: were obtained. It was then possible to 
compute C,, and Coo approximately to get the 
second approximation of a; and de, say a2 and dep, 
respectively. This led then to better values for the 
higher a’s, to more reliable values for C1; and Coo, 
etc. The process is unfortunately very tedious 
and while not converging very rapidly soon gives 
results which do not oscillate very much in 
successive stages of the calculations. 

The results are given in Table V for ay. and 
a,; together with the values of C); resulting." 


TABLE V. Potential coefficients for upper state. 


BAND VALUE OF 
Eo. CONSTANT 1 A see in Ch 
10 Y,,*=a;* a, — 0.8604 —0.8608 Cy — 53.45 
11 Voo*=x,* ae + 6.884 +6.948 Cp, — 46.22 
15 Voa*=a,. a; —11.33 —11.44 Cor + 214.40 
13 Vsot=x. ay +30.40 +31.11 Cc + 3550 
16 Y3,;* =a; a, —66.57 —67.77 Ca — 8840 
17 Y40* =X3 as +67.32 +-68.48 


The striking fact is the rapidity with which the 
coefficients increase in absolute magnitude as 
compared with the more normal ground state, 
and in particular the fact that a2 is over 12 times 
as large as |a,|. This will be discussed in more 
detail below. 

The next step was to determine whether this 
set of a values was satisfactorily consistent with 
the predictions of the remaining equations. For 
this purpose Eq. (7) was used to determine from 
Y;o*=wo* the true value, w,, of the rate of oscil- 
lation for zero amplitude, and from Ed. (8), 
Y».*=B,*, the true value of 4/82*/. These then 
gave p’=(B,/w,)* and from Eqs. (9) and (19) 
theoretical values of Yoo2*=D,.* and Yo3;*=I/,*, 
respectively, could be obtained. The results 
were Yo2*(calc.) = —5.162 10-4, Yo3**(calc.) = 
+34.2210-§, Yoo*(calc.) = —16.88 10-4 and 
Yo3*(calc.) = 202.18 K 10-°. It was next necessary 
in each case to obtain least-square curves which 


The values of C,, and Cx are the fourth approxima- 
tions since these were the values determining the correc- 
tions in Eqs. (10) and (11) for the fifth approximations to 
a, and a», respectively. 
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passed through the observed data and these 
calculated values for v= — }. Since measurements 
could not be made for v=1 on LiD (and were 
uncertain on LiH) it was necessary to assign a 
probable value to say D,, compute a least- 
square expression for D,; to D; and calculate a D,* 
by setting v= —}. After a few trials a value of 
D, could be chosen such that this extrapolated 


D.* coincided with the one calculated theo- 
retically. This gave for LiD 
B,;'=+2.64x 10-4, B2'= —0.8907 X 10-5, 


and B.=—4.7353x10-°. 


6,=+11.23x10-4 


A similar procedure in the case of the 77, values 
gave for LiD y:'= —39.54X10-%, yo'= +18.420 
and for LiH, y,;=—250.0110° and ye= 
+122.44x 10-*. Now it is to be admitted at once 
that great reliance cannot be placed on the actual 
values obtained from such an attempt to repre- 
sent by least squares functions changing as 
rapidly as are D, and J/,. The results, however, 
do represent a considerable improvement over 
those obtained empirically in Parts I and II. 
Thus in the case of LiD using Eq. (17) and the 
old value of 8;'=1.3210~-4 we find a= +3.56 
which is about half the value given in Table V. 
If we use the new value of 8;' given above we find 
adg=+5.90 a result differing by only 15 percent. 
Likewise for LiH from the old 8,(=10.1X10~*) 
we have d2=5.30 while from the new 81, a2=5.80, 
likewise a shift in the right direction. Again 
using Eq. (20) and the values of y;‘ and y; we 
have for LiD a; (old)=+3.9, a3 (new) = —5.9 
and for LiH, a3; (old) =+3.7, a; (new) = —3.0. 
Here the new values are both too small but have 
been altered 200 percent in the right direction 
and actually changed in sign by our process. 
Finally, using Eq. (18) and the values of 62‘ and 
B. we have for LiD: a, (old)=+69.0 and 
a, (new) = + 51.2, while for LiH : a4 (old) = +56.5 
and a, (new) = +53, again a shift in each case in 
the right direction. 

We therefore conclude that to the extent that 
(1) the undetermined Dunham corrections p*C;; 
and (2) the omitted Kronig and Van Vleck terms 
may be neglected—the anomalous behavior of 
the upper states of LiD and LiH can be ac- 
counted for semi-quantitatively in terms of a 


potential curve which departs in a characteristic 
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way from the usual shape. It is thus apparent 
that the anomalies of the upper state involving 
some 9 small coefficients which differ in sign or 
magnitude from the usual case, are pictured in a 
very satisfactory way in terms of a single 
potential curve and that this same curve ac- 
fraction of the observed 


counts for a sizable 


isotope discrepancy. 
RELATIONS 


Upper STATE ISOTOPE 


We are now in a position to compute the 
Dunham Aipij, 
discrepancies just as in the case of the ground 
state. The results of this with p(p?—p*) of Eq. 
(27) =0.4917 X 10-4 are given in Table III where 
the striking fact is apparent that the isotope 
discrepancies are reduced in four out of five 
cases, the Dunham term varying from 13 to 43 
percent of the total discrepancy (p:;*—p). 


contributions, to the isotope 


POTENTIAL CURVES 


With the a values determined we may now 
write approximations to the potential curves 
(convergent near the equilibrium positions only) 
of the ground 'Y and upper 'S states of LiH and 
LiD. These are for the ground state: 


V =ay'#(1 — 1.858¢+ 2.5808? — 3.45¢' 
+5.79&4—12.2¢5+20.72+---, (28) 


where for LiH ao’ =6.57078X10* and for LiD 
dy = 6.5715 X 104 and for the upper state: 


V =ao" (1 —0.8608¢ + 6.9482 — 11.443 


+31.11¢4—67.77¢+68.48¢°+---, (29) 


where for LiH ao’ =4.91646X10* and for LiD 
ay’ =5.11154x 10%. 

The ground states for the two molecules are 
thus very closely alike and follow the usual type 
lying above the osculating parabola for §<0 and 
below for £>0. On the other hand, the curves for 
the upper states differ somewhat more from one 
another and are distinctly anomalous in that 
although they rise above the corresponding 
parabolas for §<0 they quickly cross over and 
pass above again for small positive values of &. 
This behavior is more readily seen in Fig. 1 where 
the dotted curve is given by ao’ and the heavy 
curve represents Eq. (28). It is this peculiar form 
upper which 


of anharmonicity in the state 
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Fic. 1. The heavy line represents Eq. (29) while the 
parabola is aot. The rapid rise of the curve above the 
parabola on the right is to be noted. It is this which causes 
the molecule to travel more rapidly than usual along the 
right branch and hence to have an average £, =, which is 
less than £9. The corresponding average values for v=0, 1, 
etc., are indicated by circles at the appropriate height. It 
is this peculiar shape which likewise causes the observed 
vibrational rate w, to increase with v to a maximum value 
before decreasing thereafter in the normal way. The heavy 
curve of course must eventually fall away and cross the 
parabola as the molecule goes toward dissociation, though 
that point is far out of the range of validity of Eq. (29). 


accounts in detail for the anomalies to which 
attention has been called so often in the litera- 
ture. We therefore assume that although some 
L uncoupling may be present in this state it 
probably does not have much effect on the 
potential curve and is reflected only in the 
residual isotope discrepancies (p);* —p—Ajpz)). 

For purposes of reference the values of B,* and 
w.* as well as the values corrected for the 
Dunham effect for both states of the two 
molecules are given as follows: 


‘B.'*" =1.6048 B,” =1.6107 re’ =2.5795A 
w,'*" =180.711 w.*’ =181.472 


LiD aa ' . 
B.**" =4.2338 B,'"’ =4.2344 


re’ =1.5909,A 


wei” =1055.12 — w” =1055.015 


Tr. JORGENSEN, JR. 


B.*’ =2.8186 B,’ =2.8375 re’ =2.5896A 
az |@e% = 234.41 w,’ =236.225 
Li) Be 7.5131 B." =7.5149 r. =1.5915A. 
w.*” =1405.649 w.” =1405.401 


From these the equilibrium separations are found 
where (r,’—r,”)=+0.0101A and (r,’’—r,*”’)= 
+0.00055A. 

Further from Eq. (12) we see that the height 
of the potential curves at the minima should be 
different for the two isotopes. Thus for the 
ground state Yoo7*= 15.93 cm and Yoo'*=8.98, a 
difference of 6.95 cm~', while for the upper state 
Yoo*=7.39 cm and Yoo'*=0.26 cm, a differ- 
ence of 7.13 cm~. Consequently we see that the 
electronic origin of the LiH bands should lie 
0.18 cm~! towards the violet of the LiD bands 
an effect unfortunately too small to detect in 
this case. 

These facts together with the data accom- 
panying Eqs. (28) and (29) then enable us to 
draw three conclusions, viz., (1) the heavier 
isotopic molecule is the smaller, its potential 
curve (2) has the greater convexity at the 
minimum and (3) lies below that of the lighter 
isotope. These differences it will be recalled are in 
the same direction as those predicted by Kronig® 
from the interaction of nuclear and electronic 
motions (effect (2) above). 

By using Morse’s expression for the dis- 
sociation energy, D=w,?/4x, we have after cor- 
recting for the zero-point vibrational energy 
D,” =2.53 volts for LiD and 2.54 volts for LiH 
where the difference is not significant. If then the 
upper state dissociates into a normal H atom and 
a 2°P Li atom we have Dy’ =1.15 volts for LiD 
and D,’=1.18 volts for LiH. These results due 
to the small number of vibrational levels in the 
ground state can hardly be depended on to more 
than +0.2 volt. 
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Effect of Electron Pressure on Plasma Electron Oscillations 


ERNEST G. LINDER, RCA 


Manufacturing Co., Inc., 


Camden, N. J. 


(Received March 18, 1936) 


A general equation for electron motion in a plasma is 
developed which includes a term arising from electron gas 
pressure. The resulting expression is 


PE/OP + (4rne?/m)E= (RT / m)V7%E, 


where é is electron displacement, » electron density, and 
T electron gas temperature. From this it is found that 


N their investigation of electron oscillations in 

ionized gases, Tonks and Langmuir! did not 
take into consideration forces arising from elec- 
tron gas pressure. The writer has found that if 
this force is included the previously not under- 
stood variation of resonance frequency with 
electron temperature may be explained. 

In order to determine the effect of electron 
pressure, consider the equation of motion, 


mt= —eE+dp/ndx, (1) 


where the last term represents the force arising 
from pressure. The electron displacement is 
denoted by & The electric field E arises from 
electron displacements. In an undisturbed plasma 
the electron and positive ion densities are uni- 
form and equal, and the net charge in any 
element of volume is zero, hence the plasma is 
field free. However, if electron displacements 
occur, the net charge in a volume element is no 
longer zero, in general, and electric forces are 


set up. The expression for this force, is! 


E=4rneé. (2) 


Departures from a uniform distribution will 
give rise also to pressure gradients. It will be 
shown that the forces thus arising may be of the 
same order of magnitude as the force EZ, and 
that their inclusion in the theory gives a variation 
of frequency with electron gas temperature. 

If the undisturbed electron density is taken 
to be m, the change due to any displacement 
E(x) is 


bn = nAE/dx, (3) 


so that n=n(1+0/dx). (4) 


1. Tonks and I. Langmuir, Phys. Rev. 33, 195 (1929). 


fo. 


5 


the possible frequencies of free vibration form a series 
given by f;=(k7T/A2m+ne?/rm)'. The lower limit corre 
sponds to the Tonks-Langmuir value (ne?/wm)!, while the 
other frequencies depend upon the possible standing waves 
which may exist. The theory explains the observed varia 


tion of frequency with electron gas temperature 


Now p=nkT, 
Op On Ore 

hence = kT—=nykT—. (5) 
Ox Ox Ox? 


Substituting (5) and (2) in (1) yields 


" ° g mi . oe \ 6 
ar m m ox" 


A more general equation can be derived from 
Maxwell's equations, and a force equation 
similar to (1), wherein dp/dx is replaced by 
grad p. The procedure is the same as that given 
by Tonks and Langmuir,' except that the term 
involving pressure is carried through. The re- 
sulting equation is the same as (6) if d°£/dx° be 
replaced by VE. 

For a plane wave the possible frequencies are 
given by, 

fi=(RT/X2m+ne/xm)}, (7) 


where the \,’s are the possible wave-lengths 
determined by the standing wave systems which 
may exist. For a given m and 7, there is therefore 
a series of possible frequencies, the lower limit 
of which is the Tonks-Langmuir frequency 
(ne’/xm)'. This series, however, appears experi- 
mentally as a single broad band due to the 
breadth of the individual resonances and the 
closeness of spacing. 

It is that 
wave-lengths lies between a length comparable 


assumed the possible range of 
to the size of the containing vessel, and a length 
comparable to the Debye distance for the plasma. 
It will be shown that the experimental data can 
be fitted by choosing a wave-length within this 
region. 


? 
5 





74 ERNEST 


2.0 





























3 - 
T+10* 


Fic. 1. Showing combinations of m, the number of 
electrons per cc, and 7, the electron gas temperature, which 
give resonance under applied frequency of 1.6X 10°. 


In Fig. 1 are plotted data, some of which were 
taken from Tonks’ paper,? and some of which, 
previously unpublished, were kindly supplied to 
the writer by him. These represent combinations 
of m and T, which gave resonance when the 
plasma was placed between condenser plates to 
which a fixed frequency (1.6108) was applied. 
The linear relation is represented by (7) if both 
sides are squared, and f is regarded as a constant. 
The two lines, (a) and (b), correspond to two 
resonances which were observed, but for which 
there is no explanation. They apparently repre- 
sent two modes of oscillation, and may be related 
to plasma shape.’ It will be noted that there is a 
very definite effect due to electron temperature, 
the values of » required for resonance varying 
by as much as 2 : 1, over the temperature range 
investigated. 

The slopes of these lines should be, according 
to (7), 

dn/dT = —rk/e’?* (8) 


Tonks, Phys. Rev. 37, 1458 (1931). 
. Tonks, Phys. Rev. 38, 1219 (1931). 


oe 


(; 


LINDER 
To fit the data it is found that the effective 
wave-lengths must be \, =0.36 cm, and A, =0.22 
cm. The plasma diameters in the two cases 
were about 1.9 cm and 2.5 cm, respectively. 
The discharge did not fill the entire cross section 
of the tube, which was 3 cm in diameter. The 
Debye distances, Ap=(Rk7T/4rne*)', are 0.064 
and 0.038 cm, respectively. It is evident that 
the necessary \’s lie within the required range 
in both cases. 

The width of the 
estimated as follows: 
length limit to be given by a@Ap, where a is a 
small numerical constant, found by Langmuir* 
to be 3.31. Since \p = (kT /4rne*)? we may make 
use of (7) and write the limiting frequency as 


resonance band may be 


Assume the small wave- 


4r? tyne} ne» } 
moh 
ae 7m 7m 
Hence the maximum frequency is 2.14 times the 
lowest frequency, giving a band width of about 
one octave, which is in agreement with experi- 
ment. 

The contribution of the pressure term to the 
frequency may be computed by using the above- 
determined values of \. It is found that the 
frequencies predicted for the (a) and (b) reso- 
nances are 50 percent and 45 percent greater, 
respectively, for a given electron density, than 
when pressure is not considered. The available 
data do not permit a decisive check of this 
because of the unexplained fact that several 
resonances are observed. The agreement is better 
with the (a) type but worse with the (b). 

It should be noted that Eq. (8) for the slope 
is independent of any constant factor multi- 
plying f in (7), for example, to correct for 
plasma configuration as was done by Tonks.’ 
The expression for slope would be unchanged by 
any such factor. 


41. Langmuir, Proc. Nat. Acad. Sci. 14, 633 (1928 
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An Interpretation of Page’s ‘‘New Relativity” 


H. P. Ropertson, Princeton University 


(Received March 4, 1936) 


Page's proposed extension of the theory of relativity is 
examined from the standpoint of the general kinematical 
theory developed by the present author, of which it is 
found to be a special case. The contention that this 
extension renders Einstein’s theory untenable would seem 


to be without real foundation ; for although such an attempt 


N a recent issue of this Review L. Page! has 

developed a relativistic kinematics suitable 
for observers suffering constant relative accelera- 
tions, and has announced his intention of setting 
up a corresponding theory of the electromagnetic 
field. We here examine this program, in the light 
of previous related work, in order to determine 
whether or to what extent its results are incom- 
patible with the existing relativistic theories. 
As a result of this examination we are led to the 
conclusion that Page's kinematics is included as 
a special case in the general kinematical theory 
formulated some months ago by the present 
author? for a different purpose, and that as such, 
although outside the frame of the special theory 
of relativity, it may readily be understood from 
the standpoint of the general theory. It is further 
pointed out that an electrodynamical theory, 
based on transformations of the kind employed 
by Page, has been developed by H. Bateman® 
and by E. Cunningham.‘ For the sake of concise- 
ness and completeness the terminology of the 
theory of groups is at times invoked, but it is to 
be emphasized that all results are obtained by a 
straight-forward analysis, for an understanding 
of which a knowledge of the theory of groups is 
not required. 


1. KINEMATICAL PRELIMINARIES 


We begin, then, with a brief statement of the 
results of our kinematical analysis, referred to 


' Page, Phys. Rev. 49, 254 (1936). 

? Robertson, Astrophys. J. 82, 284 (1935) 

% Bateman, Proc. Lond. Math. Soc. 7, 70 (1908); 8, 223 
(1909). 

* Cunningham, Proc. Lond. Math. Soc. 8, 77 (1909). It is 
of interest to note that this paper bears the title ‘‘The 
Principle of Relativity in Electrodynamics and an Exten- 
sion Thereof,” 


~! 


wn 


to treat of observers in motions other than those contem 
plated in the special theory of relativity must lead to a 
theory transcending this latter, Page's theory is shown to 
be readily understandable in terms of a simple, though 
somewhat artificial, space time of a type met in the general 


theory of relativity. 


above,’ which will be of service in the following. 
It is there assumed that, for the purposes ir 
mind, the space-time region in question admits 
a three-parameter family of possible equivalent 
observers, in the sense that through each event 
there passes the world-line of one (and in general 
only one) such observer A, and that his view of 
the development of the whole, as obtained with 
the aid of a proper clock, light signals and a 
theodolite, is intrinsically indistinguishable from 
that of any other observer A’ in the family. 
It is then shown that under these conditions A 
can define operationally a set of coordinates 
n*=(rnv), and can introduce a Riemannian 


metric 
ds*=dr*— &(1r)du’, 
du? =dn*?+07(n)[d@+sin? ddy* J, 5) 
in which 
o(n)=sinh 7, » or sing (1.2) 


according as the Riemannian curvature & of 
the auxiliary three-space with metric dx’ is —1, 
0 or +1, respectively; the function £(7)20 
characterizes the type of relative motion between 
any two equivalent observers of the family, and 
is so far arbitrary. This metric ds* possesses the 
important properties: 


(a) The world-line of each of the fundamental particle 
observers A’ is a geodesic »*=const. (a=1, 2, 3) of 
ds*, along which the interval s measures the clock time 
rof A’. 

(b) The world-lines of all light paths are the null-geodesics 


of the metric. 


The coordinates »* of any event E in the 
region are determined operationally by A as 
follows. Let 4, be the time at which A must send 
out a light signal in order that it pass through 


wn 
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E, and let & be the time at which a signal sent 
out from E reaches the world-line of A. Then 
the coordinates 7, » assigned the event E by A 
are obtained from his clock readings th), tf by 


means of the equations 


2F (7) = F(t) + Fh), 2n = F(te) — F(t), (1.3) 


where F(r) = Sfdr/i(r) (1.4) 


is any primitive of 1/&(r). The angles @ and y 
are those made by the direction of either light 
signal with the zenith and meridian, respectively, 
of A’s theodolite. The transformation *—n’# 
from the coordinates employed by A to those 
employed by any other fundamental observer A’ 
is given by 7r’=7r and a transformation *—7’* 
of the spatial coordinates which leaves du*, and 
therefore ds*, invariant in form as well as in fact; 
that this latter transformation is a member of a 
six-parameter group Gg, of motions is accounted 
for by the fact that each of the «* particle- 
observers may orient his theodolite in any one of 
2% ways. For the case k=0, of greatest interest 
for our present purposes, this group Gg is the 
group of Euclidean motions, consisting of the 
three-parameter group T; of Euclidean transla- 
tions and the group R; of rotations ; we may then 
for most purposes ignore the latter, and confine 
ourselves to the group T; of translations between 
observers with similarly oriented theodolites. 


2. CONSTANCY OF LIGHT VELOCITY 


The kinematical theory given in résumé above 
is clearly the general solution of the equivalence 
problem stated by Page in Part 1 of his paper. 
However he, following E. A. Milne, would have 
the theory expressed in terms of a time ¢ and a 
distance r assigned to the event E as in the 
special theory of relativity, i.e. he would replace 
our 7, 9 or 4), & by their expressions in terms of 
the variables 


t=3(4+h), r=}c(te—h). (2.1) 


This in itself is of course no restriction on the 
generality of the solution, although it leads to 
the more cumbersome expression 


ds? =[#(7r)/E(t) E(t) ] {df 
—[dr?+ R2(d@+sin? Ody) ]/c?} (2.2) 


for the line element (1.1), where 
R°(t, r) =c?&(t)) E(te)07(n) (2.3) 


and +, 7, 4, fg are to be expressed in terms of ¢, 7 
with the aid of (1.3), (2.1). 

Now the path of any light signal which 
intersects the world-line of A is given by 


r=+c(t—to), #=const., y=const., 


and may therefore be interpreted as a radial 
straight in the Euclidean (xyz)-map constructed 
from the Cartesian coordinates (xyz) obtained 
from the polar coordinates (réy) in the usual 
way ; hence any light signal sent by A or received 
directly by A is propagated rectilinearly in the 
(xyz)-map with the constant coordinate velocity 
dr /dt= +c. However, it is readily seen from (2.2) 
that other light paths do not in general possess 
this property, and that to follow Milne and Page 
in requiring that they do, places a severe re- 
striction on the function £(r) characterizing the 
relative motions of the fundamental particle- 
observers. We remark in passing that this result 
in no way contradicts that obtained in section 3 
of the paper referred to above,? where it is 
shown that, at least in the open models k= —1, 0, 
coordinates (7X YZ) may always be introduced 
in which light is propagated rectilinearly with 
constant coordinate velocity c; these latter are 
defined in quite another way than the (trys) 
here employed, and are such that the transfor- 
mations from one particle-observer to another 
constitute a six-parameter subgroup of the 
general Lorentz group. 

In order to interpret Page’s work in terms of 
the kinematics outlined above, we must first 
investigate the restriction imposed by his re- 
quirement that all light signals, and not only 
those which A is in position to observe directly, 
are propagated rectilinearly with constant co- 
ordinate velocity ¢c in the Euclidean map con- 
structed from the coordinates (xyz). It is seen 
immediately from (2.2) that, in order for all null 
elements ds=0 to be null elements in the 
Minkowskian (txyz)-map, the function R(t, r) 
defined by (2.3) must be identically equal to r; 
that this condition is also sufficient follows from 
the fact that the null geodesics of a conformally- 
flat manifold are at the same time null geodesics 
of the associated flat. On replacing &(r) by its 
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value 1/F’(r) in terms of the derivative F’ of F, 
the condition R=r yields the differential equa- 
tions 

F’(a) F’(r)(r—a)?= { 20 (F(r) — F(a))/2]}*? (2.4) 


for F(r), where we have expressed all variables 
in terms of 4, tf and have then set 4;=a, ¢ =r. 
The solution F(r) is given, to within an additive 
constant, by a quadrature, and the corresponding 
forms of &(r) by differentiation. The results of 
this computation may be expressed in terms of 
the general linear fractional function 


f(r) =(Ar+B)/(Cr+D), A=AD-—BC>O, (2.5) 
in the various cases & as follows: 
F=f, §=(Cr+D)?/A, (k=0); (2.6) 


F=log f, §=(Ar+B)(Cr+D)/A, 
(k=—1); (2.6’) 


Ste 


=[(Ar+B)?+(Cr+D)*]/2A, (k= +1). (2.6’’) 
Now in a certain sense 
A=ck(r)o(u) (2.7) 


may be considered as a significant measure of 
the ‘‘distance”’ at ‘‘time’’ +r between two particle- 
observers with the fixed u-interval u, and we 
define 


v=dyX/dr, a=d*\/dr’, (2.8) 


as the corresponding measures of the relative 
“velocity’’ and “acceleration,” respectively. It 
now follows from the results of the preceding 
paragraph that the requirement of constant 
coordinate-velocity of light in the (xyz)-map 
implies that the relative motion of any two funda- 
mental particle-observers must be such that their 
relative a-acceleration is constant. 

An examination of the three sets of equations 
(2.6) shows that the only case in which the 
relative velocity v can be zero for all 7 is that 
in which k=0, C=0; we then have 


D/A: f= f, 


wre 
II 


r=Den/A, (2.9) 


and the original line element (1.1) is_ itself 
Minkowskian on change of spatial scale. The 
only other case in which the acceleration a 
vanishes throughout the motion is that in which 
k=—1, C=0, whence 


7 


t=r; t=r7 cosh 7, =crsinh y, (2.10) 


uw 
~J 
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and the new line element (2.2) is the Minkow- 
skian form on which the special theory of 
relativity is based; the original metric (1.1) is in 
this case the “‘polar’’ form implied in the work 
of Milne. 

Finally, to come to the case adopted by Page, 
we note that he requires (pp. 260, 263) that the 
relative motion of two particle-observers be such 
that at some instant they coincide spatially and 
are then momentarily at relative rest. The first 
of these conditions allows us to reject the case 
(2.6), for in it the distance function &(r) has 
only complex roots; the second condition then 
throws out the hyperbolic case (2.6’), in which 
the two roots of £(r) are real but distinct. Thus 
we find, as the only case in which two particle- 
observers are at some time momentarily in coinci- 
dence and at relative rest, the parabolic case (2.6), 
in which &(r) has a double root and the spaces 
t=const. are Euclidean. The remainder of our 
investigation will be devoted to a discussion of 
this case. 


3. TRANSFORMATION THEORY 


On taking as the origin of + that event at 
which all «* particle-observers momentarily 
coincide, the invariant metric (1.1) defined by 
(2.6) may be written 


ds*=dr? —(gr*/2c)*{ dn? +dn2+dns], (3.1) 


L 


where the », (a=1, 2, 3) are the Cartesian 
coordinates associated with the polar coordinates 
(ny) in the usual way, and ¢ is the (constant) ac- 
celeration between two particle-observers whose 
(constant) u-distance is unity. Our first task is 
that of discovering the transformation (rn,) 

»(tx,) which throws (3.1) into the conformally- 
Minkowskian form (2.2), where the x, are the 
Cartesian coordinates associated with (r@y). We 
may now take 


F(r) = —2c/er, (3.2) 
and the Eqs. (1.3), (2.1) yield as the required 
transformation 


7=(f—?r*/c*) /t, Na = 2x,/(P—Pr/c*)e; (3.3) 


the form (2.2) of the metric is then found to be 


ds*=(t t)*[ de — (dx? +dx.*+dx;") c?). (3.4) 
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The transformation (3.3) has merely the effect 
of allowing the given observer A at =O to 
introduce a new coordinate mesh, of the kind 
demanded by Page, in which A is at the origin 
r=0(0 of the new spatial coordinates x,, and in 
which the new temporal variable ¢ is also a 
direct measure of proper-time along his world- 
line. In order to find the relations between the 
coordinates (/x,) employed by A and the co- 
ordinates (¢’x,’) similarly introduced by another 
particle-observer A’ at no =8., we must consider 
the group of motions admitted by the metric 
(3.1); as indicated above, we may here confine 
ourselves to the subgroup T; of translations 
between observers with similarly oriented theo- 
dolites. Hence the transformation A—-A’ is 
given by 

, 


7 


Na’ = Na— Bay (3.5) 


and the transformation which it induces on the 
corresponding Minkowskian coordinates is de- 
fined implicitly by 

(t? —9"? /c?)/t =(P—Pr/c*)/t (=7) 
(3.6) 


, 


Xa /t =x_/t— ¢erTB./2. 


Now this transformation must leave (3.4) in- 
variant in form as well as in fact, and it must 
therefore be a member of the general conformal 
group in four dimensions; hence it must be 
compounded of at most dilatations, translations, 
rotations (including reflections), and inversions. 
The explicit form, as found directly from (3.6), is 


t'=0't/R?, 
Xa —ba=0?[ xa — (1+2E5b.x,/b*) bq ]/R’, 


(3.7) 


where R?(txa3 ba) == (Xa td.)? -—CF, 


b, = 20°Ba/B’¢, ==p.", 


this transformation may in fact be broken down 


(3.8) 


San Bh 2 
5? = 5),’; 


into the reflection 
Ea=Xa—2(1+2).x,./b*)b, (3.9) 
in the plane through x,= —0, orthogonal to the 
vector },, followed by the space-time inversion 
t’ = 5*t/R?, Xe —b,=b2(F.+6,)/R? (3.10) 
in the pseudo-sphere 


R*(t%.; ba) = R2(tx,; b,) = 8? (3.11) 


of radius b and center at the event (0, —d,). 
Thus the group I; induced by T; is that subgroup 
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of the general conformal group which consists of 
the «* transformations t(d,) of the form (3.7), 
and which possesses the differential invariant 
(3.4); the law of composition of this (Abelian) 
group I; is given by 

t(da)t (Ca) = t(bat+Ca). (3.12) 


That (3.7) is indeed the explicit form of the 
transformation derived by Page can be seen by 
comparing it with his Eqs. (46), (47), (68), on 
going over to the variables there employed. The 
particle-observer A’ at x,/=0 moves radially 
away from A in the direction of the vector } 
with a coordinate-velocity v which satisfies the 
equation 

v(1—v?/c?)-!=Bet, (3.13) 
as can be found directly from (3.7); hence the / 
rate of change of the left-hand side of (3.13), 
which Page calls the “relativity acceleration” of 
A’ relative to A, is the constant Byg—and is the 
same as their relative a-acceleration, as defined 
by (2.8). The geometrical theorems derived by 
Page in Part 4 of his paper are consequences of 
general theorems on inversions, as applied to 
the spatial inversion defined by (3.10) at time 
t=f’ =(0. 

We also see, in retrospect, that the transfor- 
mation (rn,)—+(tx.), defined by (3.3) for a given 
observer A, is itself interpretable as an inversion 
of the variables 


by? /c?r, bona (bp =2c?/ ¢) 


in the pseudo-sphere of imaginary ‘‘radius"’ id, 
about the origin. 

Our approach has been based on the existence 
of an invariant interval ds* which satisfies the 
conditions (a), (b) stated in section 1 above, 
and which therefore offers an invariantive de- 
scription of the only elements so far involved in 
the theory—namely, the equivalent particle- 
observers, their proper clocks and theodolites, 
and all light signals. In the light of this, and 
considering the rather arbitrary definition of / 
and +, it is difficult to see why Page's “physical 
interval” cdf —dx* should be expected to appear 
as an invariant. It is to be noted, in this con- 
nection, that the content of Page’s Eq. (32) is 
equivalent to the fact that our ds? is invariant 
in the two-dimensional case 6=const., ¥=const. ; 


“his expressions (31), (32) are not, however, 
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invariants of the transformation (tr)—(?’r’) in 
the technical sense, for they depend on the 
parameter ¢ of the transformation. On the 
other hand, the differential expression (69) which 
he writes down for the full four-dimensional case 
is not at all equivalent to our differential 
invariant (3.4); even if it be modified by the 
same treatment as that accorded (32) in bringing 
it into invariant form, it could at best be ex- 
pected to be invariant under a one-parameter 
subgroup of I;, and not under the full group. 


4. ASSOCIATED PARTICLE-OBSERVERS 


It remains to clarify the role of the «* possible 
particle-observers which Page associates with 
each given particle-observer A of the type dealt 
with above. These auxiliary observers have as 
their world-lines the curves x,=const.; in order 
to avoid confusion with the primary particle- 
observer A we shall refer to the former as test- 
particles. Computation shows that the world- 
lines x,.=const. are not geodesics of the metric 
ds*, unless x,=0—hence, in Page’s terminology, 
the reference systems employed by two observers 
A and A’, while equivalent, are not homogeneous. 
Now we have examined,° in connection with the 
general kinematical theory described in section 1 
above, the possible modes of motion of such a 
six-parameter family of test particles in a space 
time possessing the uniformity properties here 
implied, and have found that they must satisfy 
differential equations of the form 


d?n* m | dn’ dn? dr drdn* 
+| . -1(-. ) 694 — , (4.1) 
ds* ve) ds ds ds ds ds 


where n°=7r, the Christoffel symbols are com- 
puted from the coefficients g,, of (1.1), and I is 
some function of at most 7 and dr/ds. On 
evaluating these expressions for the metric (3.1) 
and taking »* as defined by (3.3), in which the 
X%q are constants, these equations of motion 


reduce to 
dt/ds=t/r, ('+1/r)xadr/ds=0; (4.2) 


hence the «* world-lines x,=const. associated 
with each of the particle-observers A are solutions 
of (4.1) provided we set 


I'(r, dr/ds) = —1/r. (4.3) 
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These results allow us to conclude, in accord- 
ance with the general theory developed in 
section 5 of the paper cited above,’ that the «® 
solutions of (4.1) fall into classes characterized 
by a parameter d>0. The class d=0 consists of 
the «* world-lines of the fundamental particle- 
observers, the class d= ~ of the «® light paths 
in space time, and each class d(#0, ~) of the 
«° test particles which are at the fixed x distance 
d from the particle-observer with which each of 
them is associated. Further, the motion of each 
member of such a class d is such that its A 
distance, as defined by (2.7), from the observer 
A with which it is associated, approaches asymp- 
totically the value d. 

Finally, these considerations unearth a sur- 
prising formal connection between Page’s work 
and the transformation theory of the de Sitter 
universe in the general theory of relativity. To 
exhibit this connection, we first note that in any 
case in which [ is a function of 7+ alone, the 
equatrons (4.1) may be interpreted as the 
equations of the «* geodesics of the metric 


dS=eS lads? (4.4) 


conformal to the original metric (1.1). For the 
case (3.4) of interest here, this new metric may, 
on subjecting ¢ to the transformation 


t=(bo c)ec?!>o, (bp =2c? ¢) (4.5) 


be taken in the form 

dS? = dz? — e~***!%0(dx;?+dx2?+dx;3"). (4.6) 
This invariant metric satisfies the conditions (a) 
and (b) of section 1 above, except that it is no 
longer a direct measure of clock time along the 


geodesic paths of the fundamental particle- 
observers. It is, on the other hand, the stationary 


form of a contracting de Sitter universe of constant 


Riemannian curvature —c*/b,?; the relation 
between a particle-observer A and an associated 
test particle x,=const. in the space time (3.4) 
of interest here is, to within the time transfor- 
mation (4.5) above, formally the same as that 
between two associated free observers in the 
de Sitter universe (4.6), and the relation between 





* Robertson, Astrophys. J. 83, 187 (1936); in particular, 
sections 4 and 5. 

*Cf. the discussion of the de Sitter universe by the 
present author in Phil. Mag. 5, 835 (1928). Eq. (11) of this 
discussion are in fact equivalent to those defining Page's 
transformation (68) between equivalent particle-observers. 
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two equivalent observers A, A’ here is formally 
the same as that between two free observers in 
general relative motion in the de Sitter universe. ® 


5. SUMMARY AND CONCLUSION 


This examination of Page’s work, from the 
standpoint of our general kinematics, has shown 
that it deals with one of the three possible cases 
in which there exists a set of «* equivalent 
particle-observers with respect to whom light is 
propagated rectilinearly with a constant co- 
ordinate velocity ¢c in a coordinate system (tx,) 
defined as in the special theory of relativity. 
In each of these cases there exists a significant 
Riemannian metric, of the type employed in the 
cosmological applications of the general theory 
of relativity, which is invariant under the group 
Ge. of motions which describe the equivalence ; 
it is apparent a priori that this group is a six- 
parameter subgroup of the general conformally- 
Minkowskian group of transformations ou the 
variables (¢x,), and that it is the direct product 
of the group R; of Euclidean rotations on the 
spatial variables x, and a three-parameter group 
I; which has as its minimal invariant varieties 
the spaces r=const. We have in particular 
examined this group I; in the case adopted by 
Page, and have shown how it may be generated 
from translations, reflections and simple in- 
versions. 

The three cases thus found are characterized 
by the fact that in them the fundamental 
particle-observers undergo a constant relative 
a-acceleration, which in Page’s case agrees with 
the ¢ rate of change of the “‘relativistic’’ velocity 
defined by the left-hand side of (3.13); we must 
regard this agreement as a happy coincidence, 
however, for it is difficult to see what significance 
is to be attributed to Page’s procedure outside 
the special relativistic theory of electrodynamics, 
where it leads to the usual classical expression 
for the ponderomotive force (and not to the 
relativistic four-vector by which it should even 
there be replaced). It would seem that, insofar 
as the present line of attack is concerned, this 
leaves little justification for Page’s “hope of 
finding a rational detailed description of atomic 
structure” by an investigation of “equivalent 
reference systems having other types of motion 
{than with constant acceleration ], particularly 





relative rotation” ; indeed, it has long since been 
pointed out by E. Cunningham’ that “no such 
relation [as that defined by the conformal 
group | can express anything corresponding to a 
rotational motion of the space-frame of refer- 
ence.” It would further seem, particularly in 
view of the severe restrictions which their intro- 
duction implies, that the quasi-Minkowskian 
coordinates (éxyz) are of no very great intrinsic 
significance, although possibly of convenience 
for some purposes. But, above all, it is to be 
remembered that we are here on the kinematical 
level of description, and that the question of 
whether or not two observers are dynamically 
equivalent can be decided only in terms of the 
more complete physical theory of which the 
kinematics is but a noncommittal tool; thus, 
while we are willing to entertain the view that 
a general kinematical theory such as that out- 
lined in section 1 above may play a useful role 
in comological speculations, we should not be 
disposed to attempt to reconcile it with the 
obvious dynamical inequivalence of observers in 
general relative motions in terrestial or atomic 
problems. Otherwise stated, either these more 
general types do represent possible motions of 
dynamically equivalent observers, in which case 
we should expect to have recourse to methods 
and concepts of the type met in the general 
theory of relativity, or their dynamical inequiva- 
lence should be expected to set serious limitations 
to the usefulness of the kinematics. 

Insofar as electrodynamics is concerned, we 
confine ourselves to the remark that the problem 
of transforming the electromagnetic field under 
the general conformal group has been given a 
detailed solution by Bateman and by Cunning- 
ham in the papers cited above*:‘*; we do not 
here go into the question of whether or to what 
extent their theory is equivalent to the phe- 
nomenological theory of electrodynamics based 
on the line element (3.4) as in the general theory 
of relativity. 

Finally, I should Jike to thank Princeton 
University for granting me sabbatical leave 
during the present academic year, and to express 
my appreciation to the California Institute of 
Technology for the facilities for work which it 
has placed at my disposal during this time. 


7 Cunningham, The Principle of Relativity (Cambridge, 
1914), p. 89. Cf. also p. 79 of reference 4 above. 
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The design, technique of operation, and performance of 
the Round Hill electrostatic generator are presented in 
some detail. The problems of operating wide paper belts, 
of eliminating vibration, and of controlling humidity are 
discussed. The original belt charging system, with the belts 
operated at saturation charge density, is described and 
its range of usefulness is indicated. A consideration of the 
problems of voltage control and voltage steadiness at 
reduced charging currents leads to two modified designs of 
the belt charging system. The measurement of voltage by 
means of the generating voltmeter is discussed with special 
emphasis on the precautions necessary to eliminate or to 
correct for sources of systematic error. Under the condi- 


N 1933, there appeared a general discussion! 

of the principles involved in the electrostatic 
production of high voltage by means of belt 
generators operating at atmospheric pressure. 
There was included a brief description of the 
external structure of the Round Hill electrostatic 
generator. During the intervening period a belt- 
charging system has been designed, constructed, 
and improved, so that a very steady, controllable 
voltage is now available for nuclear research. 
Complete generator performance data have been 
obtained through the utilization of the generating 
voltmeter which, with proper precautions, has 
proved to be a convenient and reliable instrument 
for voltage measurement. 


EXTERNAL STRUCTURE 


The generator structure consists essentially of 
two spherical terminals supported by cylindrical 
insulating columns on steel trucks. The spheres 
are made of aluminum alloy 15 feet in diameter 
with walls } inch thick. The insulating columns 
are cylinders of Textolite 22 feet in length of free 
insulation, 6 feet in diameter, and of $ inch wall 
thickness. The steel trucks act as heavy bases for 
the generator units and house the driving and 
An electrically-driven 


accessory equipment. 


* National Research Fellow, July, 1933 to March, 1935. 
1R. J. Van de Graaff, K. T. Compton and L. C. Van 
Atta, Phys. Rev. 43, 149 (1933). 


tions realizable with the Round Hill generator the probable 
error in voltage measurements is less than one percent on 
a relative scale and about two percent on an absolute 
scale. The magnitude of voltage fluctuations is less than 
0.1 percent. The generator performance data for the 
positive and negative terminals are given in graphical form. 
The maximum charging current is 2.1 ma and is practically 
independent of voltage. The highest voltage obtained 
consistently without sparking is 2.4 megavolts positive 
and 2.7 megavolts negative, giving 5.1 megavolts between 
the terminals. At this voltage there is 1.1 ma of current 


available for application to an accelerating tube. 





Fic. 1. Electrostatic generator and housing. 


winch, with an endless cable passing under the 
trucks, is used to transport the generator units 
along a 14-foot gauge railroad track. As shown in 
Fig. 1, a steel building 140 feet long, 75 feet 
wide, and 75 feet high houses the generator and 
the laboratory and shop rooms. Structural girders 
supporting the roof limit the unobstructed ceiling 
height to 60 feet. Further details relating to the 
external structure have been given previously.! 


BELTS AND PULLEYS 


In each generator unit electrical charge is 
conveyed to the spherical terminal by means of 
two simple belts running inside the supporting 
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Fic. 2. Sectionalized view of the interior of one generator unit showing the belt-charging system and auxiliary equip 
ment in terminal and truck. 




















ROUND HILL ELE 
column as shown in Fig. 2. The belt material is 
electrical insulation paper® 0.017 inch thick and 
47 inches wide. This material satisfies very well 
the electrical and mechanical requirements and 
has the advantage of being inexpensive as 
compared with silk or rubberized fabric. The 
seams are made on a 45-degree bias for greatet 
strength and smoother operation. The edges to 
be joined are tapered to bring all laminations of 
the paper into the joint and to preserve a uniform 
belt thickness. A celluloid-base belt cement* is 
used for the joints instead of glue because of its 
superior electrical properties and its flexibility at 
low humidity. 

Proper initial conditioning is essential to the 
successful operation of such a wide paper belt. 
After being mounted between pulleys in the 
generator column the belt first undergoes a two- 
day period of stretching under a tension of about 
1000 pounds between pulleys at a relative 
humidity greater than 75 percent. The belt is 
then dried over a period of about three days 
while the humidity is gradually reduced to 
approximately 25 percent. The belt lengthens by 
0.3 of its length during the 
stretching period and contracts by about 0.2 
percent upon drying. The pulley is not allowed to 
turn until the belt is dry since a damp belt tends 
to crease as it runs over the pulley. 


about percent 


The belts are carried on steel pulleys 6 inches in 
diameter and 50 inches long, the lower pair 
acting as drivers and the upper pair as idlers. 
The pulleys are of seamless steel tubing having 
chrome-vanadium forged end plugs which extend 
with reduced diameter to provide substantial 
shafts. Self-aligning ball bearings are shrunk 
onto the shafts and are housed in solid pillow- 
blocks. One shaft of each of the lower pulleys 
extends through the bearing for a belt drive to a 
10 hp induction motor. An endless rubberized 
fabric flat belt and a constant-tension motor 
base* are employed to obtain a quiet and 
efficient drive. Each pair of pulleys is mounted on 
a heavy square frame of channel iron so that the 
pulleys and frame can be isolated as a unit, 
mechanically and electrically. 


? John A. Manning Paper Co., Troy, N. Y. 

3 Sea Lion Waterproof Belt Cement, Chicago Belting Co., 
113 North Green Street, Chicago, Ill. 

‘The Rockwood Mfg. Co., Indianapolis, Ind. 


TROSTATI( 





763 


GENERATOR 

Because of the small stretch in dry paper belts, 
the pulleys require little crowning. The crowning 
the 6 
diameter, applied as a uniform taper on the last 


actually employed is 3); inch in inches 
4 inches of pulley at each end. Since the pulley is 
3 inches wider than the belt only part of this 
crowning ts effective when the belt is centered 
The pillow-blocks of each lower pulley are 
mounted on heavy threaded studs in the sup 
porting frame to provide a sensitive means for 
adjusting the level of the pulley, since the 
centering of the belt, as well as its tension, 
depends critically upon this adjustment. 

Two pulley speeds are provided, 3600 r.p.m. 
and 2880 r.p.m., corresponding to belt speeds of 
5650 f.p.m. and 4520 f.p.m. At the higher speed 
the two belts enter the sphere at the rate of 
44,200 sq. ft. per minute. The windage loss per 
belt is approximately 4 kw at the higher speed 
and 2 kw at the lower speed. The no-load power 
consumed by both units of the generator in 
operation at the higher belt speed is then 16 kw, 


while the full-load power consumed is 26 kw. 


VIBRATION ISOLATION 


Because of the necessity of using sensitive 
the 
chanical vibration and noise must be reduced to a 


measuring instruments in terminal, me- 
minimum. This has been accomplished by me- 
chanically isolating all moving parts, as indicated 
in Fig. 2, and by absorbing inside the terminal 
such noise as is transmitted to it. In addition the 
heavy pulleys for the paper belts were dynami- 
cally balanced until the their 
vibration when rotating at 3600 r.p.m. in free 


amplitude of 


support was less than 0.002 inch. 

For isolating the upper and lower pulley frames 
and all the heavy machinery a soft rubber® was 
used. In the case of lighter equipment and 
machinery a sponge rubber® proved more satis- 
factory. In addition, all the equipment within the 
truck is mounted on a heavy false floor, which is 
supported at the corners on soft rubber. In all 
cases a one-inch thickness of rubber was used 
with the proper area to obtain a ten percent 
deformation in the rubber under the isolated 
load. The maximum load per square inch is about 
100 Ibs. for the soft rubber and 10 Ibs. for the 


5 United States Rubber Co,, Passaic, N. J. 
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sponge rubber. Because of the marked improve- 
ment in quietness of operation gained by such 
efforts, vibration isolation has been extended to 
include every piece of machinery associated with 
the generator. 

The interiors of the metal terminals would have 
extremely poor acoustical properties if a large 
area of sound absorbing material were not 
provided. This is partly accomplished by a 
is-inch sponge rubber flooring. Additional ab- 
sorbing area may be obtained by ‘‘quilting”’ 


t-inch flameproof felt® to the metal wall. 


HuMIDITY CONTROL 


In the case of any high-voltage, low-current 
generator, insulation requirements are severe. 
The unusually high average relative humidity 
and the fact that the hangar is not completely 
weather-tight, make this problem exceptionally 
serious at Round Hill. The average relative 
humidity is approximately 84 percent during the 
summer months and 74 percent during the winter 
months. Normally the hangar is about 300 yards 
from the sea, but during extreme high tides has 
been entirely surrounded by water. 

The space inside the column and truck is 
carefully enclosed to make humidity control 
practicable. Both heating and drying are em- 
ployed to reduce the humidity. The heating 
system consists of an automatically operated 
kerosene water heater mounted on the outside of 
the truck, a water circulation system which 
includes a gear pump, an automobile radiator 
inside the truck, and a blower for circulating the 
air of the truck through the radiator and into the 
column. A heater output equivalent to 4 kw is 
sufficient to maintain the truck and column 10°C 
above the surroundings, but ten kilowatts are 
available for emergency use. 

The drying system consists of three trays each 
containing 30 pounds of Silica Gel.? Two trays 
are mounted over small blowers in the truck. The 
third tray is reactivated outside the truck over a 
combination blower and heater which circulates 
air at 160°C through the tray. The trays in the 
truck are replaced alternately at 12-hour inter- 


® Merchandise No. 3565, American Felt Co., 211 Con- 
gress St., Boston, Mass. 
7 Silica Gel Corporation, Baltimore, Md. 
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vals. Water is removed from the air at the 
average rate of one-quarter pound per hour. 

The combined heating and drying systems 
described above and indicated in Fig. 2, maintain 
the inside humidity below 35 percent with an 
outside humidity of 100 percent. Experience has 
shown that this renders the maximum working 
voltage of the generator independent of weather 
conditions. 

The resistance of the Textolite column sup- 
porting each of the terminals has been maintained 
above 10'° ohms by coating the outer surface 
with ceresin wax and by keeping the interior of 
the column warm and dry. The insulation require- 
ments for the charging belts are especially 
stringent because of the high gradients imposed 
upon them in the region of the charging mecha- 
nism and at the entrance to the high voltage 
terminal. For this reason the belts are of electrical 
insulation paper and are maintained thoroughly 


dry. 


POWER SUPPLY IN THE TERMINALS 


Both a.c. and d.c. power are supplied to the 
laboratories in the high voltage terminals. Fo 
reasons of economy, 6 kw of d.c. power are 
initially developed to care for the major heating 
and lighting demands, while a 1 kw inverted 
rotary converter supplies a.c. where this is 
essential. The long vertical drive from a motor in 
the truck to the d.c. generator in the terminal was 
accomplished, as shown in Fig. 2, by the use of a 
large-diameter, flanged paper pulley at the top, a 
rubberized fabric endless belt, and a motor base‘ 





Fic. 3. View of the interior of one of the generator trucks 
showing the lower pulley frame. 
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Fic. 4. Schematic diagram of the original belt-charging 
system. 


designed to maintain constant tension on the 
belt. 

From a centralized power panel in each termi- 
nal, the a.c. voltage and frequency and the d.c. 
voltage can be controlled and measured. Con- 
venience outlets make these voltages available at 
the work benches. 


ORIGINAL BELT-CHARGING SYSTEM 


The original belt-charging system employed 
for this generator, shown schematically in Fig. 4, 
is equivalent in principle to that indicated in 
Figs. 2 and 3 of a former publication.' The 
structural frames carrying the upper and lower 
pulley systems are insulated from the terminal 
and truck, respectively. Voltage is applied to the 
lower pulley system from a filtered 20-kv rectified 
a.c. voltage source. Whereas there are two 
charging beltsin each generator unit the following 
discussion can be simplified by considering the 
system as comprised of only a single belt. The 
charging current is supplied at each end of the 
belt by fine wires® stretched across the face of the 


*0.003-inch bright, hard drawn, high carbon steel wire, 
John A. Roebling’s Sons, 107 Liberty St., New York City. 


belt at a distance of about 3} inch from the upper 
and lower pulleys as shown in Fig. 5 and con- 
nected, respectively, to the terminal and to 
ground. 

The charge sprayed on the belt by the lower 
spray wire under the excitation of the frame 
voltage is removed from the ascending belt 
inside the terminal by a collector connected to 
the insulated upper frame. The collector consists 
of a 0.003-inch wire inside a 3-inch diameter 
cylindrical shield connected electrically to the 
wire and slotted for the passage of the belt. The 
shield around the collector wire is almost es- 
sential for reliable operation since the collector 
current is then independent of the upper frame 
voltage and dependent only on the charge 
density on the belt. 

As the upper frame voltage increases due to the 
collector current, corona is produced at the uppet 
spray wire and charge is deposited on the 
descending belt. The equilibrium voltage of the 
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Fic. 5. Isometric drawing of belt showing the location of 
spray wires and collector shields. 


Fine wires as corona sources for belt charging were first 
used on the Carnegie electrostatic generator at the Depart 
ment of Terrestrial Magnetism, Washington, D. C. 
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upper frame is reached when the leakage through 
the belt of charge deposited by the upper spray 
wire equals the collector current. In effect, both 
upper and lower spray wires must neutralize the 
charge brought to them on the belt before 
producing a net charge of the opposite sign. 
Unless the lower wire is spraying more than 
enough charge to saturate the belt, each spray 
wire current is approximately equal to the 
generator charging current delivered by the belt. 

The collectors on the lower frame enable the 
generator to operate by self-excitation due to the 
priming action of frictional charge produced on 
the paper belts in running over the pulleys. 
Because of the air conditioning inside the column, 
the generator consistently delivers its full current 
under self-excitation. In fact, when the lower 
collectors are removed the generator still de- 
velops half its full current with negative polarity. 
However, the self-excitation feature is not par- 
ticularly useful in the present case because of the 
desirability of definitely fixing the polarity and of 
controlling the charging current. 

Under normal conditions the system just 
described is operated with saturation charge 
density on both ascending and descending belts. 
This condition is realized by the proper choice of 
the lower frame voltage and spray wire spacing, 
and the upper spray wire spacing. The adjust- 
ment is not critical since excess charge deposited 
on the belt discharges to the pulley at the point 
of departure of the belt from the pulley. This 
charging system has been used considerably and 
has proved quite satisfactory so long as the belts 
are operated at saturation charge density. 


CONTROL OF GENERATOR VOLTAGE 


In order to realize the maximum usefulness of 
the generator, it is necessary to have a con- 
venient method for precise voltage control over a 
considerable range of voltage. In many cases it 
would be feasible to operate the belts at satura- 
tion charge density and to obtain voltage control 
by varying the belt speed or by imposing a 
controllable corona current drain on the terminal. 
In the present case it would not be practicable to 
control the speed of the four 10-hp induction 
motors. A corona device for limiting the voltage 
has been tested and found to be fairly satis- 
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factory electrically for controlling the voltage of 
the positive terminal. It is not satisfactory when 
used with the negative terminal because of 
persistent sparking to the corona device and is 
not convenient to use because of the dimensions 
of the apparatus and the precision required. 

The generator voltage may be controlled most 
conveniently by varying the charge density on 
the belts. For a given drain on the terminal the 
equilibrium voltage assumed by the generator 
will depend on the charging current. If the 
current drawn from the terminal is altered the 
charging current must be correspondingly altered 
to maintain a constant voltage. This adjustment 
may be made manually, or may be accomplished 
automatically by employing the output from a 
generating voltmeter to control the exciting 
voltage applied to the lower frame. 

To test this method of voltage control the 
exciting voltage applied’ to the lower frame was 
made continuously variable. As the charge 
density on the belts was reduced below satura- 
tion, the charging current and the equilibrium 
terminal voltage decreased in the expected man- 
ner. At reduced charge densities, however, there 
were present violent fluctuations in the charging 
current having the period of the belt, 1.6 c.p.s., 
and an amplitude too large to be tolerated. The 
typical curve in Fig. 6 shows the variation of 
charging current with lower frame voltage and 
the limits of the fluctuations encountered. 
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Fic. 6. Limits of the 1.6-c.p.s. current fluctuation en- 
countered with the original belt-charging system at reduced 
charge densities. 
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ELIMINATION OF CURRENT FLUCTUATIONS 


Apparently these fluctuations were initiated by 
the electrical irregularity in the belt at the glue 
joint and were amplified in effect by the inherent 
instability of the upper frame voltage. As a first 
step in the elimination of the current fluctua- 
tions, the electrical irregularity in the belt was 
greatly reduced by making the joint with a 
celluloid-base cement instead of glue. Then the 
instability of the upper frame voltage was cor- 
rected by a series of alterations in the design of 
the charging assembly. 

The instability of the upper frame voltage was 
due partly to the shape of the current-voltage 
characteristic curves of the spray wires. Current- 
voltage characteristic curves for spray wires at 
two different spacings from a grounded metal 
cylinder are shown in Fig. 7a. Two corresponding 
curves obtained with a spray comb of fine 
needles are given for comparison. The spray 
comb consists of a steel ribbon 0.010 inch 
x} inch with No. 16 beading needles soldered 3 
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Fic. 7. Current-voltage characteristic curves; a, for 
negative spray wires and combs with different corona gaps; 
b, for positive and negative combs with and without series 
resistance. 


inch apart and extending } 


the ribbon. Whereas the ideal characteristic 
curve of a corona device for belt-charging would 


inch from the edge of 


begin at the origin and rise linearly with voltage, 
the actual curves are nonlinear and may be 
expressed quite accurately by empirical relations 
of the form 


From Fig. 7a it appears that combs and larger 
spacings are to be preferred because of the lower 
corona threshold voltage of the combs and 
because of the flatter over-all characteristics 
obtained with larger spacings. The curves of 
Fig. 7b show the marked improvement obtained 
by adding resistance in series with the positive or 
negative comb. If it were not for other con- 
siderations a still flatter characteristic could be 
obtained by further decreasing the spray gap and 
increasing the series resistance, i.e., by increasing 
the ratio of ohmic to non-ohmic resistance in 
the circuit. 

Another factor contributing ‘to the instability 
of the upper frame voltage was the method of 
limiting that voltage. In this connection the 
charge density on the ascending and descending 
sides of the belts should be equal in order to 
realize maximum current from the belts and to 
prevent corona loss of unbound charge from the 
more heavily charged side. That is, the upper 
frame should be maintained at such a voltage 
that the current carried away by the descending 
belts is approximately equal to that delivered by 
the ascending belts for a wide range of current 
values. This condition is realized by placing in 
parallel with the belt spray combs a voltage 
control comb of the same design to deposit 
charge directly on the metal frame. It is desirable 
to make the proper choice of control comb length, 
corona gap, and series resistance, though these 
factors are not critical. 

The lengths of the control comb and of the 
combined belt spray combs should be in pro- 
portion to their relative currents. These currents 
would be equal but for the fact that the collector 
removes only about two-thirds of the charge 
from the ascending belt, the remainder being 
neutralized by the belt spray combs. Since the 
control comb should have approximately the 


same current-voltage characteristic as the com- 








768 VAN ATTA, NORTHRUP, VAN 
bined belt spray combs, its corona gap should be 
the same as theirs and its series resistance should 
equal their series resistances combined in parallel. 

These conditions are expressed more precisely 


by the following relations 


(R,+7,)/(Rs+7r,) =1./I,~2.0 


and 


R,/t.=R,/To; 


where R,, 7,, and J, 
resistance, the internal nonlinear resistance, and 
the corona current, respectively, of the voltage 
control corona gap, and R,, r,, and J, refer 
similarly to the combined spray corona gaps. 
The second relation is necessary in order that the 
first relation may hold for a wide range of values 
of frame voltage. The application of these con- 
ditions can be made more readily understandable 


are the external series 


by reference to Fig. 8, in which the current 
network existing in the upper charging system is 
shown schematically. 


BELT-CHARGING SYSTEM: SINGLE 


EXCITATION 


MODIFIED 


The alterations indicated in the preceding 
section have been incorporated in the _ belt- 
charging system described in this section and 
shown schematically in Fig. 9a. Point spray 
combs have replaced the spray wires and col- 
lector wires. The belt spray combs are used with 
}-inch corona gaps and 7-megohm series resistors. 
The collector combs are placed within } inch of 
the belt, this being the closest safe distance for 
mechanical reasons. The voltage control comb 
has half the length of the combined upper spray 
combs, about the same corona gap, and half the 
series resistance. In practice minor adjustments 
in the characteristic of the control comb are 
made by altering its corona gap. Let us suppose, 
for example, that saturation charge density on 
the ascending belts, as indicated by discharge at 
the point where the belts leave the lower pulleys, 
is not accompanied by saturation on the de- 
scending belts. The corona gap of the control 
increased, with a resultant 


comb should be 


increase in upper frame voltage and current on 
the descending belts, until discharge appears 
where the belts leave the upper pulleys. Because 
of the similarity of the characteristics of the 
control comb and the belt spray combs, the 
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Schematic diagram of the upper belt-charging 
system for single excitation. 


Fic. 8. 


charge densities on the ascending and descending 
belts will then be approximately equal for all 
values of charge density. 

The upper and lower frames are connected to 
the terminal and to ground, respectively, through 
200-megohm resistors to prevent the generator 
from operating self-exciting because of the slight 
frictional charge generated on the paper belts. 
These resistors together with those in series with 
the spray combs resistors of 
approximately linear response developed in this 
laboratory for use at high voltage, and will be 
described in a later publication. The voltage 
applied to the lower frame is continuously 
variable from 0 to 20 kv by means of a small 
adjustable auto-transformer® in the primary of 
the power transformer of the a.c. power supply. 

In order to write down in compact form the 
current relations involved in the operation of this 
belt-charging system, the following symbols are 


are carbon-line 


defined. 


I¢=total charging current to the terminal 
74=current carried by ascending belts 

7 p=current carried by descending belts 

’¢ =corona current at upper collectors 


* Variac, General Radio Co., Cambridge, Mass. 
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Fic. 9. Schematic diagram of the modified belt-charging 
systems; a, single excitation; b, double excitation. 


Us=corona current at upper spray combs 

Uy =corona current at upper voltage control comb 
U,=leakage current through belts at upper pullevs 

Ur =current through grounding resistor on upper frame 
Ls=corona current at lower spray combs 

L,=leakage current through belts at lower pulleys 


L 


The following exact current relations are then 
obvious from an examination of Figs. 8 and 9a. 


Te=I4tIp=Ust+Uvt+U 
UstUc—UL=Ls—-Lt. 


Since in practice Up is always relatively small 
and U, and L, are small except at saturation 
charge density, the above relations may be 
written, approximately, 


Ic=I4t+lIp= Ust Uy= UstUc=Ls. 


These relations give a concise description of the 
operation of the charging system and are useful 
in connection with initial design and adjustments. 

The belt-charging system described in this 
section was used for the major part of the current 
and voltage measurements to be described. It 
has proved capable of giving very steady currents 
over a wide range of current values. This system 
will probably be used here considerably in the 
future, in spite of certain advantages offered by 
the belt-charging system which employs double 
excitation. 


MopiIFIED BELT-CHARGING SYSTEM: DOUBLE 
EXCITATION 


The necessity for obtaining excitation at the 
upper charging system from charge on the 
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ascending belt introduces certain complications 
and adjustments which can be avoided if an 
independent voltage source for excitation is 
employed in the terminal. The resulting charging 
system is shown diagrammatically in Fig. 9b. 

It is unnecessary that the upper and lower 
exciting voltages have the same value except at 
small charge densities, when they must be kept 
the same within about ten percent in order to 
realize steady currents. Thus the current can be 
adjusted over rather wide limits by varying 
either exciting voltage independently. For most 
convenient operation, however, the controls on 
the two exciting voltages should be operated in 
tandem by means of a cord or insulating rod 
connection. 

In Fig. 10 several comparative charging cur- 
rent curves are given. These curves show that 
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Fic. 10. Comparative charging current curves; a, single 
and double excitation with negative polarity and 4520 
f.p.m. belt speed; b, positive and negative polarity with 
5650 f.p.m. belt speed and double excitation; c, 4520 f.p.m. 
and 5650 f.p.m. belt speeds with negative polarity and 
double excitation. 
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greater currents may be obtained with double 


excitation than with single excitation, that 
greater currents may be delivered to a negative 
terminal than to a positive terminal, and that the 
current increases with belt speed, though not 
quite proportionally under the conditions ob- 
taining. In addition, the amplitude of fluctuations 
in the charging current with double excitation is 
about half that obtained with single excitation. 

The maximum charge density on the belts 
realized experimentally is 30 percent of the 
theoretical limit which corresponds to a charge 
density of 2.65 X 10~-* coulomb per sq. cm on each 
surface of the belt. Actually, since the charge is 
practically limited to one side of the belt at the 
time that the belt leaves the pulley, this charging 
current represents 60 percent of the theoretical 
limit available with the present charging 
arrangement. 

The results described above indicate that a 
charging system employing double excitation is 
to be preferred when the second voltage supply 


can be obtained without inconvenience. 


VOLTAGE MEASUREMENT 


The problem of obtaining a convenient and 
reliable measurement of the high voltage has 
received considerable attention. Several possible 
methods of voltage measurement have suggested 
resistance, attracted disk, 
generating magnetic deflection of 
ions, and proton range. The high 
method requires a resistance difficult to con- 
struct, difficult to shield from corona, and almost 


high 
voltmeter, 


themselves: 


resistance 


necessarily nonlinear because of the dimensions 


—_ 
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of the apparatus and the magnitude of the 
voltage involved. An attracted disk voltmeter is 
subject to large calibrating errors because of the 
squared scale and the relatively low available 
calibrating voltage. The magnetic deflection and 
proton range methods are more complicated than 
is desirable for many purposes, but may serve to 
confirm and standardize the results of some more 
direct method. 

The generating voltmeter method has been 
used with good results for a variety of purposes.'” 
Conflicting results have been obtained, however, 
concerning the reliability of generating voltmeter 
measurements made on belt electrostatic gener- 
ators.'! Because of the advantages offered by its 
relative simplicity and its independence of the 
ion beam, the method has been experimentally 
tested in this laboratory with special attention to 
systematic errors. 

The generating voltmeter consists essentially 
of an electrostatic alternator excited by the field 
associated with the voltage to be measured, 
together with an amplifier and an output meter. 
In the alternator shown in Fig. 11, a two-bladed 
grounded rotor plate is mounted in front of a 
grounded stator plate on a shaft extending 
through a small hole in the stator plate to a 
synchronous motor. Two diametrically opposite 


10 P, Kirkpatrick and I. Mayake, Rev. Sci. Inst. 3, 1 
(1932); Ross Gunn, Phys. Rev. 40, 307 (1932); P. Kirk- 
patrick, Rev. Sci. Inst. 3, 430 (1932); G. P. Harnwell and 
S. N. Van Voorhis, Rev. Sci. Inst. 4, 540 (1933); J. E 
Henderson, W. H. Goss and J. E. Rose, Rev. Sci. Inst. 6, 
63 (1935). 

uM. A. Tuve and L. R. Hafstad, Phys. Rev. 45, 651 
(1934); M. A. Tuve, L. R. Hafstad and O. Dahl, Phys 
Rev. 48, 315 (1935); R. G. Herb, D. B. Parkinson and 
D. W. Kerst, Rev. Sci. Inst. 6, 261 (1935). 
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Fic. 11. Component parts of the generating voltmeter. 
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quadrants are insulated from the remainder of 
the stator plate and connected to ground through 
resistors ranging from 300 to 110° ohms by 
means of an eight-point selector switch. As the 
rotor alternately shields and exposes the stator 
quadrants in the presence of an electrostatic 
field, an alternating current is produced in the 
grounding resistor. 

In order to avoid harmonics of the fundamental 
frequency and to simplify calculations, the rotor 
blades are so shaped that the generated current 
is sinusoidal. This reduces the effective area of 
the rotor by a factor of 7/2. Since the alternator 
arrangement produces an electrical frequency 
twice the mechanical frequency, the rotor speed 
of 3600 r.p.m. produces a 120 c.p.s. alternating 
current. 

The peak value of current delivered by the 
alternator may be written 

]=}FAS statamperes 
or ]=2rCl\'S statamperes 
where F=electrostatic field on the voltmeter (e.s. volts 
cm) 
A =effective area of the rotor (cm?) 
S=rotor speed (r.p.s.) 
‘=capacity between area A on the quadrants and 


the voltage terminal (cm) 
V = voltage applied to the terminal (e.s. volts) 


The former expression has proved more useful 
with this type of alternator, while the latter is 
equivalent to that used by Kirkpatrick. For 
convenient use the root-mean-square voltage 
across the resistor is given by the expression 


E=3.93XK10°' FASR millivolts 


where R = grounding resistance for the quadrants 
(ohms) and F is measured in volts/cm. This 
expression has been verified experimentally. 

The voltage generated across the resistance by 
the alternator is applied to the input of a three- 
stage, resistance-capacity coupled amplifier em- 
ploying type 41 tubes operated at reduced 
voltages from self-contained B and C batteries. 
The impedance of the output stage is matched to 
that of a 10-ma thermocouple meter providing a 
measured over-all transconductance of 14.0 
milliamperes per millivolt at 120 c.p.s. A full- 
scale reading on the output meter is thus 
obtained with a field strength on the voltmeter 


ranging from 0.5 volt/em to 2000 volts/cm 
depending upon the input resistance employed. 

The aiternator and amplifier described above 
have proved quite useful in the early phase of the 
voltage measurements because of the sensitivity 
and wide range afforded. Experience with this 
voltmeter and earlier models has led to the design 
of a much simplified voltmeter which will be 
described in some detail in a future publication. 

In order to obtain a measure of the linearity of 
the voltmeter response, known voltages were 
applied to an insulated metal plate temporarily 
mounted at a known distance in front of the 
voltmeter surface. The output current from the 
amplifier proved to be accurately proportional to 
the gradient over the entire range within which 
measurements are to be taken. This calibration 
served to determine the constant in the expression 


F=k,-I/R 


for the gradient on the voltmeter in any location, 
where J is the current in the output meter and 
R is the input resistance. 

For a particular location of the voltmeter and 
high voltage terminal a known calibrating volt- 
age, separately measured by means of a precision 
high resistance™ and milliammeter, is employed 
to determine the constant in the expression 


V=hke-I/R. 


Serious errors in voltage measurement, espe- 
cially at calibrating voltages, may result from 
the large effect on the voltmeter of charged 
insulating surfaces. The first step in avoiding 
these errors is to locate the voltmeter as favorably 
as possible. Even with the voltmeter mounted in 
the top of the terminal opposite the grounded 
metal roof, the effect of charge on the surface of 
the supporting column must be taken into 
account at calibrating voltages to obtain an 
accurate calibration constant. If the current 
output from the generating voltmeter, Jo, due to 
column charge, is measured with the terminal 
grounded, it may be taken into consideration by 
modifying the above expression as follows 


V=ke-(I4]s)/R, 


the sign being determined by the sign of the 


? Tavlor high voltage resistor, Shallcross Mfg. Co., 
Collingdale, Pa. 
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charge on the column. With this correction it is 
possible to obtain results with a probable error of 
less than one percent for individual readings. 
When no correction is made for column charge, 
the apparent value of ke varies with calibrating 
voltage as is illustrated by the curves of Fig. 12, 
taken at low calibrating voltages for two different 
positions of the generating voltmeter at the top 
port hole of the high voltage terminal. An 
examination of these curves in the light of the 
expression just above indicates that in the case 
of the upper pair of curves the column charge was 
of the same polarity as the terminal, whereas in 
the case of the lower pair of curves the column 
charge was negative for both This 
difference is correlated with the past history of 


curves. 


the generator in the two cases. 

The voltmeter has been used also, high on the 
side wall of the building, directly opposite the 
terminal. Whereas this position might appear to 
be quite favorable, the column effect is con- 
siderable. For example, with 25 kilovolts applied 
to the terminal the uncorrected reading of the 
voltmeter may gradually increase from 25 kilo- 
volts to 50 kilovolts over the period of an hour. 
Since the change is very slow it is still possible to 
get accurate results by taking the value of J) by 
grounding the terminal after each voltage read wag. 
The difficulty is that this correction must be 
made even at voltages as high as one megavolt 
with the voltmeter in this relatively unfavorable 
location. 
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Fic. 12. Uncorrected calibration curves for the generating 
voltmeter, obtained with low calibrating voltages. 
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It should be mentioned that these corrections 
for the effect of charged insulators depend 
critically upon the humidity and may become 
quite negligible in damp weather. 

The presence about the high voltage terminal 
of a space charge of ions produced by corona has 
been suggested" as a possible source of system- 
atic error in generating voltmeter measurements. 
With certain simplifying assumptions, the ap- 
proximate magnitude of this effect may be 
readily calculated. If we assume as one of the two 
extreme cases that all the corona is produced at 
the surface of the terminal, and therefore that 
the space charge is all of the same sign, we can 
calculate the maximum space charge effect. 

The geometry is simplified by assuming a 
spherical terminal in a spherical metal room. 
Then the current through a spherical surface 
about the terminal is 


i=dvA 


where 7=corona current between terminal and ground 


(amperes), 


~ 8 


=density of charge (coulombs/cc), 
v=velocity of ions (cm/sec.), 


A =area of spherical surface (sq. cm). 
Then, if we make the substitutions 
v= FM=Fyr?M/r’ 
and A=4rr’, 
we obtain the expression for the density at any 
point in the intervening space : 
d=1/4rM For,’, 


where M=mobility of ions in air =1.6 cm/sec./volt /em 
F,=gradient at the surface of the terminal (volts 
cm), 
ro =radius of the terminal = 230 cm. 


At the maximum generator voltage 


F,~20,000 volts/em, i~5X10~‘ ampere, 


so that d=2.3X10-" coulomb/cc. 


Then the total charge in the intervening space 
between the terminal and the wall of the room, 
r,=800 cm, is 


QO, = (4/3) x(r3—1°)d =4.9 X 10-5 coulomb. 


1M. A. Tuve, L. R. Hafstad and O. Dahl, Phys. Rev. 
48, 315 (1935). 
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For the sake of comparing this total quantity 
of space charge with the charge, Qe, on the 
surface of the terminal with regard to their 
effects on the voltmeter in the wall of the room, 
both charges may be considered to be concen- 
trated at the center of the sphere. For the 
conditions assumed above, 


Q2= For? X1.1K10-" = 1.16 K 10? coulomb, 


so that Q:/Q2=4.2 percent. This effect would be 
smaller on a smaller generator in proportion to 
the linear dimensions of the apparatus. With a 
given apparatus the voltmeter should be as near 
the terminal as possible, since the effect depends 
upon 7;°. 

As an indication of the relative importance of 
space charge to surface charge on insulators, the 
surface charge density corresponding to a field of 
20,000 volts/cm is 1.8X10-*° coulomb/sq. cm. 
Since the total space charge, Q:, is equivalent to 
this surface charge on 27 sq. ft. of insulation and 
since the surface exposed by one supporting 
column is 400 sq. ft., space charge might be 
expected to have a considerably smaller effect 
than insulator charge. 

An experimental search for space charge is 
possible since the effects on the gradients at the 
terminal and at the wall of the room would be of 
opposite sign. With a uniform space charge of the 
same polarity as the terminal, a voltmeter at 
ground potential would read too high, while a 
voltmeter in the terminal would read too low. 
Measurements based on these considerations 
have been made to ascertain whether space 
charge produces an appreciable effect under the 
conditions existing in this laboratory. 

The curves of Fig: 13 represent the positive 
terminal voltage as a function of charging current. 
For one curve the voltmeter was located at the 
top of the high voltage terminal as indicated in 
Fig. 2, while for the other the voltmeter was 
mounted in the side port hole of the grounded 
terminal facing the high voltage terminal at a 
distance of 25 feet. Since one of these curves was 
taken several weeks before the other, the small 
difference between them may be explained as due 
to a slight change in the charging current 
necessary to realize the corresponding voltage. 
Differences as great as these have been obtained 
from day to day between curves taken under 
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Fic. 13. Voltage-current characteristic curves for the 
positive terminal measured with the voltmeter in two dif- 
ferent locations. 


apparently identical conditions. In other pairs of 
curves similar to that shown in Fig. 13 the 
relative positions of the two curves have been 
reversed. The most significant part of the curves 
is the upper end, where an increasing departure 
would be evident if there were an appreciable 
error due to space charge. 

These experimental results indicate that the 
space charge effect does not approach the ex- 
treme value calculated above. This is not sur- 
prising in view of the inherent uncertainty in the 
calculations introduced by the simplifying as- 
sumptions. In the first place it has been actually 
observed that, in addition to the corona produced 
at the surface of the terminal, corona of the 
opposite sign originates at the hangar wall 
resulting in a partial cancellation of space charge 
effects. In the second place the interior surface of 
the hangar is very irregular and far from 
spherical. 

Finally, generating voltmeter measurements 
were compared directly with sphere gap measure- 
ments up to 1 megavolt using 62-cm spheres. 
Typical results are shown in Fig. 14. The points 
represent the sparking voltages as measured 
with the generating voltmeter. The upper curve 
was calculated from Peek’s formula,'* but the 
lower curve, which is two percent lower than the 


4 F, W. Peek, Jr., Trans. A.1.E.E. 33, 923 (1914). 
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Fic. 14. Sphere gap voltage measurements with 62-cm 
diameter spheres. 


other throughout, gives the best fit for the 
experimental points to be obtained with a curve 
of this form. In view of the fact that more recent 
experimental work'® has shown that Peek’s 
formula gives results consistently high, it appears 
that the absolute magnitude of the generating 
voltmeter results has been confirmed within an 
estimated probable error of two percent. 

The considerations outlined in this section 
indicate that the generating voltmeter method 
may be made reliable and reasonably accurate, as 
well as convenient. A final comparison with 
proton range and deflection measurements will 
be made at the earliest opportunity. 


CURRENT AND VOLTAGE FLUCTUATIONS 


Preliminary current and voltage fluctuation 
measurements have been made to determine the 
usefulness of the generator for the study of 
resonance disintegration effects. Current fluctu- 
ations were measured with the terminal at high 
voltage by means of a cathode-ray oscillograph" 
across a resistor in the connection between the 
upper charging assembly and the terminal. Cur- 
rent fluctuations with the terminal grounded 
were made with the oscillograph across a resistor 
in the terminal grounding lead. 

A visual examination showed that the principal 


46 J. R. Meador, Elec. Eng. 53, 942 (1934); J. E. Hender- 
son, W. H. Gossand J. E. Rose, Rev. Sci. Inst. 6, 63 (1935). 
1% RCA Model TMV-122-B. 
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component of the current fluctuations may be 
represented with sufficient accuracy for the 
present purpose by a single loop of a 12 c.p.s. 
sinusoidal wave appearing once for each revolu- 
tion of the belt. It was found to be caused by a 
slight dip in the charging current at the time that 
the belt joint was passing over the upper pulley. 
Since the resulting regularly-spaced pulses were 
much larger than the irregular background upon 
which they were imposed, the percentage current 
Huctuation is defined in terms of the ratio of this 
a.c. component peak value to the d.c. component 
The results for singly-excited and doubly-excited 
charging systems are given in Fig. 15. 

Since the terminal voltage represents an equi- 
librium between charging current and _ load 
current, fluctuations in this voltage may arise 
from fluctuations in either current. The load 
current consists of insulator leakage principally 
down the supporting column, ion current in the 
accelerating tube, and corona current. Except at 
the highest voltages, the ion current should be 
the principal source of load current fluctuations. 

If we neglect for the moment load current 
Huctuations, the percentage fluctuation in termi- 
nal voltage is considerably smaller than the 
associated fluctuation in charging current be- 
cause of the smoothing effect of the terminal-to- 
ground capacity. To the extent that the a.c. 
component of the current can be represented by a 
sinusoidal wave form of frequency w, the 
smoothing effect can be calculated from the 
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Fic. 15. Fluctuations in charging current with single excita 
tion and double excitation. 
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eXpression 


voltage fluctuation (percent 
S= 
current fluctuation (percent) 


| a 1 
r Rac (Ruew®C?+1)! 


where R,., Rac, and C are, respectively, the 
terminal-to-ground a.c. resistance, d.c. resistance, 
and capacity. R,. and Ra. for any given voltage 
may be determined from such curves as those 
shown in Fig. 13. For example, if we assume an 
effective frequency of 12 c.p.s. and a terminal 
capacity of 300 micromicrofarads and take R,, 
and R,. from Fig. 13, the smoothing factor, S, at 
2 megavolts is 0.01 for no ion current and 0.03 
for one ma of ion current, and at 1 megavolt the 
factor is 0.003 for no ion current and 0.06 for one 
ma of ion current. In view of the assumptions 
involved in these calculations, the results are 
useful only in showing the trend of the smoothing 
effect and in indicating that this effect is 
considerable. 

A direct measurement of voltage fluctuation 
was made by exposing an insulated 15-cm hemi- 
spherical electrode to ground at a port hole of the 
terminal. With the electrode connected to the 
terminal through the input resistor of an ampli- 
fier and cathode-ray oscillograph, a small fraction 
of the terminal voltage fluctuation appeared on 
the oscillograph. The amplitude of fluctuation 
measured in this way was consistently much 
smaller than 0.1 percent, but exact values are not 
given because the calculation of the percentage 
of the a.c. terminal voltage on the insulated 
electrode involved uncertainties. 

These preliminary current and voltage fluctua- 
tion measurements have sufficed to show that the 
limitation in precision of voltage setting imposed 
by voltage fluctuations is small compared with 
other limitations. This point will be confirmed by 
an improved method for measuring voltage 
fluctuations at the time that the disintegration 
experiments are performed. 

Voltage variations requiring a time greater 
than a few seconds may be distinguished arbi- 
trarily from the fluctuations discussed above 
because of the fact that they can be considerably 
reduced by automatic or manual control of the 
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Fic. 16. Generator performance curves for the positive and 
negative terminals. 


excitation voltage. These variations are not 
large. By the use of a double Variac arrangement 
for manual control, a given terminal voltage can 
be maintained constant with a probable error of 
0.1 percent. 


(GENERATOR PERFORMANCE DATA 


The maximum measured voltage of a single 
generator unit, obtained with a charging current 
of 1.5 ma, was 2.72 megavolts positive and 3.15 
megavolts negative. The highest voltage which 
has been obtained consistently without sparking 
is 2.4 megavolts positive and 2.7 megavolts 
negative, giving a total voltage of 5.1 megavolts 
between the generator terminals. At this voltage 
there is still 1.1 ma available for use with an 
accelerating tube. The humidity control inside 
the generator column makes it possible to obtain 
this voltage with an external humidity of 100 
percent. 

A typical set of generator performance curves 
for the positive and negative terminals is given in 
Fig. 16. The full-load charging current, J, in Fig. 
9a, was measured as a function of voltage for the 
positive terminal with saturation charge density 
on the belts and the terminal voltage limited by a 
grounded corona comb mounted in the opposite 
terminal. Because of the persistent sparking to 
the corona comb limiting the negative voltage, 
the full-load charging current for the negative 
terminal was measured at zero and full voltage 
only. For both curves the charging current 
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decreases by only 4 percent between zero and full 
voltage. 

The no-load charging current curves of Fig. 16 
were obtained by varying the charge density on 
the belts and measuring the resulting voltage of 
the unloaded terminal. Below one megavolt the 
loss current is largely leakage down the column, 
as measured by a the base of the 


column, while above onemegavolt loss by corona 


meter at 


predominates. 

The difference between the full-load and the 
no-load charging current represents the available 
output current. Reference to the figure indicates 
that more than half of the full-load charging 
current is available at the maximum generator 
voltage. The reliability of the output current 
curves obtained in this way has been confirmed 
for the positive terminal by direct measurement 
of the current to an external corona comb used 
for voltage limitation. 

The power output curves in the figure represent 
the product of the output current and the 
corresponding voltage. The power available from 
the generator increases with voltage to a maxi- 
mum of 6.3 kw at a total voltage of 4.4 megavolts. 
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On the Absorption of Slow Neutrons in Boron 


Recently attempts have been made!:*? to measure the 
energies of neutrons selectively absorbed in elements like 
Ag, Rh, I. This determination is based on the assumption 
that, in the case of light elements which are disintegrated 
with emission of a charged particle, no sharp resonance 
occurs, and therefore the capture cross section follows the 
1/v law.*: 4 

It seemed worth while to obtain a direct check of the 
validity of this law, at least for neutrons of thermal 
energies, and for this purpose we performed with boron an 
experiment already used with Ag and Cd for determining 
the dependence of the capture cross section upon the 
neutron velocity.6 An aluminum disk of 50 cm diameter 
was coated with boron carbide, of such thickness as to 
transmit about 40 percent of the slow neutrons under the 
experimental conditions, and was rotated at 5500 r.p.m., 
all conditions being substantially the same as in the former 
experiments.® By rotating the disk first in one direction 
and then the other, we observed a change in apparent 
absorption coefficient of 1.3 percent, with a statistical 
error of 0.8 percent. This change is so small that it may well 
be due to experimental error, and we may conclude that 
our result is not inconsistent with the assumption that 
there is no change in transmitted intensity depending 
upon the velocity of the boron absorber. This result, 
distinctly different from that obtained with cadmium, 
where a change of as much as 8 percent was found, ap- 
pears to indicate that the absorption of thermal energy 
neutrons in boron does not depart appreciably from the 
1/v law. 

Doubt might arise whether this conclusion is still 
justified in a case like the present one, where on account 
of their zero-point energy, the random motions of the 
boron atoms have an average velocity which is higher than 
the velocity of the rotating disk. Actually, on account of 
the vibrations in the crystal lattice, the velocities of the 
boron atoms will roughly correspond to those of free 
particles of the same mass at the Debye characteristi 
temperature of the lattice, which we may assume to be 
about five times higher than room temperature. However, 
it can be shown that, even with high velocities of the boron 
atoms, the conclusion is still valid that, under the assump 
tion of the 1/v law, the transmitted intensity cannot be 
appreciably affected by the rotation of the disk. 

F. RASETTI 
D. P. MITCHELI 
Pupin Physics Laboratories, G. A. FInk 


Columbia University, ‘ 
April 15, 1936. G. B. PEGRAM 


1 Weekes, Livingston and Bethe, Phys. Rev. 49, 471 (1936) 
? Rasetti, Fink, Goldsmit! 1 
Washington Meeting, April 30, 1936 
3’ Wigner and Breit, Phys. Rev. 49, 519 (1936). 
* Bohr, Nature 137, 344 (1936 
Rasetti, Segré, Fink, Dunning and Pegram, Phys. Rev. 49, 104 
(1936) 
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The Infrared Spectra of H.S, HDS, and D.S 


Nielsen and Nielsen have just. published a preliminary 
account! of an investigation of the above spectra. We have 
been occupied with this subject for some time, and it is 
interesting to contrast the two sets of results, since we 
find some measure of disagreement between them. It 
must be confessed, however, that the peculiarities of the 
spectra coupled with their weakness leave room for doubt 
in their interpretation. 

The existence of the H.S band at 8x is confirmed, and 
for the three molecules we find we = 1236, 1078, and 898 as 
compared with 1266, 1111, and 926 cm™ (N. and N 
The band in both H.2S and D.S is peculiar in exhibiting no 
resolution into P and R branches, although possessing 
rotational structure; we do find for HDS that there is a 
strong peak at 998 cm™ which is apparently unobserved 
by Nielsen and Nielsen. If this is the P branch of this long 
wave band, we must regard it as missing in the other two 
molec ules. 

w; for H.S is given by Nielsen and Nielsen at 2703 cm 
We consider a more likely value for the center is 2632 cm 
In D.S these authors find a band at 2198 cm which they 
attribute to the same mode: This frequency is much toe 
high, and without correcting for anharmonicity gives a 
force constant of some 5.310° dynes per cm, whereas 
Badger’s rule makes it certain that it should be of the 
order of 4.2 10° dynes per cm. We find a band at 1940 
cm which fulfills this requirement, and conforms with the 
Raman spectra of H2S and D.S. 

We have made an estimate of the anharmonic coefficients 


and using the mechanical frequencies corresponding t 
infinitesimal amplitude, find in the notation of Van Vleck 
and Cross,” that ki, Ry, ks, and ky are, respectively, 4.25, 
0.07, 0.24 and 0.04 10° dynes per cm for H.S; and 4.34, 
0.07, 0.23, and 0.04 for D.S. These values are for a vertical 


angle of 110° which fits the admittedly crude deductions 


made from the vibrational spectrum, as compared with 
92° 20’ obtained by Cross* from the rotational analysis of 
the band at 9910 cm~. The difference between the values 
of k; for H.S and D.S is not unlikely, since it is only the 
energy of dissociation, D, reckoned from the level (0,0,0) 
which we may expect to be the same to a high degree of 


approximation in polyatomic molecules and their isotopes. 


C. R. BatLey 
J. W. THompson 
J. B. Haut 


The Sir William Ramsay Laboratories 
of Inorganic and Physical Chemistry, 
University College, London, 
April 8, 1936 


1A. H. Nielsen and H. H. Nielsen, J. Chem. Phys. 4, 229 (1936 
J. H. Van Vleck and P. C. Cross, J. Chem. Phys. 1, 357 (1933 
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The Disintegration of Nitrogen by Slow Neutrons 


The disintegration of nitrogen by slow neutrons was 
Chadwick and Goldhaber! 


and by ourselves? to the reaction 


attributed independently by 


7N'4+-on'—,B"+.Het+Q). 1) 


However, the value of Q; obtained by Chadwick and Gold- 
haber by an ionization method was 0.5 MEV, while we 
deduced the value 2.33 MEV from our range measure- 
ments of tracks observed in a cloud chamber. Our value 
was calculated from the 1.06-cm range of the disintegration 
particles by assuming that the disintegration particles 
were B" nuclei and alpha-particles. In a private communi- 
cation Dr. Goldhaber has suggested that this apparent 
disagreement can be eliminated if we interpret the disin- 


tegration in a different way as postulated in the reaction 


>N4+ on'—>,C4+ ,H'+Qz. 2) 

The energy of the disintegration Q2 calculated from our 
range measurements then 0.58+0.03 MEV 
which agrees with Chadwick and Goldhaber’s ionization 
value. We have now found another reason why the former 
interpretation is impossible. We have just measured the 


becomes 


energy released in the reaction 

sB' + ,H*—>,C®+ on'+Q; (3) 
and found Q; to be 13.540.3 MEV.* Combining this result 
with that of Cockcroft and Lewis‘ in their investigation 
of the reaction 

7N"+ ,H?+.C”+ .Het+Q, (4) 
one obtains Q; = Q4—Q;= —0.28+0.3 MEV. 

This is quite different from the value 1.5 MEV which 
was calculated from Bethe’s masses and which appeared 
to agree fairly well with the value 2.33 MEV. However, 
our value of Q;= —0.28 MEV agrees quite well with the 
respective values Q;=—0.3 and Q,;=—0.2, calculated 
from the recent list of masses given by Oliphant® and 
by Cockcroft and Lewis.‘ This new value of Q; necessarily 
eliminates (1) from further consideration in 
regard to the 1.06 cm particles resulting from the disin- 


reaction 


tegration of nitrogen by slow neutrons. 

Recently McMillan® has obtained evidence of a radio- 
active C™ with a half-life of about 3 months and a maxi- 
mum beta-ray energy of about 0.2 MEY. From reactions 
(2) and (5) 

C+, N4+e-+0,; 5) 
we see that om'—,H'=Q.+Q,;; hence from these data the 
mass difference between a neutron and a hydrogen atom 
appears to be 0.8 MEV. This is in good agreement with 
Aston’s new masses of ,H! and ,H?, and Feather's binding 
2.26 MEV). 

T. W. 
W. M. 


energy of the deuteron 
BONNER* 
BRUBAKER 
W. K. Kellogg Radiation Laboratory, 
California Institute of Technology, 
April 22, 1936. 
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* Cockcroft and Lewis, Proc. Roy. Soc. A154, 261 
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Isotopic Constitution of Strontium, Barium, and Indium 


Mass spectrograph analyses have been made of the ions 
of strontium, barium, and indium, emitted when oxides of 
these elements are heated on tungsten filaments. Evidence 
has been obtained for two new isotopes, one of strontium 
and one of barium. No new indium isotopes were found 
The observations reported below have been checked on 
two different instruments. 

In the case of strontium, a peak was observed at mass 
84 whose height corresponded to 0.5 percent of the total 
strontium emission. The masses of the known isotopes of 
strontium are all greater than 84 so the possibility is 
excluded that the effect is due to a strontium compound. 
Since the only known isotope of mass 84 belongs to krypton 
and the other krypton isotopes do not appear, we believe 
that this peak is due to a new isotope of strontium. A small 
impurity of rubidium was present initially but disappeared 
with continued heating of the filament. The relative 
abundance of the known strontium isotopes was found to 
be 9.6 percent of mass 86, 7.5 percent of mass 87, and 82.4 
percent of mass 88, differing slightly from the rough values 
given by Aston.' Assuming a packing fraction of —8.2 this 
gives 87.62 for the atomic weight of strontium, in good 
agreement with the chemical value of 87.63. A search fi 
other isotopes revealed that masses 82, 83, 85, 89, and 96 
were present to less than 1 in 2000, 1 in 1000, 1 in 2000, 
1 in 500 and 1 in 2000, respectively. 

The curves for barium showed a peak at mass 134 making 
up 1.8 percent of the total emission. The same considera- 
tions which applied in the case of strontium lead us to 
believe that this is due to a new isotope of barium. Accept- 
ing Aston’s relative abundance for the known isotopes 135, 
136, 137, and 138, and assuming packing fractions of —6.1 
we deduce an atomic weight for barium of 137.36 in exact 
agreement with the international atomic weight. Barium 
isotopes of masses 132, 133, 139, and 140 were shown to be 
present to less than 1 in 2000, 1 in 1000, 1 in 400, and 1 in 
1000, respectively. 

Evidence for the possible existence of a third isotope of 
indium is given by the observation? of three artificially 
radioactive periods induced in indium by neutron bom- 
bardment. However, only two of these periods are water 
sensitive. Mass spectrograph analyses revealed only the 
two known isotopes. A search over the mass range from 
110 to 119 showed that masses 110, 111, 112, 114, 116, 117, 
118, and 119 were present to less than 1 in 5000, 1 in 10,000, 
1 in 5000, 1 in 200, 1 in 5000, 1 in 10,000, 1 in 8000, and 1 
in 30,000, respectively. The ratio of the abundance of mass 
115 to that of mass 113 was found to be 21+1 in agreement 
with Aston’s observation.’ If we assume a packing fraction 
of —6 we obtain for indium an atomic weight of 114.81 
in agreement with Aston’s value. The international atomic 
weight for indium is 114.76. 

J. P. BLEwetT1 
M. B. SAMPSON 


Palmer Physical Laboratory, 
Princeton University, 
April 28, 1936. 
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